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

Окороков Фракталы в фундаменталной физике.Фракталные свойства множественного образования частиц и топология выборки 2009

.pdf
Скачиваний:
174
Добавлен:
17.08.2013
Размер:
5.73 Mб
Скачать

cov i

, j

, пусть x x1, , xn * Rn , F x n -мерная

 

i, j 1, ,n

 

функция распределения случайного вектора .

Определение П4.3. Случайный вектор n Rn имеет n -мерное гауссовское распределение с параметрами , R , если его ха-

рактеристическая функция имеет следующий вид:

f

 

exp

 

i *

 

 

1

*R

 

 

,

Rn .

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Замечание П4.2. Любая компонента k

гауссовского вектора

имеет распределение N k ; D k , где k

k -й элемент вектора

; D k k -й диагональный элемент матрицы R .

Замечание П4.3. Если матрица R 0,

то есть положительно

определена, то F x имеет плотность распределения p x , где

x Rn , определяемую по следующей формуле:

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

*

1

 

 

 

 

 

 

 

 

 

 

 

 

exp

 

 

 

 

x

 

R

x

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

p x

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 n 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

R

 

0 – определитель ковариационной матрицы.

 

 

 

 

 

 

Ниже представлены важнейшие свойства гауссовских векторов.

1. Если

cov ,

 

*

 

 

* 0, а

вектор

* , * *

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

является гауссовским, то гауссовские векторы

и независимы.

2. Пусть N ;R

и A B, где A Rm n , B Rn ,

тогда

N

 

;R

, где A

 

B, R

AR A*.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3. Пусть n – последовательность гауссовских случайных век-

с.к.

торов. Если имеется среднеквадратичная сходимость n при

n ,

то N ;R , где

 

lim

, R

limcov n , n ,

 

 

 

n

n

n

причем указанные пределы существуют и конечны.

418

4. Если n – последовательность гауссовских случайных век-

п.н.

при n , то

торов и почти наверное имеется сходимость n

с.к.

 

при n .

присутствует среднеквадратичная сходимость n

Важно отметить, что в общем случае данное утверждение не верно.

5. Теорема П4.1 (о нормальной корреляции). Пусть вектор

* , * *

является гауссовским, причем R 0, тогда:

 

а) условное математическое ожидание имеет вид

 

 

 

 

 

 

 

 

 

 

 

1

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б) , то есть

 

и являются независимыми;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

*

в) если ,

то 0, cov , R R R

R ;

г) условное математическое ожидание

 

имеет гауссовское рас-

 

пределение

 

 

 

 

 

 

 

 

1

 

*

 

где введены обозначе-

N ;R , R

R R

R ,

ния R cov , , R

cov , .

 

 

 

Замечание. Данная теорема дает явный вид среднеквадратич-

ной оптимальной оценки для по наблюдениям в

гауссовском случае. Важно отметить, что линейно зависит от . 6. Если , , составляют гауссовский вектор, причем и

некоррелированы, то , .

*

7. Если компоненты вектора 1, , n – гауссовские и не-

зависимые в совокупности, то является гауссовским случайным вектором.

Пример П4.1.

Двумерный

случайный

вектор

1, 2 имеет

нормальное

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

f

t1,t2 име-

ет следующий вид [343]:

 

 

 

 

 

 

 

 

 

1

2

 

2

 

f t1,t2

exp i m1t1

m2t2

 

c11t1 2c12t1t2 c22t

 

,

2

 

 

 

 

 

 

 

 

419

2

причем квадратичная форма Q t1,t2 cijtitj , c12 c21 неотрица-

i, j 1

тельна, то есть Q t1,t2 0 при любых действительных t1,t2 . Если ранг квадратичной формы Q t1,t2 равен двум, то есть детерми-

нант матрицы C cij , i, j 1,2 отличен от нуля, то двумерное нормальное распределение называется невырожденным или собственным. Если ранг квадратичной формы Q t1,t2 равен нулю или

единице, то двумерное нормальное распределение называется вырожденным или несобственным.

Если случайный вектор 1, 2

имеет невырожденное нор-

мальное распределение, то его плотность вероятности равна

 

 

 

 

p x1, x2

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

1 2

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

x1 m1

2

 

 

x1 m1 x2 m2

 

x2 m2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

exp

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

,

 

 

 

 

12

 

 

 

1 2

22

 

 

 

2 1 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где detC c11c22 c122 12 22

1 2 . Смысл параметров плотности

вероятности следующий: i 1,2 : m

, 2

D 2

;

коэффици-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

i

 

i

 

 

 

i

 

 

 

 

 

 

 

 

ент корреляции компонент 1

и 2 случайного вектора определяет-

 

 

 

1 m1 2 m2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ся как

 

 

 

 

 

 

 

 

 

 

 

, справедливо

 

2

 

2

m m .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 2

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Данную плотность вероятности удобно записывать в виде

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

exp

 

Q

 

 

x1

m1, x2 m2

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p x1, x2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 1 2

1 2

 

2

 

detC

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

2

 

 

1

 

 

 

 

 

1

 

 

 

 

2

 

 

 

exp

 

c11

x1 m1

 

2c12 x1 m1 x2

m2

c22 x2

m2

 

 

 

,

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Q 1 t1,t2 cij 1tit j , cij 1

– элементы матрицы C 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

i, j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

420

Необходимо отметить, что плотность вероятности невырожденного нормального распределения сохраняет постоянное значение на эллипсах

1

x1

m1

2

x1 m1 x2

m2

 

x2 m2

2

2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

,

 

 

12

 

1 2

 

 

22

 

2 1 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

называемых эллипсами равных вероятностей, причем вероятность попадания случайного вектора 1, 2 внутрь такого эллипса равна P элл. 1 exp .

Если двумерное нормальное распределение вырождено, то в случае, когда ранг квадратичной формы Q t1,t2 равен нулю, оно сосредоточено в точке m m1,m2 , то есть P m 1. Если

ранг данной квадратичной формы равен единице, то двумерное нормальное распределение сосредоточено на прямой, определяемой собственным вектором матрицы C, который соответствует ее ненулевому собственному значению [343].

Одно из основных направлений исследований в теории вероятностей составляют так называемые предельные теоремы. «Предельные теоремы» – это общее название обширной группы теорем. Важно отметить, что именно предельные теоремы несут в себе большую часть практической значимости теории вероятностей.

Пусть k k 1

– последовательность независимых одинаково рас-

пределенных случайных величин с k :

a, D

2 0. Пусть

 

 

 

 

 

 

 

k

k

 

n

 

 

 

 

 

 

 

Sn n , очевидно, что

Sn na

и

DSn n 2 0. Тогда случай-

 

i 1

n Sn

na

 

 

 

 

 

ная

величина

 

 

является

центрированной

n

n

0 и нормированной D n 1

суммой случайных величин

k kk 1n .

Теорема П4.2 (центральная предельная теорема). Последова-

тельность случайных величин n n 1 , определенных выше, схо-

дится по распределению к стандартной гауссовской случайной ве-

421

личине , то есть P n x x при n равномерно по x R1.

Замечание. Любая последовательность случайных величин, слабо сходящаяся к некоторой гауссовской случайной величине, называется асимптотически нормальной. Центральная предельная теорема устанавливает свойство асимптотической нормальности для последовательности центрированных и нормированных сумм произвольных независимых одинаково распределенных случайных величин, имеющих конечные дисперсии [344].

В частности, одну из категорий предельных теорем для марковских процессов составляют те предельные теоремы, которые являются прямыми обобщениями соответствующих предельных теорем для независимых случайных величин.

Нормальное распределение встречается в очень большом числе приложений. Теоретическое обоснование исключительное роли нормального распределения дают предельные теоремы теории вероятностей. На качественном уровне соответствующий результат может быть объяснен следующим образом: нормальное распределение служит хорошим приближением каждый раз, когда рассматриваемая случайная величина представляет собой сумму большого числа независимых случайный величин, максимальная из которых мала по сравнению с данной суммой.

Приложение 5

Теоремы Хинчина и Леви

Пусть X t является случайной величиной, X 0 0, задано

положительное число t и приведенная случайная величина

t X t t ,

распределенная по приведенному нормальному закону, имеет очень малую вероятность принять очень большое или очень малое значение.

Однако при изменении t происходит своего рода повторение эксперимента, и случайная функция t , как известно, будет поч-

ти наверное иметь бесконечную последовательность нулей. Сход-

422

ным образом могут реализоваться и очень большие значения величины , поэтому нет основания ожидать, что функция t будет

ограничена сверху.

В связи с этим возникает важная проблема определения верхней функции для траектории процесса, то есть такой функции от параметра t, которая при t была бы почти достоверной верхней

гранью траектории t . Данная проблема была решена в 1924 г.

Хинчином [345], и результат, получивший название «закон повторного логарифма», сформулирован в виде следующей теоремы.

Теорема П5.1 (Хинчин). При c 1 почти наверное существует

такое число T , что t T :

 

 

 

 

 

 

t c

 

 

 

.

 

 

 

2t log logt

(П5.1)

При c 1

почти наверное существует такое число T , что t T :

 

t c

 

.

 

 

2t log logt

(П5.2)

Другими словами, справедливо следующее равенство:

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

P lim sup

 

 

 

1

1.

(П5.3)

 

 

 

2t log logt

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Доказательство. Положим tn qn , q 1, n 0,1,2, и

Mn max X t .

t t0

Для x 0 согласно формуле

x

P M x x 2 exp 2 2t d (П5.4)

0

можно получить следующее соотношение:

P Mn x

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

tn

 

 

exp 2

2 d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

exp x2

2

(П5.5)

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

2

 

 

 

2

 

 

 

2

 

 

 

 

 

 

exp

 

2

 

 

 

 

 

 

 

 

 

.

 

 

 

 

x

 

x

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

Следовательно, справедливо следующее неравенство:

423

 

P Mn c

 

 

 

 

 

 

 

exp xn2

2

,

 

 

 

 

 

 

2

 

n

2tn 1 log logtn 1

 

 

 

 

xn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где xn c

2

log n 1 log q ,

то есть выполняется следующее не-

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1 log q c2

q

 

k

 

 

 

 

c

2

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

равенство:

n

 

 

2

 

 

 

 

 

 

 

 

n

 

,

при n ,

 

 

 

 

xn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

logn

 

 

 

 

 

 

 

где k const. Поскольку c 1,

величину q

можно выбрать такой,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

чтобы выполнялось неравенство 1 q c2 .

Тогда ряд

n схо-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 0

дится. Согласно лемме Бореля – Кантелли в данном случае почти

наверное существует такое случайное число N , что

n N вы-

 

Mn c

 

 

,

 

 

 

 

 

 

 

 

полняется

2tn 1 log logtn 1

следовательно,

t T tn

справедливы следующие соотношения:

 

X t Mn

c

 

c

 

 

 

,

 

2tn 1 log logtn 1

2t log logt

(П5.6)

ибо каждое

t T

заключено в некотором интервале

tn 1,tn , где

n N , что и требовалось доказать.

Формула (П5.6)

применима

также и для X t , следовательно, и для

 

X t

 

.

 

 

 

 

Положим

Xn

X tn X tn 1 . Величины X n не зависят друг

от друга, что позволяет использовать лемму Бореля для изучения

 

 

 

 

 

q 1

 

 

 

 

1

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

exp x

 

dx, n

 

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

n

P Xn

xn

 

tn

 

 

 

 

2

 

q

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn

 

 

 

 

 

 

 

1

 

 

 

exp xn2

2

при

xn .

Для

следующих значений:

xn

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c'

 

 

 

 

,

 

 

 

 

 

 

 

xn c'

 

2log logtn

2log nlog q

 

где

 

xn c '

 

 

при

 

 

 

 

2logn

n ,

 

 

имеют

 

место

цепочка

следующих

соотношений:

n

 

 

1

 

 

 

n log q c '2

 

 

 

k

 

 

n log q c '2

, где

k const, n .

xn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

logn

 

 

 

 

 

 

 

Таким образом, если

c ' 1,

вероятность

n

является общим чле-

424

ном расходящегося ряда и, согласно лемме Бореля, для бесконеч-

ного числа номеров n справедливо: xn c '

2

q 1

tn log logtn . С

 

 

 

q

другой стороны, полагая в (П5.6) c 2 1,

можно получить, что

для всех достаточно больших значений параметра n

справедливы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X tn 1

 

2

 

2

 

2

tn

log logtn . Отку-

неравенства

 

2tn 1 log logtn 1

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

да для достаточно больших «избранных» номеров n

следует спра-

ведливость следующих соотношений:

X tn Xn

 

X tn 1

 

 

 

 

c '

q 1

 

2

 

 

 

 

2tn log logtn .

Пусть c 1. Выберем значение

 

 

 

 

q

q

 

 

 

 

 

величины c ' так,

чтобы c c ' 1, а затем возьмем столь большое

 

 

 

 

 

2

 

 

 

 

q, чтобы было c '

 

q 1

 

 

c. Таким образом, для бесконечно

 

 

 

 

 

 

 

q

 

q

 

 

 

 

большого числа значений

n

и для t T tn

справедливо нера-

венство

 

X t c

 

 

 

 

 

 

 

2t log logt

,

(П5.7)

что и требовалось доказать. Формула (П5.7) применима также и

для X t , следовательно, и для

X t

.

Теорема доказана. ▲

Теорема П5.2 (Леви). Если c' 3 2,

то почти наверное сущест-

вует число T такое, что t T :

 

 

 

X t

 

 

 

 

 

 

(П5.8)

2t log logt c 'log log logt .

Если c' 1 2, то T

почти наверное t T :

 

X t 2t log logt c'log log logt .

Данная теорема доказывается аналогично предыдущей. Разница состоит в том, что в первой части теоремы последовательность чи-

сел tn qn следует заменить последовательностью, растущей медленнее, например, выбрать logtn nlogn . Для доказательства

425

второй части, напротив, необходима последовательность, растущая быстрее, чем в доказательстве теоремы Хинчина, например,

tn n! или logtn n log log n .

Приложение 6

Уравнение Ланжевена

Уравнение Ланжевена получено П. Ланжевеном в 1908 г. в теории броуновского движения, его используют для описания случайного воздействия на различные динамические системы, в кинетике фазовых переходов и др. С физической точки зрения справедливо следующее определение [346].

Определение П6.1. Уравнение Ланжевена – это уравнение движения макроскопического тела, взаимодействующего с частицами термостата, влияние которых учитывается при помощи согласованного включения в уравнение силы трения и случайной внешней силы.

Если без учета взаимодействия с термостатом уравнение движения имеет вид:

d 2r

m U r,t 0, dt2

где m – масса частицы, U – потенциальная энергия, то соответствующее уравнение Ланжевена принимает форму

 

 

d 2r

 

dr

 

 

 

 

m

 

 

 

U r,t F t .

(П6.1)

 

dt2

dt

Здесь dr

dt – пропорциональная скорости v dr dt

сила трения,

а F t – случайная сила. Последняя обусловлена одновременным

воздействием на тело большого числа частиц термостата, поэтому с большой точностью ее можно считать нормально распределенной. Среднее значение силы равно нулю, а корреляционная функция

Fi t1 Fj t2

ij t1 t2

зависит лишь от t1 t2. Если время

корреляции k внешней силы, совпадающее по порядку величины со временем одного соударения, k mh, то во всех соотношени-

426

ях, содержащих лишь интеграл от корреляционной функции, ее можно считать пропорциональной -функции: ij 2 ij .

Величина связана с коэффициентом трения , так как и трение и внешняя сила обусловлены взаимодействием тела с термостатом. Эту связь легче всего установить для свободного движения, U 0, тогда при t m имеют место соотношения1:

v2 t 3 m , r2 t 6 t 2 .

Из теоремы о равнораспределении энергии по степеням свободы следует, что v2 t 3kTm, здесь T – абсолютная температура,

откуда kT .

Это соотношение между интенсивностью случайной силы и коэффициентом трения является частным случаем флуктуационно-

диссипативной теоремы. Формула для r2 t соответствует зако-

ну диффузии r2 t 6Kd t, откуда получаются связь Kd 2

между , и коэффициентом диффузии Kd , а также соотношение Эйнштейна Kd kT между коэффициентом трения и коэффициентом диффузии.

Приложение 7

Интегральные преобразования

Определение П7.1. Интегральным преобразованием в общем случае называется линейное взаимнооднозначное отображение U множества распределений (или соответствующих им функций

распределений F или плотностей p , заданных на некотором пространстве Rn ) в множество функций ,U : , при-

чем элемент множества является функцией n комплексных пе-

ременных z z1, zn Cn и имеет следующий вид:

1 Необходимо отметить, что в данном случае, в отличие от гл. 5, рассматривается уравнение Ланжевена в прострастве трех измерений.

427

Соседние файлы в предмете Интегрированные системы управления и проектирования