книги из ГПНТБ / Липцер Р.Ш. Статистика случайных процессов. Нелинейная фильтрация и смежные вопросы
.pdf§ 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
652 |
ОЦЕНКА |
ПАРАМЕТРОВ И |
РАЗЛИЧЕНИЕ |
ГИПОТЕЗ |
[ГЛ 17 |
||||||
Согласно (17.25) |
оценка |
максимального правдоподобия |
|||||||||
|
|
|
|
|
1 |
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 |
|
|
|
|
|
|
2А |
|
|
|
1/2 |
|
|
|
||||||
|
|
|
(Л - |
0) {е2ХТ - |
|
1) + 2А |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
3 месте с (17.53) |
это приводит к |
следующему |
|
представлению: |
|||||||||||||||||
|
|
|
|
|
|
|
я-ѳ |
|
|
|
|
2А |
|
|
|
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). |
||||||||||||||||||
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 !РТ.
