Материал: Панков Пособие по ТВиМС часть 2

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

P 1

P 2

P 3

...

P 0

 

 

 

 

 

 

 

 

 

 

 

 

P

2

 

 

 

 

 

 

 

...

P 3

P

 

3

...

 

...

 

 

...

 

 

 

 

...

P 1

P 2

 

 

 

Следовательно, E P n . Но при этом

n 0

P n P 0 P 1 ... P m

n 0

P m 1 ... P 2m ... P 3m ...

mP 0 mP m mP 2m ...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mP km m

 

P

 

 

 

 

C

 

k

 

 

 

 

 

 

 

1

 

 

 

 

k 0

 

 

k 1

 

 

 

 

 

 

 

 

mP 1 C .

1 P 1 C

Здесь m , 1

и C взяты из доказательства предыдущей теоремы.

 

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

 

 

Рассмотрим

теперь случайную величину i - значение траектории в

 

i 1

момент принятия решения.

Теорема (тождество Вальда). Во введенных выше обозначениях

E i E E 1 .

i 1

Доказательство. Введем случайную величину i Ind A , где A - событие,

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

 

 

 

i 1 -го шага.

 

 

 

 

i

i i .

 

 

i 1

i 1

 

 

Случайная величина i зависит от случайных величин 1,..., i 1 и не зависит от случайных величин i, i 1,....

E i

P i

1 P i ,

 

 

 

 

 

 

 

 

 

 

 

 

 

E i i

 

 

 

 

E i

 

E i

i

 

 

 

 

 

 

i 1

i 1

 

 

i 1

 

 

 

E i E i E 1 E i

i 1

i 1

 

 

E 1 P i E 1 P n E 1 E .

i 1

n 0

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

Выбор границ А и В

46

Пусть и

- вероятности ошибок первого и второго рода. Пусть A и B

фиксированы.

 

 

 

x1,...,xn :Tn x1,...,xn B

 

 

 

 

1). Обозначим через Wn

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

H1 . Если x1,...,xn Wn , то мы принимаем в качестве истинной гипотезу H1 .

По определению вероятности ошибки первого рода

 

 

 

 

 

 

 

 

 

P H1 |H0 p0 x1 dx1

p0 x1 p0 x2 dx1dx2

....

 

Так как Tn x1,...,xn B , то

W1

 

 

 

 

 

 

 

 

W2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p0 x1 ....p0 xn

1

p1 x1 ....p1 xn .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p1 x1 dx1

 

 

 

p1 x1 p1 x2 dx1dx2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

...

 

 

 

 

 

 

B W

 

 

 

W

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

P H1 | H1

 

1

1 P H0 | H1

 

 

1

.

 

 

 

 

 

 

 

B

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

Мы получаем, что

 

1

. Следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2). Выведем аналогичное неравенство для A.

 

 

 

 

 

 

 

Обозначим

через

 

 

Wn 0

x1,...,xn :Tn x1,...,xn A

область

принятия

гипотезы

H0 .

Если

x1,...,xn Wn 0 ,

то

 

 

мы

 

 

принимаем

в качестве

истинной

гипотезу H0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По определению вероятности ошибки второго рода

 

 

 

 

 

 

 

 

 

P H0 | H1 p1 x1 dx1

p1 x1 p1 x2 dx1dx2

....

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W

 

 

 

 

 

 

 

 

W

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

Так как Tn x1,...,xn A, то

 

 

1

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p1 x1 ....p1 xn A p0 x1 ....p0

xn .

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

0

 

1

1

 

 

p

0

 

 

1

0

2

 

1 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

x

 

 

 

 

x

 

x

dx dx

...

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

W

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A P H0 | H0 A 1 P H1 | H0 A 1 .

Мы получаем, что A 1 . Следовательно,

A 1 .

3). Пусть нам заранее заданы и - вероятности ошибок первого и второго рода. Укажем по ним A и B .

47

Четверка A,B, , полностью задает последовательный критерий отношения правдоподобия.

Возьмем в качестве границ числа a и b такие, что a 1 , b 1 , и

рассмотрим два критерия A,B, , и a,b, ', ' .

b

B

A

a

Утверждение. Во введенных выше обозначениях

' ' .

Доказательство. Согласно доказанным нами ранее неравенствам a 1 ' ' , b 1 ' ' .

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

 

 

 

'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

' ' '

 

 

 

 

 

 

.

1

 

1 ' или

1

 

1 '

' ' '

 

 

 

 

 

 

 

 

 

 

 

 

'

 

 

 

 

 

 

Вычтем из первого неравенства второе и получим

' ' .

Отсюда следует, что ' '.

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

Оценка среднего числа шагов до принятия решения

Пусть случайная величина - число шагов до принятия решения. Согласно доказанным нами ранее теоремам

P 1, E , E i E E 1 .

i 1

E i

Следовательно, E

i 1

,

E

 

 

 

1

 

E 1 |Hk Ek 1 ak , где k 0,1 .

Будем считать, что p0 x p1 x . Тогда i

ln

p1

xi

0

при всех i 1,2,...

p0

xi

 

 

 

 

48

ln A

ln B

Делаем вывод, что при таких условиях в момент окончания процедуры либо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i lnA, либо

i

lnB .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть верна гипотеза H0 . Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

lnA

 

lnB

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда следует, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lnB 1 lnA,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E0 i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E0

lnB

lnA

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E0 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

 

 

 

 

1

 

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При верной гипотезе H1

рассуждения таковы:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

lnA

lnB

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

lnA 1 lnB ln

 

 

 

 

1 ln

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

 

 

 

 

 

1 ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

 

 

 

 

1

 

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Можно отметить также, что если удалось построить последовательный критерий с меньшей, чем у данного критерия, суммой вероятностей ошибок первого и второго рода, то среднее число шагов до принятия решения E0 и

E1 в новом критерии будет больше.

Данный результат, опубликованный в 1948 году, и уточненный в 1976 году Г. Саймонсом, получил название теоремы Вальда – Вольфовица.

49

Тема № 17 Метод статистических испытаний

Случайные величины использовались для решения различных прикладных задач достаточно давно. Примером может служить способ определения числа пи, который был предложен Бюффоном еще в 1777 году, и рассматривался в теме № 4 части I данного курса.

Датой рождения метода Монте-Карло принято считать 1949 г., когда американские ученые Н. Метрополис и С. Улам опубликовали статью «Метод Монте-Карло», в которой систематически его изложили. Название метода связано с названием города Монте-Карло, где в игорных домах (казино) играют в рулетку – одно из простейших устройств для получения случайных чисел, на использовании которых основан этот метод.

Вычислительная техника позволяет легко получать так называемые псевдослучайные числа (при решении задач их применяют вместо случайных чисел); это привело к широкому внедрению метода во многие области науки и техники (статистическая физика, теория массового обслуживания, теория игр и другие). Метод Монте-Карло используют для вычисления интегралов, в особенности многомерных, для решения систем алгебраических уравнений высокого порядка, для исследования различного рода сложных систем (автоматического управления, экологических, биологических и т.д.).

Сущность метода Монте-Карло состоит в следующем: требуется найти значение a некоторой изучаемой величины. Для этого выбирают такую величину , математическое ожидание которой равно a: E a .

Практически же поступают так: производят n испытаний, в результате которых получают выборку объема n из случайной величины , вычисляют

1 n

выборочное среднее x n i 1 xi и принимают x в качестве оценки

(приближенного значения) a* искомого числа a: a a* x.

Поскольку метод Монте-Карло требует проведения большого числа испытаний, его часто называют методом статистических испытаний. Теория этого метода указывает, как наиболее целесообразно выбрать случайную величину , как найти ее возможные значения. В частности, разрабатываются способы уменьшения дисперсии используемых случайных величин, в результате чего уменьшается ошибка, допускаемая при замене искомого математического ожидания a его оценкой a*.

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

Оценка погрешности метода Монте-Карло

50

Источник: https://studfile.net/preview/16534730/