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

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

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

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

Добавлен: 10.04.2024

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

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

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

554

 

 

ПРИМЕНЕНИЯ

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

 

[ГЛ.

14

 

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

В силу

(14.74)

 

 

 

 

 

М [(m“+])* S (/ +

1) т “+11ЗГ)и) =

 

 

 

 

 

 

 

=

М {[с (() ut +

а (t) m “+

D ( t ) l { t + 1 )]* X

 

 

 

 

 

X S ( t + 1)[c (t) ut + a(t) m“ + D (t)Ë (t +

1)]|#T } =

 

 

=

(/n“)‘ а (t) S{t +

1) а (t) m“ +

u]c* (t) S (t +

1) c (t) u( -f

 

 

+

2u\c" (t) 5 (t +

1 )a (t) m tu +

 

 

 

 

 

 

 

 

 

+

2M (ë* (t +

1) I r f ) D (t) S (t + 1) (c (t) ut +

a (t) m%) +

 

 

+

M[ë*(/ +

l)D*(t)S(t+

\) D ( t) l( t+ 1)|

 

 

=

 

 

 

=

{т"У a (i)S(t +

1)a(t) m“+

u"tc* (t) S(t +

1 ) c{t)ut +

 

 

 

+ 2u*c(t) S (t +

1) a(t) mut +

Sp S l/2(t +

1) D (t) D" (t) Sp1/2

+

1),

где мы воспользовались тем, что

М (ё(Н - 1 ) | Г И = 0 и

 

 

M [ r ( t + \ ) D ’ (t)S(t+ l)D (t)è(t+

1 )| £гГ] =

 

 

 

 

 

= М[ё*(/+ 1)D* ( t) S ( t+ l ) D ( t) z { t+ \) \ =

 

 

 

 

Sp S 1*2 (t +

1)0(0 Më(f + 1) ë* (^ +

l)/)‘(/)S,/2 ( l + l ) =

 

 

= SpS1/2( /+

l) D {t) D

{t) S lß (t).

 

 

 

(14.76)

 

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

 

6(1), . . . , 6 (Т) — последовательность

 

З а м е ч а н и е .

Пусть

 

независимых гауссовских

векторов (6 (0 = (öj (0 ,

•••, 6ДО)) с не'

зависимыми компонентами, имеющими нулевые

средние и еди­

ничные дисперсии. Рассмотрим процесс

гп,

і =

0, . . . , Т,

опре­

деляемый из рекуррентных соотношений

 

 

 

 

 

 

« 1+1 =

с (0 и, + а(0 гп, + D (06 (t +

1),

 

m0= m, (14.77)

где ut = ut {(.о)

не

зависит

от 6(/ +1) .

Как

и при доказатель­

стве соотношений (14.75), показывается, что

 

 

 

M[mi+iS(^+ l)m <+1| и„ mt\ —

 

 

 

 

 

 

 

 

 

— tnta (0 S {t + \ )a(t)mt +

u"tc* (t)S(t-\- 1) c (t) ut +

 

 

+

2u"tc* (t) S (t + l)a(t)m, +

Sp SI/2(/ + 1) D (t) D (t) S m (t +

1).

 

 

 

 

 

 

 

 

 

 

 

 

 

(14.78)

 

Свяжем теперь с введенными выше матрицами P(t) и функ­

циями p(i),

t =

0........ Т, скалярные функции

 

 

 

 

 

 

 

 

Qt(x) = p(t) + x'P(t)x,

 

 

(14.79)

г д е x < ^ R k. П о с к о л ь к у р ( Т ) = 0, а Р ( Т ) = Н ( Т ) , т о

 

 

 

 

 

 

Q T (х) — X*Н ( Т ) х ,

 

 

 

( 1 4 . 8 0 )


§ 3]

 

 

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

 

 

555

 

Л е м м а

14.3.

Функции

Qt(x), t — О, 1, . ...

Т,

удовлетво­

ряют рекуррентным уравнениям

 

 

 

 

 

 

 

 

 

Ql«

= inf {Х Н (t) X +

V'R (t) V + M [Q/+1 {xf+\)]),

(14.81)

где

 

 

x e

Rk,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x i $ = c(t)V + a{t) X + D(t)l{t + 1).

 

 

(14.82)

При этом inf в (14.81) достигается на r-мерном

векторе

 

 

 

V -=

-

[R (/) +

с* (t) Р (t + 1) с (t)}+с*(t) Р (t +

1) a (t) X.

(14.83)

 

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

Проверим,

что

функции

Q, (х) —

— p(t) + х*Р(t)X

удовлетворяют

уравнению

(14.81),

т.

е.

что

p{t) + х*Р{і)х =

 

 

 

 

 

 

 

 

 

 

 

 

=

inf [хН (t) X +

V'R (0 1/ +

p {t + 1) + M [ ( x ^ y p (t +

1) xf’+\}}.

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

(14.84)

 

Положим

 

 

 

 

 

 

 

 

 

 

 

 

 

/ (V, X) =

V'[R (0 +

c*(t) P ( t + l ) c (f)] V +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ 2VV (t) P (t

1) a (/) x.

(14.85)

Тогда

с

учетом

замечания

к лемме

14.2 находим,

что

(14.84)

эквивалентно уравнению

 

 

 

 

 

 

 

 

 

р (t) + х*Р (t) X =

 

 

 

 

 

 

 

 

 

 

 

 

 

= p ( t+

1) +

Sp PU2(t +

1) D (t ) D

( t ) P ,/2 (t +

1) +

inf / (V, X ) .

Но в силу (14.62)

 

 

 

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P (t) =

P(i +

1) + Sp P 1/2 (/ + 1) D (t) D*(t) P 1/2 (* +

1).

 

 

Поэтому

надо лишь проверить,

что

 

 

 

 

 

 

 

 

 

 

 

 

 

х*Р (t) X =

inf / {V, х)

 

 

 

 

(14.86)

для любого X е

Rk.

 

 

г

 

 

 

 

 

 

 

входящая в І{Ѵ, х),

была бы положи­

 

Если

матрица

R(t),

тельно

определена,

то тогда

J(V, х )> — °о

и inf/(E, х) дости-

гался

бы

на

векторе

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V =

-

{R (t) + с* (t) Р (t + I) с (t)]+ с*(t) P(t +

I) a(t) X

(14.87)

и

непосредственно

легко было

бы

проверить,

что

 

J(V,x) —

х Р (t) X.

 

 

 

равенства (14.86) в общем случае

рас­

 

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

смотрим

систему

алгебраических

уравнений

(относительно

ѵ = (Ѵ1г....

V,))

 

 

 

 

 

 

 

 

 

 

 

 

± V J ( V , х) — О,

( 1 4 . 8 8 )


556

ПРИМЕНЕНИЯ К ЗАДАЧАМ СТАТИСТИКИ

[ГЛ. 14

т. е. систему

 

 

[Я (*) +

с*(/)/>(/+ \)c{t)]V = -c*(t)P {t-\- l)a(t)x.

(14.89)

Согласно лемме 13.3 эта система совместна й вектор V, определяемый формулой (14.83), есть одно из ее решений.

Поэтому

минимум

квадратичной

формы J(V,x) достигается

на векторе V, и для проверки равенства (14.86) осталось лишь

установить, что х*Р (t) х = / (V, х),

т. е. что

x'P{t)x =

x*[H(t) +

a* ( t ) P ( t + \ ) a { t + \ ) —

a*{f)P(t+ \)c(t)(R(t)-\-c*(t)P(t)c(t))+c*(t)P(t+ l)a(t)\x.

(14.90)

Справедливость этого равенства вытекает из определения матриц P(t) (см. уравнение (14.59)).

5.

Возвратимся к доказательству

теоремы 14.2.

 

Рассмотрим

управление й — (й0,

. . . . йт- х), определенное

в (14.64). Тогда в силу леммы 14.3

 

 

 

— М [Q<+1 (mf+1) — Qt (mt)] =

M [m]H (t) mt + ü]R {t) üt\

(14.91)

Суммируя равенство (14.91) по 1 от

0 до

Т — І и учитывая,

что rriQ— m, находим

 

 

 

 

 

 

 

г-і

 

 

 

 

 

Q0 (m) — МQT(mr) + 2

М [т*Я (t) mt +

utR (t) üt] =

 

 

 

 

Г

 

T- 1

 

 

 

 

=

2 Mm.H (t) mt + 2

MutR (t) üt.

(14.92)

 

 

 

/=о

 

(=0

 

 

С

другой

стороны,

пусть

u = (u0, . . . ,

uT-.i) — любое из

управлений, удовлетворяющее условию (14.56). Тогда в силу лемм 14.2 и 14.3

м [Qi+i {mt+1)

Qt {mf)\ <

M

H (0 rnut +u]R {t)

U t],

откуда

вытекает,

что

 

 

 

 

 

 

 

 

 

г

 

 

r-i

 

 

 

 

Qo ( m) <

2

М (m“Y Н (t) m? + 2

M иfR (t) и..

(14.93)

 

 

 

 

<=о

 

t

 

 

Сравнение (14.92) с (14.93) доказывает оптимальность управ­

ления

 

й = (йо,

•••> «Г-і)- Формула

(14.67)

следует из

(14.71),

(14.79)

и того,

что

 

 

 

 

 

 

 

 

У (й) =

Qo (m) + 2

Sp H u\ t ) уtH42 (0.

 

 

 

 

 

 

<=o

 

 

 

 



§ 4]

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

СВОЙСТВА

ФИЛЬТРА

557

З а м е ч а н и е .

Пусть Ѳ0 — m — детерминированный

вектор,

ö(/) = 0.

Рассмотрим

задачу управления

(по полным данным)

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

процессом

Ѳ„ ^ =

0,

Т, с

 

 

Qt+i = c{t)ut + a{t)Qt,

Ѳ0s= m,

(14.94)

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

 

т

Т - 1

 

 

 

 

 

 

 

 

 

 

Ѵ ( и ) = ^ ѳ ;я (t) 0t + S

utH (t) ut.

(14.95)

 

 

 

<=0

f=0

 

 

 

В этом частном случае оптимальное управление

 

S/ =

- [ ^ W +

c*(0/>( f + 1)с(/)]+ с‘ (/)Р (^+ 1)а (/)Ѳ„

(14.96)

где

 

 

 

 

 

 

 

a

4 +І — с(/)й,+

я(/)Ѳ,,

% = m,

 

 

 

V (ü) = m*P (0) m.

 

(14.97)

 

 

 

 

§ 4. Асимптотические свойства оптимального линейного фильтра

1.Рассмотрим задачу фильтрации*) для гауссовского про­

цесса

(ѳ,І) = [(ѳ,(0, . . . . 0*(0 ), (1,(4

 

Ы))\,

t = 0,

1 , —

удовлетворяющего рекуррентным уравнениям

 

 

 

4+1 — «і4 + а21/ + 4 ei {t + 1) +

Ь2е2 (^ ~Ь 1)>

(14.98)

ков (k

I/+1 = П, 4 +

+ В,е14+ 1) +

В2ег {І

+ 1)

 

 

 

 

 

 

 

 

с постоянными матрицами

а,, а2,

Ьи Ь2,

Аи А2,

В1 и В2 поряд­

 

X k), (k X I ) , (k X k),

(k X 4

(/ X /г), (/ X /), (/ X *),

(I X /)

соответственно.

 

 

 

 

 

 

Обозначим т , = М (4|£Г |) и

у, =

М [(4 — т ()(Ѳ( — т ()’].

Тогда согласно теореме 13.4 матрица ошибок у, удовлетворяет рекуррентному уравнению

Y(+i = a,Ytax+ b ° b —

 

 

 

- \b о В +

а.ѵИП lß 0 ß + 44,yHI]+ [*0 B + aiYHif>

(14-" )

где b ° b b\b\ +

6262, b 0 В — b\B\ + 6262, B° В = BiB\-\-B2 B2 .

В этом параграфе будет исследоваться асимптотическое

поведение матриц у* при t->oo.

В предположениях,

сформули­

рованных далее в теореме 14.3,

будет показано, что

lim

уг = у°

 

 

 

f-> оо

 

существует и 0 < S p y ° < o o .

') По доводу принятых далее обозначений см. § 2 гл. 13.