Раздел III. Системы управления
УДК 519.283
Е.Я. Рубинович
РОБАСТНАЯ ЛИНЕЙНАЯ ФИЛЬТРАЦИЯ СКАЧКООБРАЗНЫХ
ПРОЦЕССОВ*
Для процессов, описываемых линейными стохастическими дофференциалъными уравнениями Ито с мерой, рассматривается задача оптимальной (в среднеквадратическом смысле) линейной фильтрации в случае, когда частично наблюдаемый векторный процесс имеет ненаблюдаемую компоненту возмущенную скачкообразным марковским процессом с конечным числом состояний, а процесс наблюдений зашумлен семимартингалом с негауссовской мартингалъной частью. Иллюстрируются робастные свойства оптимальных линейных оценок в сравнении с оптимальными нелинейными оценками в случае неадекватности математической модели реальному процессу.
Марковские процессы; фильтрация; робастность.
E.Ya. Rubinovich ROBUST LINEAR FILTRATION FOR JUMP-WISE PROCESSES
The paper deals with an optimal (in the mean square sense) linear estimation for the processes describing by Ito's linearly differential equations with measure when partially observable vector-process has unobservable component corrupted by Markov jump-wise process with finite number of states and observation process is excited by semimartingale with non-Gaussian martingale part. It is shown robust property of linear estimations.
Markov processes; filtering; robustness.
Введение. Во многих технических задачах возникает необходимость оценки координат скачкообразных процессов, наблюдаемых на фоне шума. Такая задача обычно решается с применением фильтров Калмана, что связано с построением определенной вероятностной модели описания эволюции оцениваемых координат во времени, которая учитывала бы скачкообразный характер процесса. В качестве такой модели в работе рассматривается скачкообразный марковский процесс с конечным числом состояний. Отметим, что для построения оценок состояний такого процесса в настоящее время существует хорошо разработанный аппарат теории нелинейной фильтрации [2]. Однако, с практической точки зрения, оптимальный (в среднеквадратическом смысле) линейный фильтр оказывается предпочтительнее. Дело в том, что кроме своей простоты, с вычислительной точки зрения, линейный фильтр «грубее» оптимального нелинейного фильтра в том смысле, что оценки даваемые линейным фильтром в случае неадекватности математической модели реальному процессу (например, при аппроксимации реального процесса марковским процессом с малым числом состояний, неточным знанием интенсивностей переходов между состояниями и т.п.) оказываются более точными.
*
Работа выполнена при поддержке Программы Президиума РАН «Динамические системы и теория управления».
Цель настоящей работы - построить оптимальный в среднеквадратическом смысле линейный фильтр для оценки скачкообразного марковского процесса с дискретным множеством состояний наблюдаемого на фоне шума и на базе этого фильтра синтезировать алгоритм отслеживания траекторий реальных скачкообразных процессов, в частности, фазовых координат маневрирующих динамических объектов.
Естественный подход к решению данной задачи состоит в следующем. По среднему значению и корреляционной функции исходного марковского процесса строится его гауссовский марковский аналог, т.е. гауссовский процесс с теми же двумя моментами, что и у исходного процесса. По этому гауссовскому аналогу строится калмановский фильтр, являющийся оптимальным линейным для исходного процесса [1].
Трудность построения такого аналога состоит в сложности вычисления корреляционной функции марковского процесса с большим числом состояний (в случае двух состояний - это простая задача). Мы преодолеваем эту трудность, прибегая к так называемому, семимартингальному представлению скачкообразного марковского процесса [3] и строим его гауссовский марковский аналог без вычисления среднего значения и корреляционной функции. Такой подход (в дискретном и непрерывном времени) был успешно применен в простейшем случае оценивания одномерного скачкообразного марковского процесса с конечным числом состояний, наблюдения за которым были искажены аддитивным белым шумом [1]. В настоящей работе, кроме обобщения на многомерный случай, рассматривается более широкий класс случайных процессов, являющихся семимартингалами, эволюция состояния которых описывается стохастическими дифференциальными уравнениями с мерой. Такое представление позволяет единообразно описывать дискретные, непрерывные и дискретно-непрерывные процессы.
Постановка задачи. На стохастическом базисе (О. Р. Е = (^У)о<г<т, Р) с обычными условиями задан ненаблюдаемый векторный процесс X = (Х(*))(>о е й"1, удовлетворяющий стохастическому дифференциальному уравнению Ито:
Наблюдению доступен векторный процесс V = (У(())^0 € 14.*:
Здесь X(/—) = ПтX(в), (; (} = (СЛ/))(>0 — векторный марковский про-
цесс с конечным числом состояний; 2 = (2(£))(>о — квадратично интегрируемый мартингал (относительно Г), независимый от С}, с 2'(0) = 0 и квадратической характеристикой {^)£; <!«(■) и (1.4(-) — детерминированные матричные меры по t. Все функции от t предполагаются непрерывными справа и имеющими пределы слева. Предполагается также существование матричной функции удовлетворяющей при I € [0,Т] следующим условиям согласованности:
где верхний индекс c - соответствует непрерывным компонентам мер; - знак
матричного следа; * - символ транспонирования; Е - знак математического ожидания. Задача состоит в построении оптимального (в среднеквадратическом смысле) линейного фильтра для частично наблюдаемого процесса (А", У').
Основные результаты. Известно, что оптимальный линейный фильтр для
негауссовского процесса (А”, У) является оптимальным фильтром для его гауссовского аналога (.с, у), т.е. гауссовского процесса с теми же самыми математическим ожиданием и матрицей ковариаций [3]. Для получения представления процесса (.и. у) будем использовать семимартингальное представление компонент процесса <7, где мартингальные части заменим на винеровские мартингалы. Таким же образом мы заменим мартингал Ъ на гауссовский мартингал г — (г(0) с т°й же квадратической характеристикой.
Пусть (,>„(*), п = 1 ,то (п-я компонента процесса <7) имеет вектор состояния / I .......... ' и матрицу интенсивностей переходов между состояниями
с элементами А?-(4) - известными непрерывными функциями времени. Как ив [1], определим индикаторный процесс /п(/) — 1 (^))* с компонентами
/"(О = 1{<У„(0 = а-"}, г = 1,кп - 1, где 1{ } - индикатор множества {*}. Тогда
Индикаторный процесс 1п = (/„(*))(>о допускает семимартингальное представление [1]: ^ ^
где
с ^
и Мп и) — (М" (?),..., _!(*))* - квадратично интегрируемый, чисто разрыв-
ный мартингал по отношению к семейству Е7'* = (■Р/’1 )г>о 'Т-алгебр =<т{/п(5), 5</}.
Как отмечалось выше, гауссовский аналог х процесса X удовлетворяет уравнению ^
где д(г) — (§1 ((),..., Ят{1))* ~ гауссовский аналог процесса (^(/ ). Вследствие (3) где ” '
Ш = Ш+[ (\'п(з)1п(з)+\п(з))<1з+ [ дУ2(з)а*вд,
•'[о,*) ./[о,<]
с винеровским процессом И„ (я) = (VI "(5),.... №^-1(5)) и матричной функцией
сЦЕ(М),]
Я„(<) =
СІ?
Здесь (Л/) ( - квадратическая характеристика мартингала М (і). Матрица // . ■ имеет элементы [1]:
Я£(0 = Еіа„ |А£(<)1р?(<),
Я&М = -(А5(і)р?(*) + Аг Ф где р" (/) = Р{(}п{і) = а'?} могут быть найдены из уравнения Колмогорова
pj(t) = E W(<)-
i<kr}
Далее введем следующие блочно-диагональные матрицы:
а также блочные векторы
/(£) = соЫтп^ :/т(*)],
А(<) = со1ишп[Ах(<):... :Ат(<)], и обозначим л = ,...,а^)*. Тогда для
■' ■ 11 г. I. /' / ■ имеет место представление
блочного вектора
[ R(s)dW(s),
J[0.t)
(4)
где О - матрица с нулевыми элементами, Е - единичная матрица и IV = (11’(Т)),>о - векторный винеровский процесс, размера (1ии1Г ~ к^ + ■ ■ ■ + кт — т.
Теперь заменим мартингал 2 на гауссовский мартингал г = (.'(^)),>|;1 с квадратической характеристикой (г}< = Е{^)(. После такой замены гауссовский аналог у процесса У принимает вид
(х^-) \
' (5) /(*) /
Для гауссовского процесса (и, у) справедливы уравнения обобщенного фильтра Калмана [4], которые описывают эволюцию оптимальной (в среднеквадратическом смысле) оценки 7г((, у) ненаблюдаемого процесса г/ по наблюдениям за процессом у, и эта оценка оказывается линейной. Для аналогичной оптимальной линейной оценки У) процесса Л* по наблюдениям У имеет место следующая теорема.
Теорема 1. Оптимальная линейная оценка ,Х(#,У) процесса А" по наблюдениям У равна первым т компонентам вектора тг(4, У):
где (к\ +-----кт — т) - размерный векторный процесс тг(£, У) удовлетворяет
уравнениям обобщенного фильтра Калмана с заменой процесса у на реальный процесс наблюдений У.
Пример. Приведем пример, который демонстрирует анонсированные в заголовке робастные свойства оптимальных линейных оценок. Как отмечалось, эти свойства особенно проявляются тогда, когда реальный ненаблюдаемый процесс значительно отличается от рассматриваемой математической модели. С такой ситуацией приходится, в частности, сталкиваться в задачах оценки элементов движения маневрирующих динамических объектов.
На рис. 1 и 2 представлены оценки, даваемые оптимальным нелинейным (рис. 1) и линейным (рис. 2) фильтрами, априорно рассчитанными для оценки состояний скачкообразного марковского процесса с четырьмя наблюдаемыми на фоне шума состояниями (гц, П2,«з, сц), значения которых приведены на рисунках. Однако реальные наблюдения осуществляются за процессом, отмеченным на рисунках сплошной плавной кривой. Как видно из рисунков, там, где реальный процесс мало отличается от предписанных состояний (а,і и г/з на рисунках), линейный и нелинейный фильтры дают удовлетворительные оценки. Там же, где траектория реального процесса проходит между предписанными состояниями, оценки линейного фильтра оказываются намного точнее. «Обманутый» оптимальный нелинейный фильтр «не понимает» промежуточного состояния и начинает «метаться» между предписанными состояниями (02 и о і на рисунках).
Рис. 1. Оцениваемая траектория и оценка оптимального нелинейного
Рис. 2. Оцениваемая траектория и оценка оптимального линейного
фильтра, рассчитанного на состояния фильтра, рассчитанного на состояния
(«1, 0-2* йз, «4) скачкообразного процесса («1, а 2. аз, «4) скачкообразного процесса
Для сравнения на рис. 3 представлена оценка, даваемая оптимальным линейным фильтром, рассчитанным на два состояния (Ь], 62), где
Рис. 3. Оцениваемая траектория и оценка оптимального линейного фильтра, рассчитанного на два состояния (Ь-[, Ь%) скачкообразного процесса
Приложение. В этом разделе приводятся уравнения обобщенного фильтра Калмана для рассматриваемых в работе процессов. Для простоты мы преобразуем уравнения системы (1), (2) к системе уравнений с общей скалярной мерой [4].
Обозначим тага(() = . таг (У) сумму вариаций элементов матрицы «(л)
для е [0,4 Аналогичные обозначения будем использовать для матриц А(4) и Г. / . Положим
К*) = * + уаг а(£) + уаг т + уаг Е(.£)( и пусть д{Ь) - любой из элементов матриц а((), А(1). Е(Х)%. Тогда имеет место абсолютная непрерывность мер (!<?(-) -С Далее, обозначим через Ьу(^) про-
изводную Радона-Никодима меры <1я^(-) по мере (1/4 У а через Ь - матрицу с элементами Ьу{£), *, 3 — 1,т. Тогда <1а(1) = Ь(1)Л^(1). Аналогично представляем
сМ(^) = В(£)с1^(£),
(ВД, = №М1),
с!< = 0(<)с1/хс(4).
В новых обозначениях система (1), (2) переписывается в виде
Уравнения (4), (5) принимают соответственно вид
(6)
(7)
где w = (w(f))t>o - винеровский мартингал с квадратической характеристикой «1 if , к| / Id/.M i ) и
«№ = ( — ) R№'(t) ( j Ht);
Wn =
/ b(t) I О h(t) = --------I-----------------
V О І Л (t)0(t)
H(t) = [B(0io].
Для системы (6), (7) справедливы следующие уравнения обобщенной калма-новской фильтрации [4]:
,У)Ац(8)+
і
[
о ,t]
T(s)[dY(s) - Н(з)п(з-,ГУЩз)],
где ttq — Е«о и
г(0 = [я + л(*)М*)Ь(*-)я’(*)[0(<)+
Здесь «+» - символ псевдообращения, Afi(t) = n(t) - fi(t—), а матричнозначная функция 7(f) удовлетворяет обобщенному уравнению Рикатти с мерой
7(f) =1о+[ [*(«) + h(s)~t(s) -I- 7(s)/i*(s)]d/i(s)+
J]o,t]
+ [ r(s)[p(s) + H(s)~,(s-)H*(s)A,,(s)}r*(s)dLi(s)+
+ ^2 л(т*)7(п-)Л*(г4)(Дд(п))2,
Ti<t
где 7o = cov (7T0,7Tq), a r,-. t = 1.2.,.. - детерминированные моменты скачков меры <!/*(■)•
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Krichagina N.V., Liptser R.Sh. andRubinovich E.Ya. Kalman filter for Markov processes // In: Statistics and control of stochastic processes. - New York: Publ. Div., 1985. - P. 197-213.
2. Капылов А.К. О нелинейной фильтрации случайных процессов при дискретно поступающих данных // Problems of Control and Information Theory. - 1979. - Vol. 8 (1). - P. 39-54.
3. Липцер Р.Ш., Ширяев А.Н. Статистика случайных процессов. - М.: Наука, 1974. - 696 с.
4. Miller B.M., Rubinovich E.Ya. Regularization of a generalized Kalman filter // Mathematics and Computers in Simulation. - 1995. - Vol. 39. - P. 87-108.
Статью рекомендовал к опубликованию д.т.н. Е.П. Маслов.
Рубинович Евгений Яковлевич - Федеральное государственное бюджетное учреждение науки Институт проблем управления им. В.А. Трапезникова Российской академии наук; e-mail: rubinvch@ipu.rssi.ru; 117997, г. Москва, ул. Профсоюзная, 65; тел.: 84953349111; д.т.н.; профессор; зам. директора по научной работе.
Rubinovich Evgeny Yakovlevich - V.A. Trapeznikov Institute of Control Sciences of Russian Academy of Sciences; e-mail: rubinvch@ipu.rssi.ru; 65, Profsoyuznaya street, Moscow, 117997, Russia; phone: +74953349111; dr. of eng. sc.; professor; deputy director on R&D.
УДК 681.5-192
А.Б. Филимонов, Н.Б. Филимонов
АВТОМАТ ОГРАНИЧЕНИЙ УПРАВЛЯЕМЫХ ДИНАМИЧЕСКИХ
ПРОЦЕССОВ*
Обсуждается концепция гомеостатического управления, основанная на распространении феномена гомеостаза живых организмов на технические системы. Технический смысл гомеостаза заключается в ограничении контролируемых динамических параметров системы. Данные функции возлагаются на автомат ограничений. Предлагаются механизмы функционирования автоматов ограничений на основе концепции динамических барьеров в фазовом пространстве системы и стратегий упреждающего гомеостатического управления. Развивается полиэдральная методология синтеза автомата ограничений.
Гомеостаз в технических системах; автомат ограничений, динамические барьеры; упреждающее гомеостатическое управление; полиэдральный подход.
*
Работа выполнена при поддержке РФФИ, грант № 13-08-00161.