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

Учебное пособие 2003

.pdf
Скачиваний:
3
Добавлен:
30.04.2022
Размер:
3.84 Mб
Скачать

течение одного года работы равна 0,001 и не зависит от состояния других элементов. Какова вероятность отказа двух и не менее двух электроэлементов за год?

Решение.

 

 

2

 

 

1

 

 

 

 

 

 

 

 

 

1) P

(2)

 

 

e

 

, так как

np

1000 0,001

1

 

 

 

 

 

1000

 

2!

 

 

2e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2) P

(m

2) 1

P(0) P(1) 1

1

e 1

 

1

e 1 1

2

 

0, 264 .

 

 

 

 

1000

 

 

 

 

 

 

 

0!

1!

e

 

 

 

 

 

 

 

 

 

 

7.4. Числовые характеристики пуассоновского распределения

Пусть СВ X - число наступлений события А в n-й серии испытаний. Закон распределения X рассмотрен ранее. Найдем М(Х), D(Х).

Построим производящую функцию этого распределения

 

 

k

 

z

 

 

e zk

e z 1 .

 

 

k

0 k !

 

Отсюда

z

e z 1 ,Ф'(1)=λ;

z

2e z 1 ,Ф"(1)=λ2.

Таким образом, М(Х)=λ; D(X)=λ2+λ-λ2=λ.

Вывод. В распределении Пуассона математическое ожидание и дисперсия равняются параметру λ.

Пример. Завод отправил на базу 5000 изделий. Вероятность того, что в пути изделие испортится, равно р=0,0002. Найти вероятность того, что на базу прибудет а) ровно три негодных изделия; б) не более трех негодных изделий.

Решение. Используем биномиальное распределение

P

3 C3

0, 0002 3

0,9998 4997 .

5000

5000

 

 

Очевидно, что с помощью такой формулы вычисление вероятности затруднительно. Для упрощения вычислений

71

естественно заменить (приближенно) биномиальное распределение распределением Пуассона:

 

 

n

np

5000 0, 0002 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда

а) P5000 3

 

 

 

e

 

0, 06

,

 

 

 

 

 

 

 

 

 

 

 

 

 

3!

 

1

6e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б) P

m 3

P 0 P 1

P 2

P 3

e 1 1 1

1

1

0,95.

 

 

 

 

5000

 

 

 

 

 

 

 

 

 

 

2

6

 

 

 

 

 

 

 

 

 

 

 

 

 

8. НОРМАЛЬНОЕ РАСПРЕДЕЛЕНИЕ. ЧИСЛОВЫЕ ХАРАКТЕРИСТИКИ (ПАРАМЕТРЫ). ФУНКЦИЯ РАСПРЕДЕЛЕНИЯ ВЕРОЯТНОСТЕЙ. ИНТЕГРАЛ ВЕРОЯТНОСТЕЙ. ВЕРОЯТНОСТЬ ПОПАДАНИЯ В ИНТЕРВАЛ. ПРАВИЛО 3ζ. РАСПРЕДЕЛЕНИЕ МАКСВЕЛЛА

8.1. Нормальное распределение

Непрерывная СВ X называется распределенной по нормальному (гауссовскому) закону с параметрами т R и

ζ>0, если плотность распределения вероятностей имеет вид

 

1

 

 

x

m 2

 

 

f x

 

e

2

2

, x

R .

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

Для краткости говорят, что

 

СВ

распределена по

закону N m, .

 

 

 

 

 

 

 

 

Вероятностный смысл параметров m и ζ будет выяснен позже. Нормальное распределение занимает центральное место среди всех других распределений теории вероятностей.

График плотности распределения N(m,ζ) называют кривой Гаусса. Проведем исследование функции f(x) и выясним, как изменяется форма кривой Гаусса с изменением параметров m и ζ.

1. Область определения D(f)=R, f(х)>0 x R .

72

 

 

 

 

 

1

 

 

x

m 2

 

 

2. lim f x

lim

 

e

2

2

0

,

 

 

 

 

 

2

 

 

 

x

 

x

 

 

 

 

 

 

 

т.е. ось ОХ является горизонтальной асимптотой графика функции.

3. Вычислим производную

 

1

 

 

x

m 2

f x

 

e

2

2

 

 

 

 

 

2

 

 

 

 

 

 

 

2 x

m

.

2

2

 

f x

всюду существует и обращается в нуль при x m .

 

При x<m

 

f x >0, при х>m f x <0. Это означает, что

точка

m,

1

 

есть точка максимума.

 

 

 

 

 

 

2

 

 

 

 

4.Для нахождения интервалов выпуклости, вогнутости

иточек перегиба вычислим вторую производную f (x)

 

 

1

 

 

x

m 2

f " x

 

 

e

2

2

3

 

 

 

 

 

2

 

 

 

 

 

 

 

1

2 x

m 2

 

 

2

2

 

 

 

 

1

 

 

 

 

2

x

m

2

 

x

m 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

е

2

2

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f " x 0

2 x m 2

0,

 

 

x m 2

 

2 ;

 

 

x m

f " x 0 при x

 

 

 

; m

m

;

 

, f(x) вогнута.

f " x

0 при x

 

 

 

m

; m

 

, f(x) выпукла.

Точки

m

;

 

 

1

 

 

 

и

m

;

 

1

 

 

 

являются

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

e

 

 

 

 

 

2

e

точками перегиба.

Итак, график функции f(х) для нормального закона распределения имеет вид

73

Очевидно, что при изменении параметра т кривая Гаусса, не меняя формы, смещается вдоль оси ОХ. При изменении параметра ζ кривая или растягивается в положительном направлении оси ОY (при убывании ζ), или сжимается к оси ОХ (при возрастании ζ).

8.2. Числовые характеристики нормального распределения

В силу симметричности функции плотности вероятностей сразу получаем M(Х)=m. Вычислим дисперсию. По определению

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

x

m 2

D X

 

x m

2

 

 

f x dx

 

x m

2

 

 

 

e

2

2

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x m

2

 

1

 

 

 

 

 

 

x m 2

 

x

m

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

e

2 2 dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2e 2 d .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычислим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2e 2 d

 

 

 

e 2 d

2 / 2

 

 

 

d e 2

 

 

 

 

 

e 2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e 2 dx 2 .

Таким образом, D(Х)=ζ2, то есть ζ является средним квадратическим отклонением нормально распределенной СВ

X.

Вывод. Плотность распределения СВ, подчиняющейся

74

нормальному закону, определяется двумя параметрами т и ζ, являющимися соответственно математическим ожиданием и средним квадратическим отклонением СВ.

8.3. Функция распределения вероятностей

По определению функция распределения вероятностей СB равна

 

 

 

 

 

 

x

 

 

 

x

1

 

 

e-

y-m 2

 

 

 

y

 

m

 

F x P X x

 

 

f y dy

 

 

 

 

 

2σ 2 dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

dy

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

m

 

 

 

1

 

 

 

2

 

 

1

0

 

2

 

 

 

1

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

2 d

 

 

e

2 d

 

 

 

 

 

 

 

e

2 d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

2

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Первый интеграл равен 0,5, а второй находится по таблицам.

 

 

 

1

z

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интеграл вида

z

 

 

e 2 d

называется

 

 

 

 

 

 

2

0

 

 

 

 

 

 

 

интегралом Лапласа, интегралом вероятностей, или стандартной функцией ошибок. Закон распределения в этом случае обозначают N(0,1).

Следовательно, функция распределения вероятностей нормально распределенной СВ выражается через Ф(z) следующим образом:

F x

1

 

x m

.

2

 

 

 

 

 

В этом случае непрерывная СВ распределена по нормальному закону N(m,ζ).

8.4. Простейшие свойства интеграла вероятностей

Свойство 1. Ф(0)=0.

Свойство 2. Ф(-z)=-Ф(z), то есть функция Ф(z) нечетная.

75

Свойство 3. Ф(+∞)=0,5.

График функции имеет вид

Ф(z)

0,5

 

0

 

z

 

-0,5

Свойство 4. Функция Ф(z) очень энергично стремится к своему пределу 0,5 при z→∞. Это стремление определяется следующей асимптотической формулой:

 

 

1

 

 

z2

1

1

1 3

1 3 5

 

 

z 0, 5

 

 

e 2

...

,

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

z

 

z3

 

z5

z7

 

справедливой при достаточно больших z.

8.5. Вероятность попадания в интервал. Правило 3σ

По определению

 

 

 

 

P x1 X x2 F x2 F x1

x2 m

 

x1

m

.

 

 

 

 

Вычислим вероятность отклонения случайной величины X от математического ожидания меньше чем на δ:

P

X-M X

P M X

X M X

 

 

M X

M X

M X

M X

 

 

 

 

 

 

 

2 .

Рассмотрим частный случай, когда δ=3ζ. Тогда

P

X M X

3

2 3 /

2 3 2 0, 49865 0,9973

 

 

 

 

 

Это означает, что с вероятностью 0,997, близкой к

76

единице (практически достоверно), вcе значения нормально распределенной СВ располагаются на интервале длиной 6ζ, симметричном относительно математического ожидания.

f (x)

m 3 0 m

m 3 x

Пример. Размер диаметра втулок подчиняется нормальному распределению с параметром M(Х) = 2,5 см, ζ=10-2 см. В каких границах можно практически достоверно гарантировать размеры диаметра втулок?

Решение. Р(|Х-М(Х)|<δ)=2Ф(δ/ζ)=2Ф(δ/0,01)=0,997. По таблице находим δ/0,01=2,98 и δ=0,0298. α=2,5±0,028 (см), 2,4702≤α≤2,5298.

Замечание. Часто рассматривают другую функцию:

 

 

 

 

 

 

 

 

 

 

z

 

1

1

erf

 

 

 

z

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где erf z

 

 

 

 

e t

dt

2

 

z

2

1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для erf z справедливо erf(-z)=-erf z.

 

 

 

 

Иногда табулируется функция

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

z

 

t2

 

 

 

 

 

 

1

 

 

 

1

 

 

z

 

 

 

0

z

 

 

e

2 dt

 

 

 

z

 

 

 

erf (

 

 

).

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Очевидно, что

 

 

z

1

 

 

z .

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8.6. Распределение Максвелла

Скорость u

молекул идеального газа,

находящегося в

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

77

Максвелла c плотностью вероятностей, определенной формулой

0,

u 0

 

 

fV u

 

 

 

 

1

u , u 0,

2/

3/ 2u 2e 2

где параметр распределения

 

m

 

, т – масса молекул, Т

 

kT

 

 

 

 

 

абсолютная температура, k – постоянная Больцмана.

9. СОВМЕСТНОЕ РАСПРЕДЕЛЕНИЕ НЕСКОЛЬКИХ СЛУЧАЙНЫХ ВЕЛИЧИН

9.1. Многомерные случайные величины

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

Пусть на вероятностном пространстве (Ω,

,Р) заданы

СВ X11(ω), X22(ω), …, Xпп(ω), ω

Ω. Вектор

Х(ω)={X1(ω), Х2(ω),..., Хп(ω)} называется случайным вектором или п-мерной случайной величиной. Хk(ω) (k=1,2,...,п) называются координатами, или компонентами вектора X.

Функция F(x1,x2,...,xn)=Р(X1<x1,

X2<x2,..., Xn<xn),

где

xk R (k=1,2,...,п), (x1,x2,...,xn)

Rn, называется п-мерной

ФР

случайного вектора Х=(Х1,

Х2,…,Хп)

или совместной

ФР

случайной величины 1, Х2,…,Хп).

Для наглядности и краткости записи в дальнейшем будем рассматривать в основном двумерную случайную величину. ФР двумерной СВ (двумерного случайного вектора) имеет вид

F(х,у)=Р(X<х, Y<у).

78

Геометрически ФР F(x,y) задает вероятность попадания точки (х,у) в бесконечный прямоугольник X<х, Y<у с вершиной (х,у) (заштрихованная часть на рисунке).

Y

 

y

 

x

X

 

Рассмотрим кратко основные свойства двумерной ФР.

9.2. Свойства двумерной ФР

Свойство 1. ФР F(х,у) является неубывающей функцией по обоим аргументам F(x1,у)≤F(х2,у) при х1<х2 и F(x,y1)≤F(х,у2)

при y1<у2.

Свойство 2. F(х,у) непрерывна слева по каждому из аргументов.

Свойство 3.

F(-∞,у)= lim F (x, y) =0; F(x,-∞)= lim F (x, y) =0.

x y

F

,

lim F x, y 1 .

 

 

x

 

 

y

Свойство 4.

Р{х1≤X<х2,у1≤Y<y2}=F(x2,y2)-F(x1,y2)-F(x2,y1)+F(x1,y1)

(вероятность попасть в прямоугольник).

Свойство 5.

F x F x,

lim F(x, y) ; F y F

, y lim F (x, y).

 

y

x

Это условие согласованности означает, что ФР отдельных компонент двумерного случайного вектора могут быть найдены предельным переходом из совместной ФР этих компонент (маргинальные распределения).

1. Докажем первое из неравенств свойства 1.

79

Рассмотрим событие

А={X<х2; Y<у}, В={X<х1; Y<у}, С={х1≤X<x2; Y<у}.

Y

 

 

y

B

A

 

x1

x2

 

Y

 

 

 

 

y2

B

 

 

 

y1

 

A C

 

X

BC

x1

x2

X

Из рисунков видно, что А=В+С, причем события В и С несовместны, то есть В·С=Ø.Поэтому P(A) P(B+C) P(B)+P(C) .

Но Р(В)=Р(X<х1; Y<у)=F(х1), Р(А)=Р(X<х2; Y<у)=F(х2).

Отсюда получаем

F(х2)=F(х1)+Р(х1≤X<x2; Y<у),

или

F(х2)-F(х1)=Р(х1≤X<x2; Y<у)≥0.

Неубывание F(х, у) по аргументу у доказывается аналогично.

2.Для доказательства непрерывности слева по каждому аргументу достаточно зафиксировать один из них, а для второго провести рассуждения, аналогичные случаю одномерной СВ.

3.Свойство 3 следует из определения ФР и возможности предельного перехода для неубывающей функции. Это свойство хорошо иллюстрируется геометрически.

4.Докажем свойство 4.

Обозначим А={х1≤X<x2; у1≤Y<у2}, В={X<x1; Y<у2},

Р(В)=F(х12), С={X<x2; Y<у1}, Р(С)=F(х21), D={X<x2; Y<у2}, Р(D)=F(х22).

Тогда (рис.2) В·С={X<x1; Y<у1} и Р(В·С)=F(х11). D=А+В+С и P(D)=P(A+B+C)=Р(А)+Р(В+С)=Р(А)+Р(В)+Р(C)-

-Р(ВC).

80