Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Липцер Р.Ш. Статистика случайных процессов. Нелинейная фильтрация и смежные вопросы

.pdf
Скачиваний:
44
Добавлен:
27.10.2023
Размер:
20.66 Mб
Скачать

§ 2]

ПРОЦЕСС ДИФФУЗИОННОГО ТИПА

64?

 

 

2.

Предварительно установим справедливость

следующи

двух лемм.

 

 

 

 

 

 

 

 

 

Л е мм а

17.1.

Пусть

б = 6(х)

—измеримая

функция

с

sup Меб4(Ю< оо для

любых

Ѳь

Ѳ2

(— оо < Ѳ, < Ѳо < оо)

 

ѳ,<ѳ<ѳ2

 

 

 

 

 

1

1

>

Если

т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,<SöU<02 /

a* ® dt < °°’

— 00 <

Ѳ, < Ѳ2 < оо,

 

(17.30)

 

> функция

Мѳб(|)

дифференцируема

по В и

 

 

 

 

 

 

 

т

 

 

 

 

 

 

^ -М ѳ0(£) =

М0 б(£) \ a

t { l ) d W ,

 

(17.31)

Lо

До к а з а т е л ь с т в о . Обозначим

о

(

 

т

т

 

)

ф(Ѳ, ^ ) = ^ Г - ( ^ ) =

 

ѲJ at (W) d W t - ^

Т

 

ех р

^ ^J

a]{W)dt .

w

I

о

X

‘ '

J

І

 

б

 

 

Функция ф(Ѳ,

W) дифференцируема по Ѳ, и (Р-п. н.)

 

ду (Ѳ, W)

 

г т

 

 

 

 

 

 

 

 

at (W )d W t — B \

а] (W)dt

Ф(Ѳ,

W).

(17.32)

 

<ЗѲ

 

 

 

 

 

 

 

 

 

 

 

Пусть — оо <

Ѳ, <

Ѳ2 <

оо. Тогда в силу (17.32)

 

 

М ѳ 26 ( I ) —

М ѳ , 6 ( I ) =

M è ( W ) [ ф (Ѳ 2, Г ) - ф ( Ѳ ; , W ) ] =

 

 

 

 

 

 

ѳ2~ г

 

 

 

 

 

 

 

=

Мб (W) J

^ at (W) — B ^ a](W)dt

ф(Ѳ,

W)dB.

 

 

 

 

Ѳі

L о

 

о

 

 

 

Заметим, что согласно предположениям леммы

 

 

ѳ2

т

 

 

 

г

 

 

 

 

 

М б(ІГ)

at (W )dW t — d j a){W)dt

Ф(Ѳ, W ) d B ^

 

 

ѳJ,

L.оJ/ м .

0(g)

J at (о& d l t - B

J a](l)dt

dB =

 

U2

*

 

 

 

 

о

 

 

 

 

 

 

02

Г

Т

 

1/2

 

= |М Ѳö (|)J at (l) d W t

d B ^ j

Мѳ02(£)Ме{ а 2(|)Л

dB < оо.

о* L

648 ОЦЕНКА ПАРАМЕТРОВ И РАЗЛИЧЕНИЕ ГИПОТЕЗ [ГЛ. 1?

Поэтому

по теореме Фубини

 

 

 

 

0, г г

т

 

 

Ф (Ѳ,

W) dB =

М б (Г )|

j a t ( W ) d W t — b j

a){W)dt

Ѳ,

L0

 

 

 

т

 

 

 

 

 

 

 

 

 

 

= j M eö(i) l a t {l ) dl t - % \ a ] ( l ) d t

dB =

 

ѳ,

Lo

 

 

о

-

 

 

 

 

 

 

Öj

6 ( l ) \ a t { l)d W t dB,

 

 

 

 

=

J Mg

и, значит,

ѳ. Г

/

г

 

\

 

 

 

МѲаб (|) -

М«,б ( 6 ) = J

Мѳ

б Ш J а, (I) d W t

dB. (17.33)

Отсюда следует, что Мѳб(|) является абсолютно непрерывной функцией. Покажем теперь, что в (17.33) подынтегральная функция

Мѳб(6) J а,&) d W t = Мѳб(і)

j at (t) d l t - e j a * ß ) dt

о

 

является непрерывной по Ѳ.

 

Обозначим

 

т

б, (I) = б (l) J a? (g) dt.

б, (6) = б (g) J a, (I) d l t,

Тогда

 

Мѳб(|) j at (l) d W t =

Мѳб, (I) — ѲМѳб2 (l),

 

и для доказательства непрерывности достаточно лишь устанО' вить, что

sup Мѳбг(£)<оо,

і = 1 , 2 ,

(17.34)

й,<Ѳ<Ѳг

для любых Ѳ, < Ѳ2. Действительно, при выполнении этих условий функции МѲ6Д£), г — 1, 2, как было показано, будут абсолютно непрерывными, а следовательно, и непрерывными. Имеем

4 \ 1/2

м0а д < маб д і)м ѳ f at ß ) d l t

Lo

§ 2]

ПРОЦЕСС ДИФФУЗИОННОГО ТИПА

649

 

где в силу неравенства Гёльдера (р = 4, q = 4/3)

М / ' jт

. 0

at (l)d lt

< 2 3

1 4

II

Мѳ ( J

Г

CD

4

I flt(l)

Т

-,4

 

+ ѲJ a?(i)dt

 

<

0

о

 

 

( l ) ^

+ 0 * M e( j a ? ( g ) Ä

<

 

< 2 3 36Г j Mѳа4(|)Л + Ѳ4Г3 j Mfl8(ÖÄ

(17.35)

(Здесь использована оценка

 

 

 

 

 

Мѳ( J а, (I) dlF Л < 3 6 r j

Мѳа}(&)Л,

 

доказанная

в

лемме 4.12.)

Из

(17.35)

и (17.30) следует тре­

буемая оценка

(17.34) с і = 1 .

Аналогично устанавливается

оценка (17.34) и с і — 2.

 

 

 

 

Л е м м а

17.2. Пусть 6(a) — $ тизмеримая функция и

 

 

sup

Мѳб8(!)< со

(17.36)

 

 

Ѳ і < Ѳ < Ѳ 3

 

 

 

для любых

0j <

Ѳ2. Если

 

 

 

 

sup

f ai6(g) dt < oo,

ѳ,<ѳ<ѳ;

J

e r e s i f l .

то функция Мѳ6(£) дважды дифференциируема по Ѳ и

d>М ѳ б ( | )

г / Т

 

Ч 2Т

Мѳ6(|)

^ (D d rJ -

J a\{\)dt

de2

L ѴО

J

о

 

(17.37)

(17.38)

Д

з а т е л ь с т в о . В

силу (17.31) и определения функ­

ций 6, чо/,

62(g) (см. доказательство леммы 17.1)

 

 

^ Ме6(і) =

м ѳ0,(I) - ѲМѳ62 (ё).

(17.39)

Поэтому для существования второй производной -^дгМѳ6(|)

достаточно в силу леммы 17.1 проверить, что

 

sup М ѳбг(|)<°°, і = 1 , 2,

(17.40)

ѳ,<ѳ<ѳ2

 

для любых Oj < Ѳ2.

* / г 2 1 Р . ІИ . Л и п ц е р , А . Ң , Ш и р я е в

650

ОЦЕНКА ПАРАМЕТРОВ И РАЗЛИЧЕНИЕ ГИПОТЕЗ

[ГЛ. 17

 

В силу неравенства

Коши — Буняковского

 

 

 

 

Т

 

ч 8 - , 1/2

 

 

M0öt (І) =

Mö68(|) МѲІ j

a t ( l ) d l t

 

Используя неравенство

Гёльдера (р =

8,

q — 8/7) и лемму 4.12,

находим, что

г т

 

 

 

 

 

 

 

 

Мѳ

j at (l)d l

Мр

J а( (£)<ЛР, + ѲJ

flf(g) dt

 

 

 

о

о

 

 

 

< 2 7

мѳу

ч-

 

 

 

 

 

г

 

г

 

(17.41)

 

< 2 7

284Г3 J Ма8 (I) dt + Ѳ8Г J Мѳа'6 (|) dt

Аналогично показывается, что

Тч 8 - і 1 / 2

MuflJdX Мѳ68(|) Мѳ| J а] (§) dt

1/2

Мѳб8( і ) г / M,a?($)dt

Из этих неравенств и предположений леммы получаем тре­ буемые неравенства (17.40). Для завершения доказательства осталось лишь заметить, что формула (17.38) следует из (17.39)

и (17.31).

3. Д о к а з а т е л ь с т в о т е о р е м ы 17.2. В силу (17.22) и

(17.25)

т

Ѳг(Ю = Ѳ+ -%----------

(17-42)

J a](l)dt

 

о

 

 

Поэтому смещение

т

 

 

 

 

J at (l)dWt

 

М Ѳ ) = М ѳ[ѲИІ)~Ѳ] =

Мѳ о т

(17.43)

j а] (І ) dt

о

§ 31

ГАУССОВСКИЙ МАРКОВСКИЙ ПРОЦЕСС

651

По предположениям теоремы и в силу формулы (17.31)

г

J at (l)dWt

Me f

= Ж Me

I a\(l)dt

о

что вместе с (17.43) доказывает представление (17.28). Далее, из (17.42) получаем

В т(Ѳ) = Ме [Ѳг (|) — Ѳ]2 =

'

т

т2

Мѳ J

at (l) d W t

 

 

 

 

 

 

Но по лемме

17.2

 

 

 

Г T

-| -

2

 

 

d 2 M 0 J

4(1) dt

 

 

 

 

. 0

d&2

 

 

 

 

 

 

 

 

 

=

МѲ

 

 

 

 

что эквивалентно (17.29). Теорема доказана.

З а ме ч а ни е . Более детальное исследование величин Ьт(Ѳ) и BT (Q) для случая, когда at {x) — x t, проводится в следующем

параграфе.

§3. Оценка параметра коэффициента сноса для одномерного гауссовского процесса

1.Будем предполагать, что наблюдаемый процесс £ = (!„ ^ t),

имеет дифференциал

d%t Ѳ|/ dt 4- d W t, |o = °

(17.44)

(cp. c (17.22)), где Ѳ— неизвестный параметр, — oo < Ѳ< °o.

V221*

652

ОЦЕНКА

ПАРАМЕТРОВ И

РАЗЛИЧЕНИЕ

ГИПОТЕЗ

[ГЛ 17

Согласно (17.25)

оценка

максимального правдоподобия

 

 

 

 

 

1

ät,

 

Іт-

 

 

 

 

 

 

в, (к

 

 

 

 

(17.45)

 

 

 

 

 

 

 

 

 

 

 

 

 

J

%tdt

 

dt2

J

I ?

 

 

 

 

 

 

 

 

1

 

 

 

 

поскольку в силу формулы Ито

J* \ t d l t =

— Г].

 

Найдем для

 

 

 

 

о

случая

смещение

bT (ß) =

рассматриваемого

= Мѳ(Ѳг ( |) — Ѳ)

и

среднеквадратическую

ошибку

б г (Ѳ) =

==Мѳ[Ѳг (£)— Ѳ]2.

 

 

 

 

 

 

 

 

 

Введем в рассмотрение вспомогательную функцию

 

Pr (Ѳ, а) =

 

 

 

 

2

Ѵ^Ѳ2 + 2а

 

 

1/2

 

 

 

 

 

 

 

 

 

2а + Ѳ) е- Ѵѳ2+2а т+

Ѳ2 +

2а -

Ѳ) еѴд2+2а т

 

Ѳ2 +

 

 

 

 

 

 

 

 

 

 

 

(17.46)

Т е о р е м а

17.3.

Смещение

ЬТ (Ѳ) и среднеквадратическая

ошибка В Т (Ѳ)

задаются формулами

 

 

 

 

М в) =

ОО

 

 

ехру-()Рг (ѳ>a)^da,

 

 

 

 

J

 

 

 

 

 

(17.47)

Вг(Ѳ) =

е х р ( ~ ) |

 

рг (Ѳ, а) d a + J а

{ exp (— -^ ) рт(Ѳ, а)} da.

 

 

 

о

 

 

о

 

 

 

(17.48)

 

 

 

 

 

 

 

 

 

 

 

Д о к а з а т е л ь с т в о . Для отыскания величин Ьт (Ѳ) и В т (Ѳ)

воспользуемся

представлениями

(17.28), (17.29),

полученными

в теореме 17.2. Предварительно проверим выполнение пред­ положений этой теоремы.

Процесс ! = (!„ Ѳ~і), 0

с

дифференциалом (17.44)

является гауссовским с Мѳ£, = 0

и

дисперсией Г, (Ѳ) = Мѳ|*,

удовлетворяющей уравнению (см. теорему 15.1)

==2ѲГДѲ) +

1,

Г0(Ѳ) = 0.

Отсюда находим

 

 

Г/ (Ѳ)

— 1)>

что влечет за собой условие (17.26) теоремы 17.2.

§ 3]

ГАУССОВСКИЙ

МАРКОВСКИЙ

ПРОЦЕСС

653

 

 

 

 

 

 

Для проверки условия (17.27)

и вычисления математических

ожиданий Мр

Т

- I

г Т

-,-2

 

 

[ l ] d t

, Мѳ

\ l ) d t

используемых при оты­

скании Ьт(Ѳ) и В т(Ѳ), поступим следующим образом.

 

Пусть а > 0 и

 

 

 

г

 

 

(Ѳ,

а) = Мѳехр

 

 

 

— a

j l 2( dt

(17.49)

Если предположить, что

оо

J й‘-% (Ѳ , ß)ä!fl< ОО, ~ о о < Ѳ < о о , 4 = 1, 2, ... . (17.50)

о

Тn - k

то тогда моменты Мр

J

l \ d t

,

4 =

1, 2........

можно

найти,

используя функцию фг (Ѳ,

а), по

формулам

 

 

 

 

 

'

т

-,-k

 

 

 

 

 

 

 

 

 

 

 

М£

Jl \ d t

 

 

( k -

1)1

а*"Чг (Ѳ, а) da.

(17.51)

 

 

Lo

 

 

 

 

 

 

 

 

 

 

 

В самом деле, если

для

некоторого 4 =

1, 2,

... выполнено

условие (17.50), то тогда по теореме Фубини

 

 

оо

 

 

 

оо

 

 

 

 

/

Т

\

 

 

J а4_Іфг.(Ѳ, а) d a =

J

а6-ІМѳехр( — а J

%]dt\da =

 

о

 

Мѳ J

 

о

 

 

 

 

V

о

 

/

 

 

 

=

а6-1 exp

 

а J l]d t^ da =

(k — 1)! Me^J | 2, d/j ,

 

 

 

 

 

 

4 = 1,

2, ...

 

 

 

 

 

Итак, найдем функцию фг (Ѳ,

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

неравенств (17.50).

 

 

 

 

 

 

 

 

 

 

 

2.

Л е м м а

17.3.

Функция

 

 

 

 

 

 

 

 

 

 

Фг(Ѳ,

<з) =

ехр( — - ~ ] р г (Ѳ, а),

 

(17.52)

где рг (Ѳ, а) определено в (17.46).

 

______

 

 

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

Пусть Я = У Ѳ2 -j- 2 а , 0 ^ а < оо. Обо­

значим

и

меры

на (Cr , J r),

отвечающие процессам | ѳ

и | \

имеющим соответственно дифференциалы

 

 

 

 

 

d $ = Q $ d t + d W t,

|° = 0,

 

 

 

 

 

dl) = Xl)di +

d W t,

|£ =

0.

 

 

654 ОЦЕНКА ПАРАМЕТРОВ И РАЗЛИЧЕНИЕ ГИПОТЕЗ [ГЛ. 17

Согласно теореме

7.19

меры

и

А^л

эквивалентны и

 

 

 

 

 

 

 

 

 

 

 

г

 

 

 

 

 

 

 

г

 

 

 

 

-du,ж^- (ifc) =

ехР { (Ѳ — А) JI) d l)

-

 

Ѳ- 2 J(^)2dt

 

 

 

 

 

 

 

 

I

 

 

 

о

 

 

 

 

 

 

 

о

 

 

 

Поэтому

 

 

J

I

 

 

 

Ѳj2J- А 2

 

 

|

 

 

j*J

 

|

 

 

Поскольку

 

а ^

.

 

 

а

 

=

17.53)

 

(Ѳ, а) =

Мѳ exp

 

%] dt

= М exp

 

(^)2 dt

 

 

 

.

 

 

т

 

 

 

 

 

т

 

 

 

 

 

т

ч

 

 

=

М exp

-

а

ф )2 dt +

(Ѳ-А)

I) dl)

-

 

 

 

 

 

(£))* dt

 

 

 

 

 

 

 

 

 

 

<Н----- =—

=

0,

 

 

 

 

 

 

 

 

(17.54)

то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ь

(Ѳ. «) =

 

М exp

 

[Ѳ -

А] j

I) dl)

=

М exp { ^

 

f(^)2 -

Г| } =

 

 

 

 

 

 

 

 

 

=

exp (— = -^ Г) М exp {

 

 

( ^) 2 j .

Случайная

величина

I)

 

имеет

нормальное

 

распределение

ІѴ

( о ,

 

(а—2*-7 l

J)

j, и,

значит (леммаj

1 1 .6 ) ,

 

 

 

 

 

 

 

М « Ф І - ^ ( Г г )2

 

 

 

 

 

 

 

 

 

1/2

 

 

 

 

 

 

(Л -

0) {е2ХТ -

 

1) + 2А

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 месте с (17.53)

это приводит к

следующему

 

представлению:

 

 

 

 

 

 

 

я-ѳ

 

 

 

 

 

 

 

1/2

 

 

(17.55)

 

 

 

Фг(Ѳ,

а ) — е 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(А -

Ѳ) (еш ' -

 

I) + 2А .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где согласно (17.54) А=Ѵ/ 2а + Ѳ2.

После

простых преобразо­

ваний из (17.55) получаем требуемое представление (17.52).

 

З а ме ч а н и е .

Если Ѳ=

0,

а = 1 /2 ,

 

то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т

 

л

 

 

 

 

 

 

 

 

 

 

1

Ф г(°,{

 

Мe x p { - j / r ? Ä j « p r ( 0 , i ) =

/

 

(

 

 

 

 

т

 

(Сс И Г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ет -\-е

 

/Ср. с примером из § 7 гл. 7.)

 

 

 

 

 

 

 

 

 

 

 

 

Завершим доказательство теоремы 17.3. Анализируя пред­

ставление (17.52),

находим,

что неравенства (17.50) выполнены

для любого

6 = 1 , 2,

...

Поэтому

формулы

(17.47)

и (17.48)

ледуют из представлений

 

(17.28), (17.29),

(17.51) и (17.52).

§ 3]

ГАУССОВСКИЙ МАРКОВСКИЙ

ПРОЦЕСС

 

655

3.

Т е о р е м а

17.4.

Оценка

максимального правдоподобия

Ör (É) сильно состоятельна,

т. е. для

каждого Ѳ,

о о < Ѳ <

о о , .

 

Рѳ {Нт ѲГ(І) =

Ѳ} =

1.

(17.56)

 

 

Т-*с

 

 

 

 

 

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

Из (17.49)

получаем

 

 

 

 

 

т

j =

^ ( 0 ,

1),

 

 

Мѳехр j — J f t d t

 

где

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

■Фг(Ѳ,

l) = expj ( - у -

1/2 + 62) Г IX

 

 

 

 

 

x ! ___________________ 2 V 2 + Ѳ2___________________ \ lß

 

{ (V& +

2 -

Ѳ) + (fW + ~2 +

Ѳ) exp ( -

2T рТ + ~ Р ) j

*

Поскольку lim i|)r ( Ѳ , 1) = 0, — о о <

Ѳ <

о о , то

 

 

 

7’ -> oo

 

 

 

 

 

 

 

 

P , ( l « * = “ ) = !■

( 1 7 . 5 7 )

Ясно,

что

 

 

IT lt dWt

 

 

 

 

 

 

(17.58)

 

 

Ѳр ( I ) — ѳ +

 

 

 

J « dt

Поэтому для доказательства (17.56) достаточно показать, что

 

 

I

It dWt

 

 

 

 

lim

 

 

=

1,

— о о < Ѳ < о о .

 

Г •*> ОО

J

e dt

 

 

 

Вытекает это из следующего общего утверждения.

Л е мм а

17.4. Пусть на

некотором вероятностном простран­

стве задан винеровский процесс W =

(Wt, &~t), t ^ O , и случай­

ный процесс

f —

@~t),

0 такой,

что

1) P^J f t d t < ooj = l,

0 < Т <

с»;

2) p [ Jf t d t = o o j = l.

656

 

ОЦЕНКА ПАРАМЕТРОВ

И РАЗЛИЧЕНИЕ ГИПОТЕЗ

[ГЛ. 17

Тогда

случайный процесс

z — (zs,

$ s),

О,

с zs =

J f t d W t,

$ s =

T %s,

где

ts — inf

j*fl

du — s^j

,

является винеровским

и с

вероятностью единица *)

 

 

 

 

 

 

 

 

 

 

 

 

dW„

 

 

 

 

(17.59)

 

 

 

 

lim

-2-

• =

0.

 

 

 

 

 

 

t->ОО

du

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Если ft > 0

 

(Р-п.

н.)

для всех t > О,

то с вероятностью единица xs будет монотонно

возрастающей

непрерывной функцией от s. Отсюда следует,

что случайный

процесс 2^ =

ft dW t

также имеет (Р-п. и.)

непрерывные тра-

ектории.

о

 

 

 

то по свойствам

стохастических

Далее, если s2^ s t,

интегралов

 

 

 

 

 

 

 

 

 

Следовательно, процесс z = (zs, &s),

0, является

квадра­

тично

интегрируемым мартингалом со свойством (17.60). Зна­

чит,

по определению (§ 1 гл. 4) этот процесс является вине­

ровским.

 

 

 

t

 

 

*)

Под J /„ dWa подразумевается непрерывная модификация

стохасти*

 

о

 

 

ческого интеграла, существующая согласно (4.47) и обобщению этого свой­ ства для функций / G !РТ.

Соседние файлы в папке книги из ГПНТБ