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

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

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

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

Добавлен: 10.04.2024

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

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

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

394

ФИЛЬТРАЦИЯ МАРКОВСКИХ ПРОЦЕССОВ

[ГЛ. 9

 

 

Л е м м а

9.5.

Пусть

pß (0) > 0. Тогда процесс paß(s,

t), где

а е £ ,

 

 

допускает стохастический дифференциал

ФРар (s>t) =

ТГШ S

Пу (t) [Рау (S, t) — p„ß (s, 0] dt

(9.64)

 

 

|3

' у <=Е

 

 

и paß(s,

s) = 6 (а,

ß).

 

 

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

В силу условия pß (0) > 0 и леммы 9.3

Р (jtp (і) > 0) — 1.

Поэтому справедлива формула (9.61). При­

меняя к

правой

части

(9.61) формулу Ито и учитывая, что

s),

na(s, t)

и Jtp(0 допускают представления (9.43), (9.54)

и (9.17)

соответственно,

после несложного подсчета приходим

кискомой формуле (9.64).

3.Займемся теперь выводом обратных уравнений интер­ поляции, рассматривая при этом лишь случай конечного мно­ жества Е.

Т е о р е м а 9.5.

Пусть множество Е конечно и ра(0) > 0 для

в с е х а ^ Е . Тогда условные вероятности na{s, t) — Р (Ѳ5 =

а |

s < t, с і ё £, удовлетворяют

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

 

 

 

ds

 

 

а v

 

\

я а (s)

 

 

Яд (s, t)

2 ‘яд (s),

(9.65)

 

 

 

 

 

 

Яд (s)

 

 

 

где

 

 

Яд (S, t)

\ _

у , л

/

^

Яѵ (S, t)

 

 

 

 

 

2

 

 

(9.66)

 

 

Яд (s)

)

~

l i

Is )

я Г ( І Г

 

 

 

 

 

 

 

 

 

ye=E

 

 

 

 

 

 

 

 

 

%*na (s)—

2

 

 

 

 

 

(9.67)

 

 

 

 

 

 

 

у

 

 

 

 

 

 

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

Прежде всего заметим, что

 

 

Яд (s, t) =

М [б (Ѳ „

а) I Г Ц

=

М [М (б (Ѳ8, а) | Р ) ,

Ѳ,) 1

=

 

 

 

=

М[рд0 (s,

t ) \ ^ \ \

=

2

Рау(s,

/)Jty(/).

(9.68)

 

 

 

1

f

 

 

J

 

ye£

 

 

 

Поэтому,

если

установить, что

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

« .(.).

(9.69)

то (9.65)

будет

следовать

из

(9.68).

paß(s, t)

 

 

 

В силу леммы

9.5

у вероятностей

существует про­

изводная

по t\

 

 

 

 

 

 

 

 

 

 

 

 

^Pgß (s, t) _ __I_

<3/

яДО ye£

 

^Yß W л у W [Рау ($> 0 ' Paß (^> 0 ]-

(9.70)


§ 21

ПРЯМЫЕ И ОБРАТНЫЕ УРАВНЕНИЯ

395

 

 

Пусть R (s, t) — II paß (s,

t) II,

 

 

 

Матрица R (s, t) является

фундаментальной: R(s,

s) — единичная матрица, и

 

 

dR (s, t)

R(s,

t)C(t,

со),

(9.71)

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

где C (t,

со) — матрица

с элементами

 

 

 

 

 

 

 

^aa

 

(t)

2

^YO (t) Яѵ

 

 

Пт У, И) =

----------

 

у е £

 

 

 

 

 

5ta (0

 

 

 

 

 

 

>-aß(0 ла (0

 

 

 

 

caß (^. w) —

 

 

 

 

 

 

 

«b(0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

являющимися Р-п. н.

 

непрерывными

функциями,

поскольку

лѵ(0, Лѵа(^) (у, а(=Е)

непрерывны по /, а множество Е конечно.

Если

s < и < t, то

в силу свойств фундаментальных матриц

 

R(s,

t) =

R(s,

и) R (и,

t).

 

Поскольку матрица R(s, и) (Р-п. н.) невырождена,

 

 

R(u,

t) = R~l (s,

и) R(s,

t).

(9.72)

Из (9.71) и очевидного тождества

 

 

 

 

 

° =

-^ -U (s ,

u)R-'(s,

и))

 

следует,

что

 

 

 

 

 

 

 

 

 

 

 

д R 1(s,

и )=

~

С (и,

а>) R

1(s, и).

 

Поэтому

ди

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

___ д_

 

 

 

u)R(s, О —

 

 

 

 

ди R {и, t) — ~ R ~ l {s,

 

 

 

 

 

С (и, со) R

1(s,

и) R (s, t) — С (и,

со) R (и, t)

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

s <

t),

 

 

 

 

 

 

 

R(s, t) = C(s, a)R(s, t).

Покоординатно расписывая эту систему, приходим к системе уравнений (9.69), из которой, как уже отмечалось, вытекают

уравнения (9.65).

Если

в (9.1) коэффициенты At (Qt, I)

не за­

З а м е ч а н и е .

висят от Ѳ*, то рар (5, t) =

Р (Ѳ5

= a IѲ/ =

ß, &~})=P (0s= a 10 f= ß )=

= paß (s, t). При

этом если

множество Е конечно и Pß(0)>0,

ß e £ , то

 

 

 

 

 

дРаЪ(«. *) = Ра (S) Й

Pa ß i s ’ О

Paßt3’ V и Pais)-

(9.73)

ds

 

 

Ра (s)

Ра Н)

 


396

 

 

 

ФИЛЬТРАЦИЯ

МАРКОВСКИХ ПРОЦЕССОВ

 

 

[ГЛ.

9

 

§ 3. Уравнения оптимальной нелинейной экстраполяции

 

 

 

1. Для

 

 

обозначим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Яр(/, s) =

p(0/ =

ß |0 l) ,

 

 

 

 

 

 

 

 

Знание этих вероятностей позволяет решать разнообразные

задачи,

связанные с прогнозированием значений Ѳ( по наблю­

дениям

^ = = { |ц,

 

 

Так,

если

МѲ, < оо,

то

2

ßnß(/,

s)

является

оптимальной (в среднеквадратическом

ß^E

 

оцен­

смысле)

кой Ѳ, по I*.

 

яр(£, s)

 

 

 

 

 

 

 

 

 

 

 

 

 

Для

вероятностей

 

можно

получать уравнения как

по t (при фиксированном s), так и по s (при фиксированном /).

Первые

из

этих

уравнений (которые естественно назвать пря­

мыми уравнениями) дают возможность понять, как ухудшается

прогноз

значений

Ѳг

по

 

когда t растет. Из уравнений

по

s

(t

фиксировано)

можно судить о степени улучшения

прогноза

значений 0* с ростом «числа наблюдений» (т. е. при

s f t).

 

9.1

 

2.

 

Т е о р е м а

9.6.

Пусть

выполнены

условия

леммы

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

(9.2) — (9.6).

 

Тогда для каждого фиксирован­

ного s условные

вероятности {яр(?, s),

t ^ s ,

 

ß е

Е) удовлетво­

ряют (прямым) уравнениям

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Яр(Т s) =

Hp(s)+

I

й*Яр(«, s)du,

 

 

 

(9.74)

где

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

2*яR(u,

s)—

2

XY|3 (ü) яѵ (и,

s)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Система

уравнений

(9.74)

имеет

единственное

решение

классе

 

неотрицательных

непрерывных

решений)

xß(t,

s)

С SUp 2

 

Хп (t, s)

< оо (Р-п. н.).

 

 

 

 

 

 

 

 

 

 

 

ß

 

1

 

 

 

 

 

 

 

 

 

 

 

 

s),

s ^ . t ,

ß е

При фиксированном t условные вероятности {яр(/,

Е} допускают представление

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

яр(*’ s) =

Яр (t, 0 ) + І

Д73(g){>]

pßY(C

и)я ѵ (и)

(y,

g) —

 

 

 

 

 

 

 

О

 

уе^Е

 

 

 

 

 

 

 

 

 

 

 

-

 

2

Аа (у, I) яѵ (и)] } \diu-

 

2 Аа (у,

ю Я ѵ (и) du] .

(9.75)

 

 

Y

 

J

 

 

J L

 

 

Y

 

 

V

 

J

 

 

 

 

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

Для вывода уравнений (9.74)

восполь­

зуемся тем, что при

t ^ s

 

 

 

 

 

 

 

 

 

 

 

 

 

nß(t, s) =

P(0<=

ß |5 r|) =

M[P(0, =

ß |^ 6 ) |^ 6 ] ==

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= М [я„(*)|іГ|]

(9.76)



УРАВНЕНИЯ ЭКСТРАПОЛЯЦИИ

397

и согласно (9.17)

(s) + J* («) du + J Яр (и) - - - | ; (£)Лв (l) dWu. (9.77)

Тогда, беря от обеих частей (9.77) условное математическое ожидание М [ • | #"|], получим

(0 = (s) + М ^ J 2*яр (и) du\&~lj +

Но

+

М fJ яр (и) 4*.ІМ> ~ ^ (6) d W u \ ^ l \

(9.78)

М М 2*яр (и) du 1

=

J

^ A,Yl3 (и) М [яѵ (и) I

du =

 

Vs

 

J

s

у евЕ

 

 

 

t

 

 

t

 

 

= J

Яѵр(и)яѵ(«, s)du — J 2*яp(ü,

s)d«.

(9.79)

5

Y

 

 

 

 

 

Далее, при выводе основной теоремы фильтрации (см. за­ мечание к теореме 8.1) было установлено, что случайный про­ цесс

Au(Р, I) - Au(1)

dWu,

<

7\

Bail)

 

 

 

является квадратично интегрируемым мартингалом. Следова­ тельно,

М ^ J Яр (и) -Л-Ц-(Р-’

dWa

= 0

(Р-п. н.),

что вместе с (9.78), (9.79)

доказывает

справедливость

уравне­

ний (9.74).

лгр (/, s) и Xß(/,

s) — два решения

системы

уравне­

Пусть

ний (9.74).

Тогда

t

 

 

 

 

 

 

 

 

Aß(^,

s) — Aß (t, s ) = J

^yP («) іхУ(u>s) — Xy(u, s)]du

s у e £