Файл: Медич Д. Статистически оптимальные линейные оценки и управление.pdf

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

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

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

Добавлен: 11.04.2024

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

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

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

Вычитая почленно из уравнения (8-15) уравнение (8-16), получаем:

 

 

£ { * ( М О І * ( М О - 2 * ( ' , і О ] ' } = = 0 -

 

 

 

(8-17)

Точно

таким

же

 

 

образом

получим

 

соотношения

E [x* (t, I t) x'

(t, 101 0

и

E [x* (t,

I t)x*'

(t, \ f)\=

0,

так

что

 

E {x* (tt \t)(x{tt\

 

t) -x*

{t, \ t)}'}

= 0.

 

 

(8-18)

Вычитая

почленно

из уравнения (8-18) уравнение

(8-17),

имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Е {[**(*, 10-

x

 

| 0]

(f, \t)-x*

(t,

1 01} =

0.

 

Ho x*

{t1

\ f)—x(tt\t)

 

=

x (t, 11) ~ Jc*(f, 10.

откуда

сразу

следует утверждение

леммы.

 

 

 

 

 

 

 

 

8-2. ОПТИМАЛЬНАЯ ФИЛЬТРАЦИЯ

 

 

 

 

 

 

 

 

Решим

уравнение

Винера — Хопфа

при

t\ = t,

т,

е.

для задачи

оптимальной

фильтрации,

где

(x(t),

t^U}

гауссовский марковский процесс, описываемый уравне­ нием (7-1). Здесь предполагается в отличие от постанов­

ки задачи § 7-1 и

задачи

оптимальной фильтрации

в § 7-3, что процессы

 

 

 

{w(t),

t^Q

и

{v(t),

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

E[w(:t)v'(x)]=S(t)ö(t-x)

для всех t, x^to,

где S(t)—непрерывная

матрица

раз­

мера

рХт.

 

 

 

 

Уравнение Винера — Хопфа

для, оптимальной

филь­

трации

в случае,

когда S(t)=0,

было

решено Калманом

и Бьюси (Л. 8-1]. Решение в настоящей главе отличается от их решения только несущественными деталями.

Алгоритм оптимальной фильтрации, излагаемый ни­ же, обычно называют фильтром Калмана — Бьюси.

Фильтр

 

 

 

 

Рассмотрим уравнение

Винера —Хопфа

на

интер­

вале tQ^a<t, где ti = t:

t

 

 

 

 

 

 

 

E [x (0 z'(p)]-JA

{t,

z) E[z(z) z' (o)] dx =

0.

(8-19)

22—85

 

 

 

337


 

Важно заметить, что здесь рассматривается

изме­

нение а только на

полуинтервале

вида

[to, t).

Случай

a=t

будет рассматриваться

отдельно по

причинам, яс­

ным из дальнейшего.

 

 

 

 

(8-19) по t,

 

Дифференцируя

левую

часть уравнения

получаем два слагаемых:

 

 

 

 

 

 

 

^fE[x(t)z'

(a)] = E[x{t)z'

(а)] = Е {[F (t) X {t)

+

+

G (t) w (01 z' Щ =F(t)E

[x (f) z' (,)] + G(t)E

[w (t) z' (,)];

 

 

t

 

 

 

 

 

 

(8-20)

 

 

T)£[z(x)z'(3 )]rfx =

 

 

 

 

~^А(І,

 

 

 

t

 

 

 

 

 

 

 

 

= = |ІеЬ _ 1я [2 (х)г'(с»)] Л + Л(/, t)E[z(t)z'(a)]t

 

(8-21)

t.

 

 

 

 

 

 

 

 

где

г 0 < а < г .

 

 

 

 

 

 

 

 

Рассмотрим последнее

слагаемое в уравнении (8-20).

Из (7-2) имеем:

 

 

 

 

 

 

 

так что

z' (а) =х'{а)Н'{а)

+ ѵ' (о),

 

 

(8-22)

 

 

 

 

 

 

 

 

 

Е [w

 

{t)z'{e)]=E{w(t)x'{o)}H'(a)+E{w{t)v'{o)l

Поскольку

a<t,

то E[w(t)v'(а)]

= 0. Поэтому

 

 

E[w{t)z'{a)\=E{w{t)x'{o)\H'{o).

 

 

 

(8-23)

Решение уравнения (7-1)

можно представить

в виде

 

X (а) =

Ф (a, Q X (/„) + [ Ф (а,

x) G (x) w (х) Л ,

(8-24)

 

 

 

 

І

 

 

 

 

 

где

Ф'(а, т)—переходная

матрица

состояния.

 

Используя этот результат,

получаем:

 

 

 

 

Е [w {t) X' (о)} = E[W (0 x' (t0)} Ф' (a,

Q

+

 

 

 

а

 

 

 

 

 

 

 

 

+ jEV)ш'(х)1G'С")ф'(°« т ) d % -

 

Первое слагаемое в правой части этого соотношения

обращается в нуль,

так как x(t0)

не зависит от

{w(t),

338


t^t0}. Во втором слагаемом Е [w (t) w'(х)] = Q (t) à (t—т.). Однако поскольку t находится вне интервала интегри­ рования, это слагаемое также обращается в нуль. Сле­ довательно,

для

 

tQ^a<t.

E[w{t)x'(ü)]

= 0

 

 

 

 

 

(8-25)

 

 

 

 

 

 

 

 

(8-23), по­

 

Подставляя этот результат в уравнение

лучаем:

 

Я ( ш ( 0 г , ( о ) 1 = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

для

t0^a<t.

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

последнее слагаемое

в вы­

ражении (8-20) равно нулю, и

 

 

 

 

 

 

 

 

 

 

-^Е[х

(0 г'

(,)] = F(t)E

[x (t) z'

(з)]

 

 

(8-26)

для

t0<^a<t.

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим второе слагаемое в правой части урав­

нения

(8-21). Используя уравнение (7-2),

получаем:

 

 

Е [z (0 z' (,)] = £ {[Я (0 x(t)

+

o (01 z' Щ

 

=

 

 

 

 

= Я (0 Я (0 г' (?)] - f Е [и (0 z'

(з)].

 

(8-27)

 

Из

уравнения (8-22) следует, что

 

 

 

 

 

 

Е [V (0 z' (з)] =

Е[ѵ (0 x' (-,)] Я ' (о) +

£ [у (0 о' (01.

(8-28)

 

Так

как з < < / ,

то Е [ѵ (t) ѵ' (з)] =

0.

 

 

 

 

 

 

Используя уравнение (8-24), получаем:

 

 

 

 

 

 

Е [V (t) X' (з)] =

£ [о (0 X ' (f.)] Ф' (a,

f0) - f

 

 

 

 

 

а

 

 

 

 

 

 

 

 

 

 

 

 

+

f £ [у (0 te»' (x)] G' (т) Ф' (з,

т) d

t .

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

Первое слагаемое в правой части этого соотношения

обращается

в нуль, поскольку

x(U)

не зависит

от

{v(t),

i^to).

Хотя

E[u(t)w'(x)]

= S(t)ô(t—т.),

второе

 

слагаемое

также

обращается в нуль, так как t

находится вне ин­

тервала интегрирования. Следовательно,

 

 

 

 

 

 

 

 

E\[v(t)xf(a)]

= 0

 

 

 

 

 

 

для

to^a<t.

В результате уравнение

(8-28)

приводится

к виду

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£ И 0 2 ' ( а ) ] = 0

 

 

 

 

 

 

22*

 

 

 

 

 

 

 

 

 

 

 

339



для to^a<4, а это означает, что уравнение (8-27) при­ нимает вид

 

 

E[z{t)z'{a)]

= H{t)E[x{t)z'{o)}

 

 

(8-29)

для

to<^a<t.

 

 

 

 

 

 

 

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

уравнение

(8-21)

можно

записать

в виде

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

-§f^A{t,

 

х ) £ [ г ( х ) г ' ( з ) ] ^

=

 

 

 

t

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= Г дА (t,

x) Е [ z ( т ) z , ( з

) ]

^ 0 / /

( 0 £ ( ѵ ( 0

г ,

( 8

3 0 )

 

Объединяя уравнения (8-26) и (8-30), получаем част­

ную производную по t

от (8-19):

 

 

 

 

 

 

 

 

t

 

 

 

 

 

F (t) E [x (0 г' (a)] -

j"

£ [г (x) г' (,)]

rfx

-

 

 

— A(t,

 

t)H(t)E[x(t)z'(•>)]=()

 

 

для

tn<3<^t.

Подставляя сюда

E [x (t) z' (т)} ИЗ уравне­

ния

(8-19),

получаем:

 

 

 

 

 

 

t

[F(t)~A{t,

t)H{t)]^A{t,

x) E [z (x) Z' ( I ) ] rfx -

Ç[F(0i4(/,

x)

дА [t, x)

 

 

 

 

 

или, что то же самое,

 

 

 

 

— A(t, t)H(t)A{t,

4)]E[z(%)z'(i)]dt

= 0

(8-31)

для t0^a<t.

Для того чтобы уравнение (8-31) удовлетворялось, Достаточно, чтобы матрица A (і, т) удовлетворяла урав­ нению в частных производных вида

F (t) A (t, x) - д Л § т ) - A (t, t) H (t) A (t, x) = 0 (8-32)

340