Файл: Липцер Р.Ш. Статистика случайных процессов. Нелинейная фильтрация и смежные вопросы.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 10.04.2024

Просмотров: 221

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

§ 21

а с и м п т о т и ч е с к и е с в о й с т в а ф и л ь т р а

617

ненаблюдаемых состояний Ѳ,. При этом, как и в аналогичной задаче для случая дискретного времени, возникает важный для приложений вопрос об условиях, когда матрицы y t стаби­ лизируются при t I оо. Исследованию вопроса о существовании предела lim yt и способах его вычисления и посвящен настоя-

t - > оо

щий параграф.

2, Прежде чем переходить к точным формулировкам, за­ метим, что, полагая

а = а ] — ( Ь о В ) ( В о В ) - 1 Л „

Ь =

[(ЬоЬ) — (Ь° В) (В о ß)“1(b о ß)*]1/2,

(16.43)

ß =

[ßoß]1/2,

Л = Л„

 

уравнение (16.42)

можно переписать в несколько более удобном

виде:

аУі + Y<о* +

 

(16.44)

Уі =

ЬЬ* — y tA*{BB*)~l A y t.

Это уравнение совпадает с уравнением для ковариации при рассмотрении гауссовской пары процессов (Ѳ, |), удовлетворяю­ щих системе

dQt = аѲ, dt + b d W x (t),

(16.45)

d \ t — AQt dt Ar В dW 2 (t).

Так что с точки зрения исследования поведения матриц yt при

оо достаточно вместо системы (16.39) рассматривать более простую систему (16.45).

Т е о р е м а

 

16.2.

Пусть

рассматривается

система (16.45),

для которой выполнены следующие условия:

 

(I) ранг блочной матрицы

 

 

 

 

 

 

 

 

/ А

\

 

 

 

 

 

G\

Аа

 

(16.46)

 

 

 

 

 

 

 

 

 

 

 

Aak~l

 

размерности

{kl X k)

равен k\

 

 

(II) ранг блочной

матрицы

 

 

 

 

 

 

G2 = (b

ab . ,.

ak~lb)

(16.47)

размерности {k X Ik) равен k\

 

 

(III) матрица B B *

не вырождена.

 

Тогда для

у. =

М (Ѳ, — m t) ф, m tY существует lim у ( = у.

 

у

 

 

 

 

 

<-»оо

Этот предел

не

зависит от начального значения у0 и является


618

 

ПРИМЕНЕНИЕ К НЕКОТОРЫМ ЗАДАЧАМ УПРАВЛЕНИЯ

 

[ГЛ. 16

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

единственным (в классе положительно определенных матриц)

решением

уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ау + Уа*+

bb* у А* (В 'В )~ 1Ау =

0.

 

(16.48)

 

Доказательству этой теоремы предпошлем ряд вспомога­

тельных утверждений.

 

 

 

 

 

 

А — матрицы

 

 

 

3.

(k

Л е м м а

16.2.

Пусть

D

 

и

размерностей

(I X k),

X k).

Образуем

блочную

 

матрицу

(порядка (nl X k))

 

 

 

 

 

 

 

 

 

 

( D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

VDA'1"'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т

 

 

 

 

 

 

 

 

 

Тогда

матрицы

D*nDn

и

[ e~AHD*De~At dt,

0 < Г < о о ,

одно-

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

временно либо вырождены, либо невырождены.

 

 

 

 

Д о к а з а т е л ь с т в о .

Согласно

лемме 14.4 матрицы

D*nD n

и D \D k ,

n ^ k ,

 

одновременно

или

вырождены

или невыро­

ждены.

Если матрица D \D k

вырождена, то по этой же лемме

найдется ненулевой вектор x =

( x t ,

. . . ,

х п) такой,

что D A ’ x

— Q,

} — 0, 1, .. .,

k,

k

+

1, .. .

Но тогда

 

 

 

 

 

 

 

 

 

 

 

 

D e - ^ x = \ ( . f (Е>АД) = 0,

 

 

 

 

и,

следовательно,

 

т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

J e - ^ D ' D e - v

d t x = 0,

 

 

 

(16.49)

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

что и

доказывает вырожденность

 

матрицы

|

e~A*tD’,De~At dt.

 

Наоборот,

 

пусть

выполнено

(16.49).

0

 

очевидно,

 

 

Тогда,

x ' e - ^ D ' D e - M х =

0,

0 < / < 7 \

Поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

De~Atx = 0

 

 

 

 

(16.50)

и (после дифференцирования по t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D Ае~АІ X =

0,

 

 

 

 

(16.51)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

DAk~^e~st =

0.

 

 

 

 

 

 

Из

(16.50)

и

(16.51)

при

f = 0

 

вытекает,

что

DAjx = 0,

/ =

0, ... ,

k — l,

что эквивалентно

равенству

x'DIDkX — 0.

 

Лемма доказана.

 

 

 

 

 

 

 

 

 

 

 

 

 


§2J

 

 

Ас и м п т о т и ч е с к и е с в о й с т в а ф и л ь т р а

619

Сле д с т в и е .

Пусть

D k =

(D

/S.D ...

Аk~ l D ) — блочная

ма­

трица

порядка (k X kl),

где

D

и

А — матрицы размерностей

 

 

 

 

 

 

 

 

 

 

т

 

(k X I)

и

( k y , k ) .

Тогда

матрицы

ЪкЪ\

и

\ e~htDDte - h'’t dt

 

 

 

 

 

 

 

 

 

 

ö

 

одновременно либо вырождены, либо не вырождены.

 

Л е м м а

16.3.

Если матрица

G2

имеет ранг k, то при t > О

матрицы

y t,

определяемые

из

уравнения

(16.44), являются

положительно определенными.

 

 

y t является матрицей

ко­

Д о к а з а т е л ь с т в о .

Матрица

вариаций

условно-гауссовского

распределения

 

Если это распределение имеет (Р-п. н.) плотность, то тогда, очевидно, матрица y t будет положительно определенной. Рас­

сматривая систему уравнений (16.45) и принимая во внимание следствие 1 теоремы 7.23 (п. 5 § 9 гл. 7), получаем, что рас­ пределение Р (Ѳ, а |^"|), t > 0, имеет плотность (Р-п. н.), если

плотностью обладает распределение Р(Ѳ ,^а), что эквивалентно условию положительной определенности матрицы Г/ = соѵ(0/, Ѳ,).

Согласно теореме 15.1 матрицы Г, являются решениями дифференциального уравнения

Г, = аТ,

+ Г,а* + ЪЪ\

(16.52)

Отсюда находим

 

 

Г; = eatTüeaH + eat

e~asbb*e~a*s ds

0a*t

Но в силу следствия леммы

16.2 матрицы Г„

t > 0, положи­

тельно определенные, поскольку таковой же является и ма­

трица G2G2

(rang G2 = k).

Лемма доказана.

 

Л е м м а

16.4. Если

ранг матрицы G{ равен k, то элементы

всех матриц у ( равномерно ограничены.

Д о к а з а т е л ь с т в о .

Рассмотрим вспомогательную задачу

управления детерминированным процессом xf = {x{ (t), . .. , xk (t)),

Оудовлетворяющим уравнению

 

dx

A*ut, х0 = х,

(16.53)

 

— a*xt +

с функционалом

 

т

 

 

 

 

V (и;

Т) х"Гу0хт+

С[x]bb*xt + и)ВВ*и^ dt.

 

 

 

о

 

Управления ut,

выбираются из класса допустимых

(см предыдущий параграф).

 

 


620 ПРИМЕНЕНИЕ К НЕКОТОРЫМ ЗАДАЧАМ УПРАВЛЕНИЯ [ГЛ. 16

Согласно (16.37) оптимальное управление ü t существует и

задается формулой

 

ü t = (ВВ*)~' Аут-tXf,

(16.54)

где

xt— решение уравнения (16.53) с ut = ü t , 0

Г. При

этом

V (Гг, Т) = х*утх. Поскольку элементы матриц y t

являются

непрерывными функциями, для доказательства леммы доста­

точно показать, что равномерно ограничены все

элементы

матриц ут при Т >

 

1.

 

ß,

матрица

GJGi не вырождена, и по

Поскольку rangG1=

 

лемме 16.2

не вырождена матрица

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

0

 

e - aHA*Ae-atdt.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Возьмем теперь

специальное управление

 

 

 

 

J

 

 

 

 

 

 

 

 

 

 

Ае~аЧ f е - ^ Л ’А е - 03 | х,

0 < Н < 1 ,

 

üt — '

 

 

 

0,

 

/

 

t >

1,

 

 

и пусть £( — решение уравнения

(16.53) с

ut =

üt. Решая это

уравнение,

находим,

что

£t = 0,

Ѵ ^ \ .

Но тогда

в

силу

опти­

мальности управления üt,

O ^ t ^ T ,

Т >

1,

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

х*утх ^ . I \%tbb*%t +

utBB*üt]dt

< оо,

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

что и доказывает лемму.

 

 

 

 

 

 

 

с у°==

Л е м м а

16.5. Пусть

 

y°t решение уравнений (16.44)

= у0 = 0 и

rangG, =

/j.

 

Тогда существует

lim y°t

и

у0 —

limy?

 

 

 

 

 

 

 

 

 

 

/- » о о

 

 

/- » о о

является неотрицательно определенной симметрической матри­ цей, удовлетворяющей уравнению

ауо + у0а*+

ьь* _ ѵоЛ* (ßß*)'1Луо =

0.

(16.55)

Если к тому же rang G2 — k,

то у0 является

положительно

определенной матрицей.

В силу

предположения

rang G{ — k

Д о к а з а т е л ь с т в о .

из предыдущей леммы следует, что элементы всех

матриц у\,

t ^ O , равномерно ограничены.

 

 

 


§ 2]

а с и м п т о т и ч е с к и е с в о й с т в а ФИЛЬТРА

621

Покажем, что при у0 — 0 функция х*у\х монотонно не убы­

вает (по Т). Пусть 7\,>7Ѵ Тогда, обозначая мДД) и x t (T{)

оптимальные управления и отвечающие им процессы во вспо­ могательных задачах управления, і — 1, 2, находим, что

г,

х у тх = J [(*, (Т2)У bb*xt (T2) + (Т2)У ВВ*щ (Т2)} dt >

о

> JГ,[(*, (Т2)У bb*x, (Т2) + (Щ (Т2)У ВВ 'щ (Гг)] dt >

о

ТI

> J [(*, (Г,))* bb*xt (Т,) + (Щ (Т})У ВВ' ЩІТУ) dt = х у Т[х.

о

Из ограниченности и монотонности функций х*у°тх вытекает существование матрицы у0 = lim у°т с указанными свойствами.

Г - > ОО

Если же дополнительно rang G2 — k, то по лемме 16.3 не вырождены матрицы y°t, а следовательно, таковой же является

и

матрица у° = lim у?.

 

 

 

 

 

 

4.

/

оо

 

те оре мы 16.2.

Обозначим у° =

=

Д о к а з а т е л ь с т в о

lim у? при уо =

0 и

положим

 

 

 

 

 

 

t->oo

 

 

 

 

 

 

 

 

 

 

 

 

Щ = - ( В В Т '

Ay°xt,

 

 

(16.56)

где xt — решение уравнения (16.53)

с ut =

üt, х0 =

х.

Покажем,

что лу—>0, t —>оо.

Для этого достаточно,

например,

показать,

что

 

 

lim #у°л:г == 0,

 

 

(16.57)

 

 

 

 

 

 

 

 

 

 

t со

 

 

 

 

 

поскольку матрица у0

является симметрической и положительно

определенной.

 

 

 

 

 

 

 

 

В силу (16.53), (16.55) и (16.56)

 

 

 

 

— (x*ty°xt) = х\уо [а* -

А" ( В В Т '

Лу°1 x t +

 

 

 

 

+ xt [а — у° А * ( В В Т ' Иу°] X, — #уМ ’( В В У ' А у % =

 

=

- xtbb'xt -

xty°Л* ( B B T '

{ВВ*) ( В В Т 1Ау°х( =

 

 

 

 

 

 

 

=

[x*fbb*xt + ГцВВ*щ\ч

/