Лекция 9
Браунов мост

В тази тема:

Нашата цел е да покажем, че много асимптотични резултати в непараметричната статистика следват от свойствата на брауновия мост.

9.1  Браунов мост

9.1.1  Определния

Нека с W(t), (0 Ј t) означим стандартния винеров процес. Това е гаусов процес с независими нараствания, E W(t) = 0, D W(t) = t. Лесно се извеждат крайномерните разпределения на този процес - всички те са гаусови. Съгласно теоремата на Колмогоров съществува вероятностна мярка в пространството RT от всички функции определени на интервала [0,T]. Мярката на подпространството от непрекъснати функции C[0,Т] е 1, така че траекторите на процеса лежат в него с вероятност 1, Затова обикновено се казва, че той е процес с непрекъснати траектории.

Определение 1 Браунов мост наричаме процеса: w(t) = W(t)-W(1)*t, 0 Ј t Ј 1.

Определение 2 Да означим с Fn(t) извадъчната функция на разпределение, построена по n независими наблюдения над равномерно разпределена сл.в. в интервала (0,1). Извадъчен браунов мост наричаме процеса: wn(t) = Цn(Fn(t)-t), 0 Ј t Ј 1.

Траекториите на този процес естествено са непрекъснати от ляво и лежат в пространството D[0,1] на непрекъснати от ляво и ограничени функции с краен брой точки на прекъсване.

9.1.2  Конструкция на извадъчен Браунов мост

Да разгледаме едномерния случай.

Да означим с b(k,p) следната функция на векторите k, p О Rs (n = еi = 1s ki, ki і 0, ki - цели числа, еi = 1s pi = 1, pi > 0):

b(k,p) = ж
з
и
n
k1,k2,...,ks
ц
ч
ш
p1k1p2k2... psks = n! s
Х
i = 1 
p1ki
ki!
(9.1)

Нека t0 < t1 < t2 < ... < ts са фиксирани числа в интервала (0,1), t0 = 0, ts = 1. Да разгледаме съвместното разпределение на бройките попадения в интервалите Di = [ti-1,ti) на n независими еднакво разпределени равномерни сл.в. Това, разбира се, е добре известното полиномно разпределение, което се задава с формулата (9.1). За разпределенията на извадъчната функция на разпределение можем да използуваме този резултат по следния начин:

Нека сега t0 < t1 < t2 < ... < ts, i = 1,2,...,s са фиксирани числа в интервала [0,T], t0 = 0, ts = T. Да разгледаме поасонов случаен процес p(t), p(0) = 0 с интензивност l зададен на реалната полуправа - t О R+. Ще го разгледаме в интервала [0,T) Неговите съвместни разпределения в точките t0 = 0 < t1 < t2 < , ..., < ts = T се определят от (съвместните разпределения на) независимите нараствания на процеса в интервалите: Di = p(ti)-p(ti-1). Имаме, че li = l(ti - ti-1) и, следователно:

P (
З
i 
{Di = ki}) = s
Х
i = 1 
e-li piki
ki!
= ln e-lT s
Х
i = 1 
p1ki
ki!
= ln e-lT
n!
b(k,p),
(9.2)
където сме означили pi = (ti-ti-1)/T . Условното разпределение на тези нараствания при условието: p(T) = n се получава същото полиномно (достатъчно е да разделим формула (9.2) на условната вероятност

P (p(T = n)) = ln e-lT
n!
.

Нека p(u), u О [0,1] е еднороден поасонов процес с параметър l. Предполага се, че траекториите му са непрекъснати отляво. Да разгледаме трансформирания поасонов процес h(t) = p(F(t)). Така получаваме следната теорема:

Теорема 1 Ако ф.р. F е непрекъсната, то разпределението на процеса n Fn(t) съвпада с условното разпределение на поасоновия процес h(t) = p(F(t)) при условието h(Ґ) = n, (p(1) = n).

Доказателство: Наистина, ако F(t) = t (t О (0,1)), теоремата е доказана. Нека разгледаме произволна непрекъсната (без скокове) F и x е такава, че P (x < x) = F(x). Тогава разпределението на сл.в. h = F(x) е равномерно в интервала (0,1). Q.E.D.

Следователно, за непрекъснатите ф.р. можем да се съсредоточим върху разглеждането на тези крайно-мерни разпределения.

9.1.3  Крайно-мерни разпределения

Теорема 2 Крайномерните разпределения на wn клонят към крайномерните разпределения на w.

Доказателство: Доказателството е почти очевидно. Да означим с Dk = tk-tk-1, U = еsk = 1 ukDk и разгледаме лог-характеристичната функция на нарастванията wn:

yn(u) = lnE exp(i s
е
k = 1 
uk(wn(tk)-wn(tk-1) ) =
lnE exp(i s
е
k = 1 
uk (p(tk)-p(tk-1)-nDj)/Цn ) =
-i UЦn+n ln( s
е
k = 1 
Dk ei uk/Цn) =
-i UЦn+n ln(1+ s
е
k = 1 
Dk(i uk/Цn-uk2/n +O(n-3/2))) =
-i UЦn+n (i U/Цn+ U2
2n
-
е
u2kDk

2n
)+O(n-1/2) =
- 1
2
( е
u2kDk-U2) +O(n-1/2).
Сега остава да покажем, че точно това е х.ф. на w.
s
е
k = 1 
uk(wn(tk)-wn(tk-1)) =
s
е
k = 1 
uk(W(tk)-W(tk-1)- W(1)Dk)) = s
е
k = 1 
(uk-U)(W(tk)-W(tk-1)
E exp(i s
е
k = 1 
uk(W(tk)-W(tk-1) ) =
exp(- 1
2
s
е
k = 1 
(uk-U)2Dk) = exp(- 1
2
( s
е
k = 1 
uk2Dk - U2)).
Q.E.D.

От тази теорема, разбира се, следва сходимост по разпределение на произволни ограничени функции зависещи от краен брой стойности на процеса в фиксирани точки по времето. За по-общи функционали това не е достатъчно.

Определение 3 Модул на непрекъснатост CD(y) О D[0,1] наричаме функцията:

CD(y) =
sup
0 Ј t Ј 1 -D 

sup
0 < h < D 
|y(t)-y(t+h)|

Ще докажем едно вспомагателно твърдение за модула на непрекъснатост CD(wn) .

Лема 1


lim
D® 0 

limsup
n®Ґ 
P (CD(wn) > e) = 0

Доказателство: Да разгледаме двоично-рационални D = 2-l, l > 3. Имаме при l < m

CD(wn) Ј C[m]D +2
max
k Ј 2m 
C[m](k) =
(9.3)

sup
|k-j| < D2m 
|wn( k
2m
)-wn( j
2m
)|+2
max
1 Ј k Ј 2m 

sup
[(k)/( 2m)] Ј u Ј [(k+1)/( 2m)] 
|wn( k
2m
)-wn(u)|.
Да оценим сега първото събираемо (|k-j| Ј 2(m-l))
Зmr = l З2rk = 1{|wn( k+1
2r
)-wn( k
2r
)| < e
r2
} М {C[m]D Ј e}
P (C[m]D і e) Ј P (Иmr = l И2rk = 1|wn( k+1
2r
)-wn( k
2r
)| > e
r2
).
За да оценим вероятностите на всичките тези събираеми, ще използуваме неравенството на Чебишов за 4 момент (E (Sn-np)4 Ј np+3n2p2) на Биномно разпределение с вероятност за успех p = 2-r:

wn( k+1
2r
)-wn( k
2r
) = Цn (Fn( k+1
2r
)-Fn( k
2r
) - 1
2r
) = Sn- n 2-r
n1/2
,
P (|Sn- n 2-r| > eЦn
r2
) Ј (np+3n2p2)r8
e4 n2
= r8
e4
(n-12-r+3*2-2r)
m
е
r = l 
r8
e4
(n-1+3*2-r) Ј c
e4
( m9
n
+ l8
2l
)
Като положим m = 3 ln2 n и отчетем произвола при избора на l виждаме, че изразът може да се направи произволно малък.

Сега (9.3) можем да препишем във формата

P (CD(wn) > 3e) Ј P (C[m]D > e)+2m P ( C[m](k) > e).
(9.4)

Да се заемем с втората част на (9.4). Събитието {sup|wn([(k)/( 2m)])-wn(u)| > e} под знака на вероятността означава (при фиксираното m: 2m = n3), че в интервала ([(k)/( n3)],[(k+1)/( n3)]) е изпълнено |n(Fn(u)-u) -Fn([(k)/( n3)])| > Цne. Тъй като при достатъчно голямо n имаме Цne і 3, значи в интервала са се случили поне две събития: Sn і 2. P (Sn і 2) = O(n2 p2) . Тъй като p = n-3 от това следва,че n3 O(n-4) = O(n-1) ® 0. Q.E.D.

Тази лема е основното оръжие при изучаването на близостта на извадъчния Браунов мост към теоретичния. Сега лесно следва теоремата.

Теорема 3 Нека F е функционал определен на D[0,1] и F(wn) и F(w) са сл.в. Нека F е непрекъснат в метриката на C[0,1] в всяка точка на D[0,1] . Тогава F(wn) ® F(w)

Доказателство:  Да означим за всяко x О D[0,1]:

||x|| =
sup
0 Ј t Ј 1 
|x(t)|,   F+e(x) =
sup
||x-y|| < e 
F(y),   F-e(x) =
inf
||x-y|| < e 
F(y).

Ако сега разгледаме непрекъснатата частично линейна функция през точките (kD, x(kD) и я означим с xD(t) , ще видим, че

||x-xD|| < CD(x)
(9.5)

Твърдението на теоремата лесно следва (ref3) от връзките:
(wn - wnD ) (по разпределение и п.с.);
(wnD - wD) - (вж. теорема 9.2);
(w - wD) (по разпределение и п.с.) и разпределенията на стохастично наредените F-e, F , F+e върху тях. Q.E.D.

9.2  Статистики свързани с извадъчното разпределение

Сега можем да сумираме знанията си по непараметрична статистика. Разглеждат се два основни типа статистики свързани с извадъчното разпределение.

Определение 4 Статистика от интегрален тип наричаме:

S(mn) = h( у
х
g(x) dmn(x))

Лесно се вижда, че за изучаването на граничното поведение на тези статистики е достатъчно да се изучава сходимостта по разпределние. При това интегралния (сумиращ) характер на функционала води до асимптотично нормално поведение. Това остава верно и при произволна размерност на елемента на извадката. Ще цитираме една теорема от този тип. Втората и част позволява лесно да се изучава асимптотичното поведение на максимално- правдоподобните оценки и самото правдоподобие в точката на максимума.

Теорема 4 Нека функцията h(x):Rs® R1 е диференцируема в т.a = тg(x) dm(x) О Rs и матрицата от втори моменти: s2 = т(g(x)-a)(g(x)-a)ўdm(x) е крайна.

Вторият тип статистики е от ''точков '' тип. Нека F е функционал определен на множеството от разпределения и той е непрекъснат относно някакво разстояние върху това множество.

Определение 5 Статистика от точков тип наричаме:

S = F(mn)

За изучаването на граничното поведение на общи функционали от този тип вече не стига слабата сходимост. Налага се да се използуват теореми от типа на Гливенко-Кантели и Брауновия мост (виж теорема 9.4).

Двата класа, разбира се, имат непразно сечение.


Начало на лекцията | Съдържание | Индекс


File translated from TEX by TTH, version 2.10.
On 4 Jun 1999, 15:57.