Научная статья на тему 'ФИЛЬТРАЦИЯ СИГНАЛОВ ПРИ СКАЧКООБРАЗНЫХ МУЛЬТИПЛИКАТИВНЫХ ПОМЕХАХ'

ФИЛЬТРАЦИЯ СИГНАЛОВ ПРИ СКАЧКООБРАЗНЫХ МУЛЬТИПЛИКАТИВНЫХ ПОМЕХАХ Текст научной статьи по специальности «Физика»

CC BY
29
6
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ИНФОРМАЦИОННАЯ СИСТЕМА / СЛУЧАЙНАЯ СКАЧКООБРАЗНАЯ СТРУКТУРА / МУЛЬТИПЛИКАТИВНАЯ СКАЧКООБРАЗНАЯ ПОМЕХА / БАЙЕСОВА ФИЛЬТРАЦИЯ / ДВУХМОМЕНТНАЯ ПАРАМЕТРИЧЕСКАЯ АППРОКСИМАЦИЯ / ИНФРАКРАСНЫЙ ПЕЛЕНГАТОР ОБЪЕКТОВ

Аннотация научной статьи по физике, автор научной работы — Бухалёв В.А., Скрынников А.А., Болдинов В.А.

Решена задача построения приближенно-оптимального байесова алгоритма фильтрации выходного сигнала линейной стохастической динамической системы. Измерена нелинейная смесь выходного сигнала и мультипликативной помехи. Помеха представляет собой непрерывнозначный случайный процесс с неизвестным законом распределения в заданных пределах [c(1), c(2)], c(1) > 0, c(2) > 0, в диапазоне частот [0, ∆ω]. Синтез алгоритма фильтрации осуществляется методом приближенно-оптимального оценивания сигналов, основанным на теории систем со случайной скачкообразной структурой и методе двухмоментной параметрической аппроксимации вероятностных распределений фазовых координат. Метод состоит в приближенной замене неизвестных плотностей вероятностей фазовых координат известными законами распределения с неизвестными математическими ожиданиями и ковариациями, определяемыми в результате решения задачи. Предлагаемый приближенно-оптимальный алгоритм фильтрации основан на замене непрерывнозначной мультипликативной помехи случайным скачкообразным процессом --- марковской цепью с двумя состояниями c(1), c(2) и равными друг другу интенсивностями переходов из одного состояния в другое q′ = 2∆ω. Условные плотности вероятности выходного сигнала при фиксированных состояниях марковской цепи c(1), c(2) аппроксимируются гамма-распределениями, зависящими от условных математических ожиданий и дисперсий сигнала x(t) при фиксированных c(1), c(2) и измерений x(t) в смеси с марковской скачкообразной помехой. Рассмотрен пример построения алгоритма оценивания дальности летательного аппарата до объекта по измерениям облученности объекта инфракрасным пеленгатором. Облученность равна отношению силы излучения к квадрату дальности. Сила излучения --- мультипликативная помеха --- случайный процесс с неизвестным непрерывным распределением в ограниченном диапазоне. Заменой этого процесса марковской цепью с двумя состояниями и аппроксимацией плотности вероятности дальности рэлеевским распределением построен приближенно-оптимальный фильтр для оценки дальности

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Бухалёв В.А., Скрынников А.А., Болдинов В.А.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

SIGNAL FILTERING AT JUMP MULTIPLICATIVE INTERFERENCE

The problem of constructing the approximately optimal Bayesian algorithm of filtering the output signal for a linear stochastic dynamical system is solved. Non-linear mixture of the output signal and of the multiplicative interference was measured. This interference is a continuous-valued random process with the unknown distribution law within the given limits of [c(1), c(2)], c(1) > 0, c(2) > 0, and in the frequency range of [0, ∆ω]. The filtering algorithm is synthesized by the approximately optimal signal estimation method based on the theory of systems with random jump structure and the method of two-time parametric approximation of the phase coordinates probability distribution. The method consists in approximate replacement of the unknown probability densities of phase coordinates by the known distribution laws with the unknown mathematical expectations and covariances determined as a result of solving the problem. The proposed approximate-optimal filtering algorithm is based on replacement of the continuous-valued multiplicative interference by a random jump process, i.e., Markov chain with two states c(1), c(2) and equal intensities of transitions from one state to another q′ = 2∆ω. Conditional probability densities of the output signal for fixed states of the Markov chain c(1), c(2) are approximated by gamma distributions that depend on mathematical expectations and variances of the x(t) signal at the fixed c(1), c(2) and measurements x(t) in a mixture with the Markov jump interference. An example of constructing an algorithm for evaluating the distance from an aircraft to the object based on object irradiance measurement by the infrared direction finder was considered. Irradiance was equal to the radiation strength ratio to the square of the distance. Radiation strength, i.e., multiplicative interference, was the random process with unknown continuous distribution in a limited range. The approximate-optimal filter was constructed for distance evaluation by replacing this process with a two-state Markov chain and approximating the distance probability density by the Rayleigh distribution

Текст научной работы на тему «ФИЛЬТРАЦИЯ СИГНАЛОВ ПРИ СКАЧКООБРАЗНЫХ МУЛЬТИПЛИКАТИВНЫХ ПОМЕХАХ»

УДК 681.5

DOI: 10.18698/1812-3368-2023-2-4-16

ФИЛЬТРАЦИЯ СИГНАЛОВ ПРИ СКАЧКООБРАЗНЫХ МУЛЬТИПЛИКАТИВНЫХ ПОМЕХАХ

В.А. Бухалёв1

A.А. Скрынников2' 3

B.А. Болдинов3

1 МНИТИ, Москва, Российская Федерация

2 ГосНИИАС, Москва, Российская Федерация

3 МАИ, Москва, Российская Федерация

vadim.bukhalev@yandex.ru

a1260@mail.ru

viktorboldinov@mail.ru

Аннотация

Решена задача построения приближенно-оптимального байесова алгоритма фильтрации выходного сигнала линейной стохастической динамической системы. Измерена нелинейная смесь выходного сигнала и мультипликативной помехи. Помеха представляет собой непрерывнозначный случайный процесс с неизвестным законом распределения в заданных пределах [с(1), с(2)], ^ > 0, 62) > 0, в диапазоне частот [0, Дш]. Синтез алгоритма фильтрации осуществляется методом приближенно-оптимального оценивания сигналов, основанным на теории систем со случайной скачкообразной структурой и методе двухмоментной параметрической аппроксимации вероятностных распределений фазовых координат. Метод состоит в приближенной замене неизвестных плотностей вероятностей фазовых координат известными законами распределения с неизвестными математическими ожиданиями и ковариациями, определяемыми в результате решения задачи. Предлагаемый приближенно-оптимальный алгоритм фильтрации основан на замене непрерывнозначной мультипликативной помехи случайным скачкообразным процессом — марковской цепью с двумя состояниями № и равными друг другу интенсивностями переходов из одного состояния в другое q' = 2Дю. Условные плотности вероятности выходного сигнала при фиксированных состояниях марковской цепи ^ c(2) аппроксимируются гамма-распределениями, зависящими от условных математических ожиданий и диспер-

Ключевые слова

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

сий сигнала x(t) при фиксированных с(1), с(2) и измерений x(t) в смеси с марковской скачкообразной помехой. Рассмотрен пример построения алгоритма оценивания дальности летательного аппарата до объекта по измерениям облученности объекта инфракрасным пеленгатором. Облученность равна отношению силы излучения к квадрату дальности. Сила излучения — мультипликативная помеха — случайный процесс с неизвестным непрерывным распределением в ограниченном диапазоне. Заменой этого процесса марковской цепью с двумя состояниями и аппроксимацией плотности вероятности дальности рэлеевским распределением по- Поступила 04.06.2022 строен приближенно-оптимальный фильтр для Принята 15.02.2023 оценки дальности © Автор(ы), 2023

Работа выполнена при финансовой поддержке

РНФ (грант РНФ № 22-29-00708)

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

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

Состояние системы ССС в каждый момент времени характеризуется вектором ^(0,s(t)] (или в дискретной форме ^,sk]), где x(t) — вектор фазовых координат (переменных), принадлежащих непрерывнозначному (континуальному) множеству; s(t) — вектор состояний структуры, принадлежащих конечному множеству.

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

Другая особенность математических моделей ССС состоит в их применении в задачах фильтрации сигналов, наблюдаемых в смеси с мультипликативными скачкообразными помехами. Например, случай, когда в системе измеряется произведение или отношение двух случайных процессов ) = X ^ ^ ^) или ) = X(t)/ Y (t), из которых один является полезным сигналом, а другой — помехой. Если оба процесса континуальны, то их апостериорные оценки будут мало отличаться от априорных. Если хотя бы один из процессов описывается ССС, то, как здесь показано, с использованием методов теории систем ССС и метода двухмомент-ной параметрической аппроксимации [4, 6-8, 10-13] можно построить алгоритм фильтрации, обеспечивающий сходимость оценок.

Постановка задачи. Рассмотрим стохастическую систему:

где к — дискретный момент времени, k = 0,1,2,...; п — целое число; Х0 — случайная величина (Х0 е (0, да) с математическим ожиданием Х0 и дисперсией Я0); ак — известная функция времени, 0 < ак < 1; <^к — возмущение (последовательность независимых при различных к случайных величин с математическим ожиданием и дисперсией Ок, к е (0, да)); Хк — сигнал измерения; с* — мультипликативная помеха (случайный процесс с неизвестным непрерывным законом распределения в пределах [с(1),с(2)], с(1) > 0, с(2) > 0, в диапазоне частот [0, Аш], с-1).

Требуется найти оптимальную оценку Хк, основанную на измерениях Х0к - (Х0, Х1,..., Хк ).

Алгоритм фильтрации. Запишем этапы алгоритма.

1. Аппроксимируем непрерывнозначный марковский процесс с* марковской цепью с^к) с двумя равновероятными состояниями (с(1),с(2)). Процесс Ск описывается равными вероятностями перехода qк из состояния Sk в состояние Sk +1, Sk +1 ^ Sk, Sk, Sk +1 = 1,2: фк) = с(г) при sk = г. Найдем вероятность перехода qk, исходя из равенства частотных спектров процессов с* и с^к). Вероятность перехода связана с интенсивностью перехода qк приближенной формулой (при малом шаге счета At:

qk = q'k At, интенсивность перехода зависит от эффективной ширины спектра Аю [6-8, 10, 11, 13]):

xk+i = (1 - at) xi + ak ;

* n zk = ckxk,

(1) (2)

q'k = 2Аш.

(3)

2. Приближенно заменим в (2) помеху с* ее дискретной аппроксимацией с(вк):

гк « с(вк )хП. (4)

Обозначив ук = г1/п, Ь(вк) = [с(вк)] _1/п, преобразуем (4) к виду

Ук = Хк / Ь(вк). (5)

3. Апостериорное распределение вероятностей вектора состояния системы (Хк+1, вк+1) ]к+1(хк+1, $к+1| у ОТкТГ) определим формулой Байеса, обобщенной на класс систем ССС [6, 10]:

л _ ]к+1(хк+1, Ук+1| Вк+1, у 00~к)рк+1^к+1)

!к+1( хк+1, Вк +11 у 0ТкТ1) =--:—-, (6)

//+1( Ук+11 у

где _/к +1(хк+1, ук +1| Вк+1, у 0"^) — условное прогнозируемое на шаг вперед распределение (хк+1, ук+1) при фиксированных вк +1 и у 0 к; рк + 1(вк +1) = Р[вк +1| у 0~к] — прогнозируемая на один шаг вероятность

Вк +1; /к +1(ук +1| у — условная плотность вероятности ук +1 при фиксированных наблюдениях у на отрезке [0, к]. Из (6) следует

<х>

рк + 1(вк +1) = I¡к + 1(хк +1, вк +1 | у 0ТкТГ) йхк +1 =

0

~cy

рk + 1(sk +1)U+1(Ук + 1I Sk +1, y ~k) ~2 I :

X ~рк + 1(sk +1)//+1(Ук +11 sk +1, y 0äк

sk +1

(7)

где рк + 1(вк +1) = Р[вк +1| у 0ТкТГ] — апостериорная вероятность вк +1;

2

рк + 1(вк +1) = Е Як (вк +1| вк)рк(вк) = (1 - 2qk)рк(вк) + Як. (8)

Вк +1

4. Плотность вероятности, согласно (5), определим по формуле

йхк +1

f{+1( yk+11 sk+ь У 07k) =

dyk

+1

fk + 1(Xk +1 | Sk +1, y 0, k)

xk +1 = bk + 1(sk +1) yk +1 (9)

Аппроксимируем fk(•) плотностью вероятности гамма-распределения [13-15]:

fk (•) =

[ßk(Sk)]ak) xkiak(sk)" 1]exp[-ßk(sk)xk] при xk > 0; T(ak (sk)) k (10)

0 при xk < 0;

ak(sk) > 0, ßk(sk) > 0,

где Г(-) — гамма-функция (интеграл Эйлера), Г(а) = 11а~1в^^. При а

0

целом положительном Г(а) = (а -1)!. Подставив (10) в (9), получим

f! (■)=

^)]ak^) yk[ak(sk)"1] exp[-Yk(sk)yk] при yk > 0; T(ak (sk))

0 при ук < 0,

где у к ^к) = Ьк ^к )Рк ^к).

Условные прогнозируемые математическое ожидание у к (sk) и дисперсия Яу ^к) при фиксированном Sk гамма-распределения связаны с его параметрами ак ^к), Ук (sk) формулами [13-15]:

~ I ч «к(sk) 5у, ч ак(sk)

ук (sk) = ——; Ц ^к) = > ук ^к) у к (sk)

откуда следует

ак ^) = = ~ХМ; Кук^к) &к^к)

У к ^к) = Щ = ,

где Х к ^к), Я к (sk) — условные прогнозируемые математическое ожидание и дисперсия Хк при фиксированном Sk.

5. Прогнозируемые Х к+^к+О и Як+^к+О, согласно (1), связаны с апостериорными Хк ^к), &к ^к) на предыдущем этапе формулами [6, 10]:

X qk^к+11 Sk)рк^к)(1 - ак )Хк(sk)

Х к + l(Sk +1) = —-=-;---+ (11)

р к + l(Sk +1)

Я к+^¡к+1) = р 1\^¡к+1)Х Як (¡к+1 | ¡к ) рк (¡к ) X

¡к

х|(1 - йк )2 Як (¡к) + [ Хк+1(Вк+1) - (1 - ак )Хк (¡к) - ак!]21 + а^к, (11)

где

. . Г Як при ¡к +1 ф ¡к; Як (¡к +1| ¡к) = -!, „

[1 - 2як при ¡к +1 = ¡к.

Как следует из (5),

Хк (¡к) = Ьк (¡к)ук; Як (¡к) = 0. (12)

Подставив (12) в (11), получим

(1 - йк )X Як (¡к +1 | ¡к )рк (¡к )Ь(^к )Ук

Х к + 1^к +1) =---=-;---+ йк \к;

р к +1 (¡к +1)

Як + ^¡к +1) = Р11+ 1(^к + 1)Е Як (¡к +1 | ¡к )рк (¡к ) X (13)

¡к

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

X [ Хк+1^к+1) - (1 - йк )Ь($к)Ук - йк 1 ]2 + ак£к. Безусловная апостериорная оценка определяется по формуле

3

Хк = X Хк (¡к)рк (¡к), безусловная апостериорная дисперсия — по фор-

¡к =1

3

муле Як = X [ Як (¡к) + Х^к) ] рк (¡к) - Х2, откуда с учетом (12) следует

¡к =1

Хк = bkyk; Rk =

I b2(sk) Рк (sk) - b\

. sk

y2, (14)

где Ьк = Е Ь(5к)рк (¡к); Ук = %к'п; Ь^к) = Ш )]"Ш.

¡к

Апостериорная оценка

Ск = Х Ск (¡к )рк (¡к ). (15)

В целом формулы (7), (8), (13)-(15) образуют замкнутую систему рекуррентных уравнений, описывающую алгоритм фильтрации.

Пример. Процесс сближения летательного аппарата с объектом описывается уравнением [6, 7, 10, 12]: 1^+1 = Ьк + , где Ьк — расстояние от летательного аппарата до объекта; к — последовательность независимых при различных к случайных величин с математическим ожиданием к, £,к < 0, и дисперсией Ск.

Инфракрасный пеленгатор измеряет облученность Ек, связанную с силой излучения 1к и расстоянием Ьк формулой [16]:

Ек (16)

Ьк

где X — коэффициент, учитывающий ослабление излучения при прохождении через атмосферу, 0 < X < 1.

Сила излучения I* — случайный процесс, ограниченный по величине пределами 1(1), I(2) и лежащий в полосе частот [0, Лш], с-1.

Облученность Ек измеряется практически без ошибки Хк ~ Ек . Требуется построить алгоритм оценивания расстояния Ьк.

Алгоритм оценивания расстояния Ьк. Запишем этапы алгоритма.

1. Приближенно заменим непрерывнозначный случайный процесс 1к случайным скачкообразным процессом 1к ^к) с двумя возможными состояниями 1к ^к = 1) = I(1) и 1к ^к = 2) = I(2), где $к — марковская цепь, Sk = 1, 2, заданная равными друг другу вероятностями переходов q из одного состояния в другое. Согласно (3), вероятность переходов на одном шаге Дt: q = 2AюAt.

2. Находим апостериорную вероятность состояния структуры, согласно (7), (8):

р к+1^к+1) = (1 - 2q) рк ^к) + q; , рк + l(Sk +1)/к + 1(Хк +1| Sk +1)

рк+l(sk+1)=—2-;

X рк + 1^к +1)/к + 1(Хк +1| Sk +1) (17)

Sk +1 = 1

1

fk+1(zk+11 sk+1) = 2

Xlk + 1(sk +0 г ~ . -|

|_ Jk + 1(Lk +11 sk +1) J

Z3 , L J Lk +1 =ф( zk +1)

k +1

XIk + 1(sk +1)

где ф(гк +1) = ,

Хк +1

3. Определяем прогнозируемую и апостериорную оценки Ьк, согласно (13):

2 л - —

X q($k+1|$k)рк^к)[Ьк^к)+£,к]

Ьк+l($k+1) = —-г—-—--; (18)

р к+l($k+1)

Lo(1) = Lo(2) = Lo; Lk + i(sk +1) = .

XIk + 1(Sk +1)

(18)

Zk +1

Lk +1 = X Lk + 1(sk +1)pk + 1(sk +1).

Sk = 1

4. Аппроксимируем неизвестную плотность вероятности fk (Lk | Sk) рэлеевским распределением [13-15]:

fk (Lk | Sk) =

Lk

°k(sk)

exp

J2 > Lk

2^2 (sk)

Lk > 0,

(19)

где Gk(sk) = J—Lk(sk). V л

5. Находим из (19):

r , . . hnlk+1(sk+1) ^ fk + 1(zk +11 sk +1) = —=2—--—2—exp

4Lk + 1(sk + 1)zk +1

Xnlk+1(sk+1) 4 Lk + 1(sk + 1)zk +1

(20)

Это означает, что случайная величина Хк = 1/ Zk имеет условное экспоненциальное распределение [13-15]: fk(Хк |¡к) = ^к(¡к)ехр(-Щ;(¡кХ), где

Mk (sk) =

Д ^nlk (sk)

4 L2k (sk)

(21)

Подставив (17)-(21) в (16), получим

p k+1(2)1 (2)L2k+1(1)

pk+1(1) =

1 +

p k+1(1) I(1) Lk+1(2) pk+1(2) = 1 - pk+1(1); p k +1(1) = (1 - 2 q) pk (1) + q; p k +1(2) = 1 - p k +1(1);

po(1) = p1o; po(2) = 1 -p1o,

exp (-Ak +1)

-1

(22)

. д Xn где A k = -— 4zk

I (2) I(1)

L L2k (2) L2k (1)

Следовательно, алгоритм приближенно-оптимального оценивания расстояния Ьк по измерениям облученности гк описывается рекуррент-

ными уравнениями (18), (22). Апостериорная оценка силы излучения приближенно определяется по формуле

Ik -I(1)Pk(1) +I(2)pk(2). (23)

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

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

ЛИТЕРАТУРА

[1] Pierce B., Sworder D. Bayes and minimax controllers for a linear system with stochastic jump parameters. IEEE Trans. Automat. Contr., 1971, vol. 16, iss. 4, pp. 300-307. DOI: https://doi.org/10.1109/TAC. 1971.1099732

[2] Sworder D.D. Bayes controllers with memory for a linear system with jump parameters. IEEE Trans. Automat. Contr., 1972, vol. 17, iss. 1, pp. 119-121.

DOI: https://doi.org/10.1109/TAC. 1972.1099877

[3] Mariton M. Jump linear systems in automatic control. Taylor & Francis, 1990.

[4] Казаков И.Е., Артемьев В.М., Бухалёв В.А. Анализ систем случайной структуры. М., ФИЗМАТЛИТ, 1993.

[5] Пакшин П.В. Дискретные системы со случайными параметрами и структурой. М., ФИЗМАТЛИТ, 1994.

[6] Бухалёв В.А. Распознавание, оценивание и управление в системах со случайной скачкообразной структурой. М., Наука, 1996.

[7] Бухалёв В.А. Основы автоматики и теории управления. М., ВВИА им. Н.Е. Жуковского, 2006.

[8] Бухалёв В.А. Оптимальное сглаживание в системах со случайной скачкообразной структурой. М., ФИЗМАТЛИТ, 2013.

[9] Zhang C., Zhu H., Zhou H., et al. Deterministic and stochastic differential games. In: Non-Cooperative Stochastic Differential Game Theory of Generalized Markov Jump Linear Systems. Cham, Springer, 2017, pp. 17-29.

DOI: https://doi.org/10.1007/978-3-319-40587-2_2

[10] Бухалёв В. А., Скрынников А. А., Болдинов В. А. Алгоритмическая помехо-защита беспилотных летательных аппаратов. М., ФИЗМАТЛИТ, 2018.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

[11] Бухалёв В.А., Скрынников А.А., Болдинов В.А. Игровое управление системами со случайной скачкообразной структурой. М., ФИЗМАТЛИТ, 2021.

[12] Себряков Г.Г., Красильщиков М.Н., ред. Управление и наведение беспилотных маневренных летательных аппаратов на основе современных информационных технологий. М., ФИЗМАТЛИТ, 2005.

[13] Бухалёв В.А., Болдинов В.А., Прядкин С.П. и др. Двухмоментная параметрическая аппроксимация распределений в информационно-управляющих системах навигации и наведения. Вестник компьютерных и информационных технологий, 2016, № 8, с. 8-15. DOI: https://doi.org/10.14489/vkit.2016.08.pp.008-015

[14] Джонсон Н.Л. Одномерные непрерывные распределения. Ч. 1. М., Бином. Лаборатория знаний, 2010.

[15] Королюк В.С., ред. Справочник по теории вероятностей и математической статистике. М., Наука, 1985.

[16] Лазарев Л.П. Инфракрасные и световые приборы самонаведения и наведения летательных аппаратов. М., Машиностроение, 1976.

Бухалёв Вадим Алексеевич — д-р техн. наук, главный научный сотрудник МНИТИ (Российская Федерация, 105094, Москва, ул. Гольяновская, д. 7а, стр. 1).

Скрынников Андрей Александрович — канд. техн. наук, начальник сектора ГосНИИАС (Российская Федерация, 125319, Москва, ул. Викторенко, д. 7); доцент кафедры «Теория вероятностей и компьютерное моделирование» МАИ (Российская Федерация, 125993, Москва, Волоколамское ш., д. 4).

Болдинов Виктор Александрович — канд. техн. наук, доцент кафедры «Авиационные робототехнические системы» МАИ (Российская Федерация, 125993, Москва, Волоколамское ш., д. 4).

Просьба ссылаться на эту статью следующим образом:

Бухалёв В.А., Скрынников А.А., Болдинов В.А. Фильтрация сигналов при скачкообразных мультипликативных помехах. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки, 2023, № 2 (107), с. 4-16. DOI: https://doi.org/10.18698/1812-3368-2023-2-4-16

SIGNAL FILTERING AT JUMP MULTIPLICATIVE INTERFERENCE

У.А. Bukhalev1 vadim.bukhalev@yandex.ru

А.А. Skrynnikov2' 3 a1260@mail.ru

У.А. Boldinov3 viktorboldinov@mail.ru

1 Moscow Scientific-Research Television Institute, Moscow, Russian Federation

2 State Research Institute of Aviation Systems, Moscow, Russian Federation

3 Moscow Aviation Institute, Moscow, Russian Federation

Abstract

The problem of constructing the approximately optimal Bayesian algorithm of filtering the output signal for a linear stochastic dynamical system is solved. Non-linear mixture of the output signal and of the multiplicative interference was measured. This interference is a continuous-valued random process with the unknown distribution law within the given limits of [c(1), c(2)], c(1) > 0, c(2) > 0, and in the frequency range of [0, Ara]. The filtering algorithm is synthesized by the approximately optimal signal estimation method based on the theory of systems with random jump structure and the method of two-time parametric approximation of the phase coordinates probability distribution. The method consists in approximate replacement of the unknown probability densities of phase coordinates by the known distribution laws with the unknown mathematical expectations and covariances determined as a result of solving the problem. The proposed approximate-optimal filtering algorithm is based on replacement of the continuous-valued multiplicative interference by a random jump process, i.e., Markov chain with two states c(1), c(2) and equal intensities of transitions from one state to another q' = 2Ara. Conditional probability densities of the output signal for fixed states of the Markov chain c(1), c(2) are approximated by gamma distributions that depend on mathematical expectations and variances of the x(t) signal at the fixed c(1), c(2) and measurements x(t) in a mixture with the Markov jump interference. An example of constructing an algorithm for evaluating the distance from an aircraft to the object based on object irradiance measurement by the infrared direction finder was considered. Irradiance was equal to the radiation strength ratio to the square of the distance. Radiation strength, i.e., multiplicative interference, was the random process with unknown continuous distribution in a limited range. The approximate-optimal filter was constructed for distance evaluation by replacing this process with a two-state Markov chain and approximating the distance probability density by the Rayleigh distribution

Keywords

Information system, random jump structure, multiplicative jump interference, Bayesian filtering, two-time parametric approximation, object infrared direction finder

Received 04.06.2022 Accepted 15.02.2023 © Author(s), 2023

The work was financially supported by the Russian Science Foundation (grant no. 22-29-00708)

REFERENCES

[1] Pierce B., Sworder D. Bayes and minimax controllers for a linear system with stochastic jump parameters. IEEE Trans. Automat. Contr., 1971, vol. 16, iss. 4, pp. 300-307. DOI: https://doi.org/10.1109/TAC. 1971.1099732

[2] Sworder D.D. Bayes controllers with memory for a linear system with jump parameters. IEEE Trans. Automat. Contr., 1972, vol. 17, iss. 1, pp. 119-121.

DOI: https://doi.org/10.1109/TAC. 1972.1099877

[3] Mariton M. Jump linear systems in automatic control. Taylor & Francis, 1990.

[4] Kazakov I.E., Artemyev V.M., Bukhalev V.A. Analiz sistem sluchaynoy struktury [Analysis of systems of random structure]. Moscow, FIZMATLIT Publ., 1993.

[5] Pakshin P.V. Diskretnye sistemy so sluchaynymi parametrami i strukturoy [Discrete systems with random parameters and structure]. Moscow, FIZMATLIT Publ., 1994.

[6] Bukhalev V.A. Raspoznavanie, otsenivanie i upravlenie v sistemakh so sluchaynoy skachkoobraznoy strukturoy [Recognition, evaluation and control in systems with a random jump-shaped structure]. Moscow, Nauka Publ., 1996.

[7] Bukhalev V.A. Osnovy avtomatiki i teorii upravleniya [Fundamentals of automation and control theory]. Moscow, VVIA im. N.E. Zhukovskogo Publ., 2006.

[8] Bukhalev V.A. Optimalnoe sglazhivanie v sistemakh so sluchaynoy skachkoobraz-noy strukturoy [Optimal smoothing in systems with a random discontinuous structure]. Moscow, FIZMATLIT Publ., 2013.

[9] Zhang C., Zhu H., Zhou H., et al. Deterministic and stochastic differential games. In: Non-Cooperative Stochastic Differential Game Theory of Generalized Markov Jump Linear Systems. Cham, Springer, 2017, pp. 17-29.

DOI: https://doi.org/10.1007/978-3-319-40587-2_2

[10] Bukhalev V.A., Skrynnikov A. A., Boldinov V.A. Algoritmicheskaya pomekhoza-shchita bespilotnykh letatelnykh apparatov [Algorithmic noise protection of unmanned aerial vehicles]. Moscow, FIZMATLIT Publ., 2018.

[11] Bukhalev V.A., Skrynnikov A.A., Boldinov V.A. Igrovoe upravlenie sistemami so sluchaynoy skachkoobraznoy strukturoy [Game control of systems with a random jump structure]. Moscow, FIZMATLIT Publ., 2021.

[12] Sebryakov G.G., Krasilshchikov M.N., eds. Upravlenie i navedenie bespilotnykh manevrennykh letatelnykh apparatov na osnove sovremennykh informatsionnykh tekhnologiy [Control and guidance of unmanned maneuverable aircraft based on modern information technologies]. Moscow, FIZMATLIT Publ., 2005.

[13] Bukhalev V.A., Boldinov V.A., Pryadkin S.P., et al. Two-moment parametrical approximation of distributions in filtering. Vestnik kompyuternykh i informatsionnykh tekhnologiy [Herald of Computer and Information Technologies], 2016, no. 8, pp. 8-15 (in Russ.). DOI: https://doi.org/10.14489/vkit.2016.08.pp.008-015

[14] Johnson N.L. Continuous univariate distributions. Wiley, 1994.

[15] Korolyuk V.S., ed. Spravochnik po teorii veroyatnostey i matematicheskoy statis-tike [Handbook of probability theory and mathematical statistics]. Moscow, Nauka Publ., 1985.

[16] Lazarev L.P. Infrakrasnye i svetovye pribory samonavedeniya i navedeniya le-tatelnykh apparatov [Infrared and light devices for homing and guidance of aircraft]. Moscow, Mashinostroenie Publ., 1976.

Bukhalev V.A. — Dr. Sc. (Eng.), Chief Researcher, Moscow Scientific-Research Television Institute (Golyanovskaya ul. 7a, str. 1, Moscow, 105094 Russian Federation).

Skrynnikov A.A. — Cand. Sc. (Eng.), Head of the Sector, State Research Institute of Aviation Systems (ul. Viktorenko 7, Moscow, 125319 Russian Federation); Assist. Professor, Department Theory of Variables and Computer Modeling, Moscow Aviation Institute (Volokolamskoe shosse 4, Moscow, 125993 Russian Federation).

Boldinov V.A. — Cand. Sc. (Eng.), Assist. Professor, Department of Aviation Robotic Systems, Moscow Aviation Institute (Volokolamskoe shosse 4, Moscow, 125993 Russian Federation).

Please cite this article in English as:

Bukhalev V.A., Skrynnikov A.A., Boldinov V.A. Signal filtering at jump multiplicative interference. Herald of the Bauman Moscow State Technical University, Series Natural Sciences, 2023, no. 2 (107), pp. 4-16 (in Russ.). DOI: https://doi.org/10.18698/1812-3368-2023-2-4-16

i Надоели баннеры? Вы всегда можете отключить рекламу.