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

Основы теории вероятностей. Что следует знать студенту-математику

.pdf
Скачиваний:
2
Добавлен:
07.09.2026
Размер:
2 Мб
Скачать
☆
2.8. Производящие функции
71
пусты, а на шестую нанесены точки в количестве от 1 до 6, при чем число точек на всех костях различно. Чему равны математиче ское ожидание и стандартное отклонение суммы очков, выпавших на шести таких костях за один бросок? (Ответ: 3,5 и 3,555.)
Задача 2.52. Вернитесь к задаче 2.46. Используя результаты за дачи 2.45, вычислите стандартное отклонение для количества ко робок, которые нужно купить, чтобы собрать полный набор кар точек. (Ответ: 61,951.)
Задача 2.53. Вернитесь к задаче 2.47. Чему равно стандартное отклонение для числа бросков, которые нужно сделать, чтобы по лучить все шесть возможных исходов? (Ответ: 6,244.)
2.8. Производящие функции
1
В примере 2.10 мы использовали подсчет, чтобы найти функцию распределения вероятности для суммы очков на верхних гранях при броске двух костей. Такой подход ведет к громоздким вычис­лениям, когда нужно решить задачу для трех, не говоря уже о пяти или десяти костях. Почти магический подход к решению подоб­ных задач состоит в том, чтобы использовать производящие функ­ции. Производящие функции изобрел швейцарский гений Лео­нард Эйлер (1707—1783) для облегчения вычислений в сложных задачах. Суть их очень проста и является примером того, что часто простыеидеи—самыелучшие.
Производящая функция
неотрицательной случайной ве
()
Gz
X
личиныX, принимающей только целые значения, — это степен ной ряд вида
¥
Gz PX kz
() ( )==
å
X
k
=
0
k
для
||z £ 1
.
-
-
-
-
-
-
-
-
Производящая функция
единственным образом опреде
Gz
()
X
ляет функцию распределения вероятностиX:r-я производная
вточке
()
Gz
X
() ( )1 =
GEX
'
X
.ЕслиXиY— независимые случайные величины
z = 0
равна
rPX r!( )=
(проверьте это). Также
с неотрицательными целыми значениями, то
1
При первом чтении этот раздел можно пропустить.
-
72 Глава 2. Основы теории вероятностей
GzGzGz
XY X Y+
=() () ()
для
то есть производящая функция для суммы
||z £ 1
XY+
,
— это произведе ние производящих функцийXиY. Это очень важный результат. Чтобы его доказать, заметим сначала, что по правилу подстановки из раздела 2.6
¥
поэтому
Xk
Ez zPX k
() ( )==
() ( )=
Gz Ez
X
X
å
k
=
0
для любого
для любого
. Далее, еслиXиY—неза
||z £1
||z £ 1
.
висимые случайные величины, то для любых функцийfиgслу чайные величины доказанное, тогда из сказанного следует, что случайные величины
X
Y
и
z
независимы для любого
z
будет доказано, что
fX()иgY()
EVW EV EW() ()()=
величинVиW. Поэтому для любого
XY XY X Y
GzEz Ezz EzEz GzGz
XY
+
=== =() ( ) ( ) ( ) ( ) () ().
также независимы. Примем это как
. В приложении к этой главе
||z £1
для независимых случайных
,
||z £1
XY+
Теперь применим метод производящих функций для нахожде­ния функции распределения вероятности суммы очков на верхних гранях при одном броске трех костей. Пусть ков на верхней граниi-й игральной кости для щая функция
Gzzzzzzz i
X
задается выражением
X
i
161
() ,,.=+++++ =
i
1
1
23456
6
6
1
6
6
— количество оч-
X
i
i = 12 3,,
1
для
6
.Производя-
123
-
-
-
Поскольку
,
и
X
X
1
2
независимы, производящая функция
X
3
для суммы очков на верхних гранях равна
161
G z zzzzzz
++
XX X
123
æ
=+++++
() для| | .z £ 1
ç è
1
1
23456
6
6
1
6
6
3
1
ö ÷
6
ø
Вычислив произведение многочленов с помощью стандартных компьютерных программ, получим
G z zzz zzz
XX X
++
123
1
1
10 11 12 13 14
+++ + + +
zz z z z z
8
1
=+++++() ++
216
25
8
216
1
345 678
72
7
12
1
365108
5
72
5
72
5
15 16 17 18
108
+++zz z.
36
7
72
1
1
72
25
216
9
z
1
216
2.8. Производящие функции
73
Функцию распределения вероятности суммы можно считать прямо с этого выражения, поскольку в нем коэффициент при это вероятность того, что сумма для
k = 318,,K
.
XX X
++
12 3
примет значение
k
z
Задача 2.54. Есть 8 симметричных шестисторонних игральных костей. Пять из них имеют по шесть очков на двух из шести гра ней, а оставшиеся три кости — по шесть очков на трех из шести граней. Какова функция распределения вероятности для числа шестерок, выпадающих при броске таких 8 костей? (Ответ:
4
æ
,,,,,, , ,
ç
24322243522432381251087764873194413194411
è
944
ö
.)
÷ ø
Задача 2.55. Случайная величинаX, имеющая пуассоновское распределение с параметромl, обладает функцией распределения вероятности
=
Gz e
X
() /!==-ll
PX k e k
z()()
--l 1
— производящая функция дляX.Каковафункция
k
для
k = 01,,K
. Докажите, что
распределения вероятности для суммы двух независимых пуассо­новских случайных величинXиYс параметрамиaиb? (Ответ: это функция распределения вероятности для пуассоновского рас­пределения с параметром
Задача 2.56.Пусть
нулли с
PX
()==
i
1
XX
1
i
+
1
функция распределения вероятности для суммы
16137
æ
вет:
,,,,,)
ç
360516171441481720
è
.)
ab+
,,K
15
и
PX
ö ÷
— независимые величины Бер-
()==-
01
i
.)
ø
1
для
i
+
i = 15,,K
1
XX
15
++L
.Какова
?(От
—
k
-
-
Приложение: доказательства для математического ожидания и дисперсии
В этом приложении представлены доказательства свойств, данных в разделах 2.6 и 2.7 для математического ожидания и дисперсии, а также некоторые дополнительные результаты. Доказательства даются для случая, когдаXиY— дискретные случайные величи ны, которые могут принимать только конечное число значений. ПустьI— множество всех возможных значенийX,аJ—множе ство всех возможных значенийY. На данный момент примите без доказательства следующий основной результат, который дальше
-
-
74
Глава 2. Основы теории вероятностей
будет доказан: для любой функции ние случайной величины
EgXY gxyPX x Y y
(( , )) (,)( )===
Двойную сумму
å
n i
++
L
()aa
iim
1
1
=
. Можно переставлять порядки суммирования,
gXY(,)
åå
yJxI
ÎÎ
n i
gxy(, )
задается выражением
m
a
åå
j
==
11
если количество членов каждой суммы конечно:
n
m
=
åå
j
.
a
ij
i
==
11
математическое ожида
и
следует читать как
ij
n i
åå
Теперь легко можно доказать следующие свойства:
Свойство 1. Свойство 2. Свойство 3.
,еслиXиYнезависятдруготдруга.
aиb
EaX bY aEX bEY()()()+= +
EXY EXEY() ()()=
var( ) var( ) var( )aX bY a X b Y+= +
,еслиXиYнезависятдруготдруга.
22
для константaиb.
для любых констант
Доказательство свойства 1. Применив основной результат для EgXY(( , ))
кфункции
gxy ax by(,)=+
, вы получите, что
EaX bY()+
равно
()( )
ax by P X x Y y
+===
åå
yJxI
ÎÎ
(
===
===+
axP X x Y y
åå
yJxI
ÎÎ
ax PX xY y b
åå
ÎÎ
xI yJ
()
и
и) ( )
и PX x Y y
+===
byP X x Y y
åå
ÎÎ
yJxI
(),==
åå
xIyJ
ÎÎ
и
и
-
m
=
a
ij
j
==
11
где во втором члене последнего выражения порядки суммирова ния перставлены местами. Далее используйте формулу
PX x PX x Y y
() ( ).== = =
Эта формула срезу следует из аксиомы 3 раздела 2.1. Анало
å
yJ
Î
гичное соотношение имеет место для
EaX bY a xPX x b yPY y aEX bEY
() () ()()().+= =+ == +
Доказательство свойства 2. По определению независимых слу чайных переменныхXиY,
åå
xI yJ
ÎÎ
PX x Y y PX xPY y()()()=====и
и
PY y()=
. Это дает
для всех возможных значенийxandy. Далее, применив основной резу льтат для
EgXY(( , ))
при
gxy xy(,)=
, вы получите, что
EXY()
равно
-
-
-
-
2.8. Производящие функции
xyP X x Y y xyP X x P Y y
()()()
=== = =
åååå
yJxI yJxI
ÎÎÎÎ
xP X
=
(
åå
ÎÎ
xI yJ
и
xyPYyEXEY
) ( ) ()().
===
75
Доказательство свойства 3. Для простоты обозначений запи шем деление
как
EX()
var V E V() ( )=-
m
X
,а
EY()
22
как
m
. Применив альтернативное опре
m
Y
для дисперсии случайной величины с математическим ожиданиемm, и воспользовавшись свойством 1, вы получите
varaXbY EaXbY EaXbY
()(())(())
+= + - + =
22
aEX abEX
() (
=+
2 YbEY a b
Дальше используйте независимость
EXY EXEY() ()()=
22 2222 22
aEX ab bEY a ab b
2
aE
[(
=
Важный частный случай свойства 3 состоит в том, что
var aX b a var X() ()+=
по свойству 2, поэтому
22() ()
++---=
mm m mm m
XY X XY Y
222222 2
XbEYavarXbvarY
)][()] () ().-+ -= +mm
XY
2
для любых константaиb.Этоможноуви-
22
22 2
)()( ).+-+
mm
XY
иY.Этодаствам
X
var aX bY()+
равна
деть, если взять в качествеYвырожденную случайную величину с
PY()==11
,длякоторой
var Y()= 0
(доказательство этого частно-
го случая дайте сами).
Остается доказать формулу
EgXY(( , ))
. Прием, который мы для этого применим, состоит
gxyPX x Y y
å
(,)( )
xy
(,)
в том, чтобы ввести случайную величину
(,)( )
gxyPX x Y y
å
xy
(,)
=
å
EZ EgXY
==
==
и
é ê ê
ë
z
() ((,)
gxyPX x Y y
åå
xy g xy zz
(, ): (, )
zP
åå
xy g xy z z
(, ): (, )
можно записать как
(, ) ( )
=
(
==
и
XxYy zPZz
=== ==
и) ()
=
).
ù
=
ú ú
û
==
и
ZgXY= (,)
для
.Тогда
-
-
V
Если положить
YX=
и
gxy gx(, ) ()=
, этот результат служит до казательством формулы подстановки из раздела 2.6 как частного случая.
На бис приведем здесь два знаменитых неравенства, в которых используются понятия математического ожидания и дисперсии. Для неотрицательной случайной величиныYс конечным матема тическим ожиданием неравенство Маркова утверждает, что
-
-
76 Глава 2. Основы теории вероятностей
EY
()
PY a
()
³£
для любой константы
a
a > 0
.
Доказательство этого факта удивительно простое. Для фикси рованной
, и 0 в противном случае. Тогда
Ya³
цательностиYполучается
, пусть индикаторная величинаIравна 1, если
a > 0
EI PY a() ( )=³
,поэтому
YaI³
EY aEI aPY a() () ( ),³=³
. Из-за неотри
что дает нам искомое неравенство. Красота в простоте! Русский математик Андрей Андреевич Марков доказал это неравенство в 1889 году. Оно легло в основу многих других полезных нера венств из теории вероятностей. Для любой случайной величины с конечным математическим ожиданиемmи конечной дисперсией
2
, неравенство Чебышёва утверждает, что
s
2
s
PX c
(| | )-³ £m
для любой константы
2
c
c > 0
.
Оно названо в честь русского математика Пафнутия Львовича Чебышёва (1821—1894), доказавшего это неравенство в 1867 году. Позднее его знаменитый ученик Андрей Марков дал более про­стое доказательство: положив в неравенстве Маркова и
2
,выполучите
ac=
PX c PX c
(| | ) (( ) )
22
EX
[( ) ]
-
ms
2
cc
2
=mm
YX=-()m
2
.-³ = - ³ £
2
Ковариация и коэффициент корреляции
-
-
-
X
2
Как вычислить
var X Y()+
,еслиXиYне являются независи мыми случайными величинами? Чтобы это сделать, вам понадо бится понятие ковариации. Ковариация двух случайных величин иYзадается формулой
cov( , ) [( ( ))( ( ))]XY E X EX Y EY=- -
Обратите внимание на то, что cov X Y(,)
можно записать в эквивалентном виде
cov(,) ( ) ()()XY EXY EXEY=-
cov( , ) var( )XX X=
.
.Формулудля
,
-
-
X
2.8. Производящие функции
77
если записать +
EXEY()()
и использовать свойство линейности оператора матема
( ( ))( ( ))XEXYEY--
ввиде
XY XE Y YE X--() ( )
+
тического ожидания. Из определения ковариации немедленно сле дует, что
cov( , ) cov( , )aX bY ab X Y=
для любых константaиb. Также, используя тот факт, что EXY EXEY() ()()=
для независимыхXиY(см. свойство 2), можно
получить
cov( , )XY = 0
,еслиX и Y — независимые случайные величины.
Однако обратное к этому результату не всегда верно. В качест
ве контрпримера возьмемX, которая принимает значения-1и1 с одинаковыми вероятностями
1
,иположим
YX=
2
.Тогда,по-
2
3
скольку
cov X Y E X E X E X(,) ( ) ()( )=- =
EX()= 0
32
и
EX()
при том, чтоXиYне являются
0
,выполучите
0=
независимыми.
Теперь мы можем написать формулу для дисперсии двух лю-
бых случайных величинXиY:
var( ) var( ) cov( , ) var( )XY X XY Y+= + +2
.
Для доказательства вам потребуется провести некоторое коли
чество алгебраических вычислений. Обозначим для краткости
m
X
EX= ()иm
Y
. Применив определение дисперсии и ли
EY= ()
нейность математического ожидания, вы получите
var( ) [( ( )) ]
XY EXY
+= +- + =
[( ) ] [( )(
EX EX Y
=-+ ­=+ +
XXYY
var( ) cov( , ) var( ),
2
mm
2 -+ - =
XX
2
mm
XY
2
mm
EY
)] [( ) ]
YY
2
-
-
-
-
-
чтоинужнобылодоказать.Формулудля
var X Y()+
можно обоб
щить на случай любого конечного числа случайных величин:
n
n
2++ = +
å
i
1
-
.
ji
ij
n
var( ) var( ) cov( , )XX X XX
L
1
åå
ni
111
i
==+=
В качестве иллюстрации, чему равна дисперсия количества це
лых чисел, которые останутся на своих местах при случайной пе
-
-
-
Глава 2. Основы теории вероятностей
78
рестановке чисел от 1 доn? Эту величину можно записать в виде
å
i
in=
1
,где
X
ло свое исходное положение, и Поскольку
=
PX X
()
var X
()=-
i
. Это приводит к ответу (сюрприз!)
ij¹
EX EX PX
===
11
ij
11
n
, если при перестановке целое числоiсохрани
X
=1
i
в противном случае.
= 0
X
i
2
() () ( )
====
iii
1
и
для всехi,и
2
n
nn
()
-
1
для всех
cov X X
1
1
для всехiи
n
(, )
=
ij
var X
,выполучите
ij¹
1
nn
()
-
1
æ ç
å
i
ç
1
in=
è
EX X
()
1
для всех
-
2
n
ö ÷
(докажи
1
= ÷ ø
ij
те это!).
Еще одно очень важное понятие в статистике — это понятие коэффициента корреляции. Единицы измерения cov падают с единицами измерения величин
EX()иEY()
.Поэтомуча-
(,)XY
не сов-
сто более удобно использовать коэффициент корреляции величин иY. Этот статистический показатель определяется как
cov( , )
XY
r
(,)
XY
=
ss
()()
XY
,
-
=
-
X
если
s()X > 0иs()Y > 0
. Коэффициент корреляции — это безраз-
мерная величина, для которой выполняется свойство
-£ £11r(,)XY
.
Алгебраическое доказательство этого мы здесь опускаем. Слу чайные величиныXandYназываются не коррелирующими,если
r(,)XY = 0
. Независимые случайные величиныXиYвсегда явля ются не коррелирующими, так как их ковариация равна нулю. Од нако обратное не всегда верно. Тем не менее
r(,)XY
часто
используется как мера зависимостиXиY.
Задача 2.57. Прямая линейной регрессии методом наименьших
квадратов для
yEY bxEX=+-() ( ( ))
мизирует выражение
минимизирующее значениеb? (Ответ:
по отношению к
Y
задается соотношением
X
, где коэффициент наклона прямойbмини
2
EY EY bX EX[( ( ) ( ( ))) ]---
bXY
=r
. Чему равно это
Y
s
(,)
()
s
()
X
.
)
-
-
-
-
ГЛАВА 3
ПОЛЕЗНЫЕ РАСПРЕДЕЛЕНИЯ ВЕРОЯТНОСТЕЙ
Эта глава познакомит вас с тремя важными дискретными рас пределениями вероятностей и с тремя важными непрерывными распределениями вероятностей. Дискретные распределения — би номиальное, гипергеометрическое и пуассоновское, а непрерыв ные — нормальное, равномерное и экспоненциальное. Дается представление о практической значимости этих распределений.
3.1. Биномиальное и гипергеометрическое распределения
Случайная величинаX,принимающая значения
01,, ,K n
биномиальное распределение с параметрамиnиp,если
n
æ
ö
где
PX k
() ()==
n
æ
ö
— биномиальный коэффициент. Это распределение возни
ç
÷
k
è
ø
knk
pp k,,,n
ç
÷
k
è
ø
-
-=
101для K
,
кает в вероятностных задачах, которые можно представить в виде последовательности физически независимых испытаний, где каж дое испытание имеет два возможных исхода — «успех» (S)и«не удача» (F). Успешный исход случается с вероятностьюp,анеудач ный — с вероятностью
. Случайная величинаX, определенная
1 - p
как общее количество успехов вnиспытаниях, имеет биномиаль ное распределение. Этот результат легко вывести. Вероятность того, что заданная последовательность изkуспехов и будет иметь место, равна
knk
()1 -
pp
-
;например,при
вероятность того, что получится последовательность
, имеет
неудач
nk­n = 5иk = 3
,равна
SSFSF
-
-
-
-
-
-
-
-
Глава 3. Полезные распределения вероятностей
80
pp p p p p p´´- ´´- = -() () ()111
32
можных последовательностей сkуспехами и
n
æ
ö
но
, так как биномиальный коэффициент
ç
÷
k
è
ø
. Общее количество всех воз
неудачами рав
nk-
n
æ
ö
—эточислоспо
ç
÷
k
è
ø
собов выбратьkразличных элементов изnдоступных элементов, см. раздел 1.1.
Математическое ожидание и стандартное отклонение для
задаются формулами
EX np()=
и
s() ( )Xnpp=-1
.
Проще всего вывести эти формулы, используя индикаторные
случайные величины. ЗапишитеXввиде
XI I
=++
L
1
,где
n
I
равно 1, еслиk-е испытание было успешным, в противном случае
равно 0. Случайная величина с распределением Бернулли
I
k
имеет
()=
EI p
k
и
s() ( )Ipp
=-1
k
, см. задачу 2.44. Все
I
k
I
не зави­сят друг от друга. Применив свойство линейности математическо­го ожидания и
EI EInp I p
-закон для квадратного отклонения, вы получите
n
n
æ ç
åå å
ç
== =
k
11 1
è
n
ö ÷
k
÷ ø
() иs () .1 -´pn
==
k
k
n
æ
ö
ç
÷
=
k
ç
÷
k
è
ø
Биномиальное распределение широко распространено и имеет
множество приложений.
-
-
-
X
k
k
Пример 3.1. Военное устройство раннего обнаружения постро ено в пустыне. Устройство состоит из семи детекторов, включая два резервных. Если работает меньше пяти детекторов, устройство перестает функционировать. Каждые шесть месяцев устройство проверяется, и все детекторы заменяются на новые. С вероятно стью 0,05 любой заданный детектор откажет между проверками. Все детекторы задействованы, и их работа не зависит друг от дру га. Какова вероятность того, что в какой-то момент между провер ками система перестанет функционировать?
Решение. Пусть случайная величинаXобозначает число детек торов, которые прекратят работать между двумя проверками. Тог да случайная величина с параметрами
n = 7иp = 005.
имеет биномиальное распределение
X
. Вероятность того, что система пре
кратит функционировать между проверками, равна
-
-
-
-
-
-
-
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]