Файл: Липцер Р.Ш. Статистика случайных процессов. Нелинейная фильтрация и смежные вопросы.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 10.04.2024
Просмотров: 246
Скачиваний: 0
§ 4] |
|
АСИМПТОТИЧЕСКИЕ СВОЙСТВА |
ФИЛЬТРА |
|
|
559 |
||||
где |
т (t, t + |
1) = |
М [Ѳ( | |
Следовательно, |
|
|
|
|||
|
+ О — mt+\ = |
\ai — (bo В) (BoB)+ |
— |
t + |
1)) = |
|
||||
и |
|
|
|
|
|
|
= a(Qt — m(t ,t + |
1)), |
||
|
|
|
|
|
|
|
|
|
|
|
M[(m6((^+ |
1, t) — mt+l) (m6i(t+ |
1, t) — m, + 1)*| = ay (t,t + 1)a , |
||||||||
где |
y (^ t + |
1) = |
М[(Ѳ<— m(t, t + |
1))(0/ — m(t, t + |
1))*]. |
Но |
со |
|||
гласно (13.110) |
|
|
|
|
|
|
|
|
||
|
Y (*, t + |
1) = yt — ytA\ [ßoß-f- H[Y^I]+ ^iYr |
|
|
||||||
|
Значит |
для |
yt, |
t > 0, |
также справедливо |
рекуррентное |
||||
уравнение |
|
|
|
|
|
|
|
|
|
Y/ + 1 = k - ( b o В ) ( В о В ) + Л ,] yt [а, ~ ( Ь о В ) ( В о В ) + Л, f +
+[Ь° b — (Ь о В) (В о В)+ (Ь оЯ)*] —
—[а, — (ft о Я) (Я ОВ)+ (Ь оß)*] Ѵ(л; X
X [Я о В + Л . у ( Л | ] + Л , у , [ а , - (Ь о Я) (Я о fi)+ (b о В)*]\
Поэтому в дальнейшем будет рассматриваться лишь система (14.100) и соответственно изучаться асимптотическое поведение матриц yt, удовлетворяющих рекуррентному уравнению
Yt+i = ayta* + bb* - aytA* [BB* + AytA*)+Ayta. (14.103)
Т е о р е м а 14.3. Пусть выполнены следующие условия'.
(I) Ранг блочной матрицы
..-(V )
|
|
|
|
|
|
\ |
Aak- 4 |
|
|
размерности |
(kl X k) |
равен |
k\ |
G2 — (b ab ...ak~xb) |
размерности |
||||
(II) |
ранг |
блочной |
матрицы |
||||||
(k X Ik) |
равен k; |
не вырождена. |
|
|
|||||
(III) |
матрица B B * |
Y0- |
При этом |
||||||
Тогда |
существует и |
не зависит от уо Ншу< = |
|||||||
Sp у0 < |
00 |
и |
матрица |
у0 |
|
t -> ОО |
|
решением |
|
является единственным |
(в классе симметрических положительно определенных матриц) матричного уравнения
у = ауа*+ bb* - ауА* (ВВ* + А у А Т 1 Луа*. |
(14.104) |
2. Доказательству этой теоремы предпошлем ряд вспомо гательных утверждений.
|
|
АСИМПТОТИЧЕСКИЕ СВОЙСТВА |
ФИЛЬТРА |
561 |
||||
Пусть |
і — max [/<;&: ct Ф 0]. |
Тогда |
из (14.109) |
получаем |
||||
|
|
|
г—1 |
|
|
£/ |
|
|
|
|
|
Уі= /=0 |
|
= |
|
|
|
|
|
|
|
Ci’ |
|
|||
и, следовательно, |
і-1 |
|
і-\ |
' |
|
|||
|
|
|
|
|
||||
|
|
Dk — dk |
^ i — ^,i°jdk |
1у j = |
|
iHk- і+і. |
|
|
|
|
|
і—о |
|
/=о |
|
|
|
Поэтому в силу (14.108) |
|
|
|
|
|
|||
|
|
Dd*x = D y * J j i c ' Dyk_ |
|
= 0 . |
|
|||
|
|
|
|
і=О |
|
|
|
|
Отсюда по индукции устанавливаем, что |
Dd‘x — 0, j ^ k , что |
|||||||
вместе |
с |
формулой |
(14.105) |
доказывает |
утверждение |
леммы. |
||
С л е д с т в и е . Пусть D = |
D(kxi), d — dikXk)—некоторые матри |
|||||||
цы и |
|
|
|
|
|
|
|
|
|
|
D n — ( D d D , . . . d * - ' D ) |
|
|
||||
—блочная |
матрица |
порядка (k X til), |
n ^ k . Тогда |
матрицы |
DnDn и DkDk одновременно либо вырождены, либо не вырождены.
Л е м м а |
14.5. |
Пусть |
Ѳ= [Ѳ, (/),..., |
ОД/)], |
/==0,1, .... |
||||||||
—гауссовская последовательность, удовлетворяющая |
рекуррент |
||||||||||||
ному |
уравнению |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
Ѳ/+і=аѲ<+ 6 е (* + 1). Ѳо = |
0, |
|
|
(14.110) |
||||||
где а и |
Ь—матрицы |
размерностей |
(k X k) |
и (k X k) и |
е(/) — |
||||||||
последовательность |
независимых гауссовских векторов |
в(/) = |
|||||||||||
= (8j (/), |
.. ., |
гк(/)) |
с независимыми |
компонентами, |
Мб; (/) = 0, |
||||||||
Мву (/) == 1, / = |
1, . . . . |
А, |
/ = |
0,1, ... |
|
|
|
|
|
||||
Если |
блочная матрица |
G2 — {b ab ... |
ak~ xb) |
размерности |
|||||||||
{ky^lk) имеет ранг k, то матрица Г<=М0/0? при |
|||||||||||||
является положительно |
определенной. |
|
|
|
|
||||||||
Д о к а з а т е л ь с т в о . |
Из (14.110) находим |
|
|
|
|||||||||
Гі+, = |
МѲі+іѲ;+1 = |
М[аѲі + |
6 е(/+ 1)][аѲ< + |
&в(/+ l)f = |
|
||||||||
Поэтому |
|
|
|
|
|
= аМѲіѲіѴ + frMe (/ + |
1) в (/ + 1) b . |
||||||
|
Г/+І = |
аГ,а* + |
&&\ |
Г0 = |
0. |
|
|
(14.111) |
|||||
|
|
|
|
|
|||||||||
Следовательно, |
|
|
|
|
|
|
|
|
|
|
|||
|
|
Г, = ЬЬ\ |
Г2 = |
М* + |
abb*а"........ |
|
|
|
|
||||
|
|
Гt = bb' + abb'a + ... |
+ а*-'ЬЬг {аГ)*~\ |
|
|
562 ПРИМЕНЕНИЯ К ЗАДАЧАМ СТАТИСТИКИ [ГЛ. 14
Пусть t — k. Тогда, |
очевидно, |
Yt — GiGl и при t > k |
|||||||
|
|
|
г* - |
g gi |
t - |
1 |
|
||
|
|
|
2 |
a’bb" (аУ. |
(14.112) |
||||
|
|
|
|
2 |
+ |
І=к |
|||
|
|
|
|
|
|
|
|
||
Поскольку ранг матрицы G2 по предположению равен к, то |
|||||||||
ранг матрицы |
G2G2 также |
|
равен |
к. Тогда |
из (14.112) следует, |
||||
что при |
к |
матрица Д |
не вырождена. |
|
|||||
Л е м м а |
14.6. |
Пусть |
(Ѳ, І) = |
( [Ѳ,, . . . , |
Ѳ„], [h, ■■■, l N]) ~ |
гауссовский вектор с положительно определенными матрицами *)
соѵ (Ѳ, Ѳ) = |
М[(Ѳ — МѲ)(Ѳ — МѲ)*], |
(14.113) |
соѵ (і, і |б) = |
М [(І — М (і 1Ѳ)) (І — М (II Ѳ)Л. |
(14.114) |
Тогда матрица |
|
|
соѵ (Ѳ, Ѳ ||) = |
М[(Ѳ — М (Ѳ IІ)) (б — М(0| І ) Л |
(14.115) |
также положительно определенная. |
|
|
Д о к а з а т е л ь с т в о . В силу невырожденности |
матриц |
(14,113) и (14.114) у гауссовских |
распределений**) Р (0 ^ а) и |
|||
Р (І <1 b 10 == а) существуют плотности f. (а) |
и f. |
.(b\ а). Отсюда |
||
легко выводится, что существует |
ѳ |
і I |
ѳ |
|
и плотность |
f (b). Поэтому |
|||
из формулы Байеса вытекает, |
что у распределения Р ( Ѳ ^ а ||) |
|||
также существует плотность f |
(а| 6), причем |
|
||
f_ _ (а I b) — |
fu yb\a)Uya) |
|
|
|
|
f-уь) |
|
|
|
'ѳ I5v 17 |
|
|
|
Из факта существования этой (гауссовской) плотности следует,
что отвечающая ей матрица ковариаций соѵ (Ѳ, Ѳ11) не выро ждена, а следовательно, является положительно определенной.
Л е м м а 14.7. Пусть yüt, t — О, I, . . . , — решение уравнения
Y«+i = ayta' + bb* — aytA* {ВВ* + AytA*)+ Ayta* (14.116)
с начальным условием у® = 0 (0 — нулевая матрица порядка
(ky^k)). Если матрица ВВ* положительно определенная, а ранг матрицы G2 равен k, то матрица у°к положительно определена.
Д о к а з а т е л ь с т в о . Пусть Ѳ®, t — 0, 1, . . . , является реше
нием уравнения |
|
Ѳ/+і = aQt |
Ьг{(/ + |
1) |
(14.117) |
|
|
|
|
||||
_*)_По теореме о |
нормальной |
корреляции (теорема 13.1) |
матрицы |
|||
соѵ (1> |
І I Ѳ), соѵ (Ѳ, |
О 11) |
не зависят |
от 0 и | соответственно. |
|
|
**) |
Запись {Ѳ |
а] обозначает событие {Ѳj |
at.........Ѳп^ап}- |
|