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

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

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

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

Добавлен: 10.04.2024

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

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

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

558 ПРИМЕНЕНИЯ К ЗАДАЧАМ СТАТИСТИКИ [ГЛ. И

Факт существования такого предела имеет важное значение для приложений, поскольку в этом случае оптимальная в средне­

квадратическом

смысле оценка

т„

 

О, «отслеживает» вели­

чины Ѳ*,

f ^ O ,

с конечной

ошибкой

даже и тогда,

когда

 

 

k

 

 

 

 

 

 

 

 

 

2 М Ѳ ^ (/)->

со,

/ - >

оо.

 

 

 

 

/=і

 

 

 

 

 

 

 

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

щих существование предела

у ° = ІfітОоу„ заметим,

что вместо

системы

(14.98)

достаточно

рассматривать

систему уравнений

 

 

Ѳ/+і =

аб, +

Ьгі (t +

1),

 

(14.100)

 

 

\t+ 1 =

Aüt +

Be2(t +

1)

 

 

 

 

 

с Ѳ0 — Ѳ0,

|о —■| 0,

 

 

 

 

 

 

 

a = al — ( b o B ) ( B o B ) + A l,

 

 

 

Л =

Л„

(14.101)

b = [(b°b) — (b°B)(BoB)+(b°B)*\>/2,

B = (B°B)'ß, (14.102)

поскольку уравнения для yt как в случае (14.98), так и в слу­ чае (14.100) будут совпадать.

Действительно, если

тѳ (t + 1, /) =

М(Ѳ,+ 1|

О , ) ,

то

 

Y/+i = М [(0,+, — mt+l)(Qt+ l — mt+lY] =

 

 

 

 

 

 

=

М [(Öt+i — m^{t +

1,t) + më<( /+ 1 ,t) — mi+|j X

 

 

X (0<+,

m bt (( +

1, () +

(t +

1,0 —

 

+

~

 

 

=

M [(0t+1 -

(/ +

1, 0) (0t+1 -

möf ( t + l ,

/))*] +

 

 

 

+

M [(V i^(/+

1, () — w ,+1)(m ö^

+

1, () — «;_,)*].

В силу (13.91)

 

 

 

 

 

 

 

 

 

 

м [(Ѳі+і

!> 0) (ѳ/+і —

1, f))] =

Y(*+ 1, 0 =

 

 

 

 

 

=

b о b — (6 о ß ) о ß) + (b о В ) "

=

bb',

а из определения шб ((-J- 1, ^) в силу замечания

к теореме

13.4

следует,

что

 

 

 

 

 

 

 

 

 

 

Щ, V +

1’ t). =

« А + a2l t + (b о В) (В о ß )+ (І<+І -

Л.Ѳ, -

A2l ty

Поскольку mt+x — М

 

(t +

1 , ( ) | ^ +1J, то из рекуррентного

уравнения для

/ и ^( ( +1, / )

получаем

 

 

 

 

 

 

mt+1=

(А ^ + 1) + Й2І* +

 

 

 

 

 

 

 

 

+ (öoß) ( ß ° ß ) ' (|/+ 1_

(/, ( + 1) — Л2|Д


§ 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. Доказательству этой теоремы предпошлем ряд вспомо­ гательных утверждений.


560

ПРИМЕНЕНИЯ

К ЗАДАЧАМ СТАТИСТИКИ

 

[ГЛ. 14

Л е м м а

14.4. Пусть

D

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

(I X k)

и {k X k) соответственно,

и

пусть

 

 

 

 

 

 

D

 

 

 

 

D„

 

Dd

n ^ k

 

 

 

 

 

 

 

 

 

 

Ddn~ ]

 

 

 

блочные

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

(nl X k).

 

либо

Тогда матрицы DID*

и D*nDn,

п > к, одновременно

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

 

блочных

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

Из правила перемножения

матриц вытекает, что

 

 

 

 

 

 

DlDn =

DlDk +

у D»DdK

(14.105)

 

 

 

/<=й

 

 

Отсюда видно, что вырожденность матрицы D*nDn влечет за

собой вырожденность

матрицы

DlDk.

DkDk. Покажем,

что

Пусть теперь вырождена

матрица

тогда вырождены и матрицы D*nDn,

п > к.

 

 

не

Обозначим

х ~ (у,

. . . , xk)

некоторый вектор-столбец,

равный тождественно нулю и такой, что

 

 

 

 

 

 

 

 

x*D*kDkx = 0.

 

 

(14.106)

Установим, что в этом случае

Dd'x =

0 для

всех / ^

к.

Поскольку

 

 

 

k - \

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

DlDk = ^ ( d y

D'Dd

 

 

 

 

 

 

 

 

і= о

 

 

 

 

 

 

 

то в силу (14.106)

 

 

 

 

 

 

 

 

 

 

 

=

0,

Ddx = 0, . . . ,

Ddk~[x = 0.

(14.107)

Положим

 

 

 

 

 

 

 

 

 

 

 

Уо=

X,

уі =

dx = dy0,

уj+\ == dyj,

/ < / г — 1.

 

Тогда

Dy0= 0, Dy{= 0, . . . ,

Dyk

= 0.

(14.108)

 

Но система

векторов

(у0, у {,

. . . ,

 

yk),

каждый

из которых

имеет размерность k, линейно зависима. Поэтому найдутся числа с0> .. •, ск, не все равные нулю, такие, что

к

, 2 ctyi = 0.

( 1 4 . 1 0 9 )


 

 

АСИМПТОТИЧЕСКИЕ СВОЙСТВА

ФИЛЬТРА

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.........Ѳп^ап}-