Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Планирование и обработка эксперимента / Spirin - Metodi planirovaniya i obrabotki 2004.pdf
Скачиваний:
601
Добавлен:
13.05.2017
Размер:
7.84 Mб
Скачать

2.КРАТКИЕ СВЕДЕНИЯ ИЗ ТЕОРИИ ВЕРОЯТНОСТЕЙ

ИМАТЕМАТИЧЕСКОЙ СТАТИСТИКИ

2.1.Случайные величины и параметры их распределений

Поскольку из-за влияния неконтролируемых факторов отклик – это всегда случайная величина, при обработке результатов эксперимента широко используется аппарат теории вероятностей и математической статистики, поэтому напомним некоторые основные понятия и определения этого раздела математики.

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

Например, реализацию такого события, как пять остановок доменной пе-

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

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

Случайная величина может быть дискретной или непрерывной.

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

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

Если при фиксированном наборе уровней всех контролируемых факторов

провести n измерений отклика X, то в результате будет получен ряд хотя и близких, но отличающихся друг от друга значений:

xi (i =1, 2, ..., n),

(2.1)

где xi – i -е измерение величины X;

x1, x2,..., xn – реализация случайной величины X.

17

2. КРАТКИЕ СВЕДЕНИЯ ИЗ ТЕОРИИ ВЕРОЯТНОСТЕЙ …

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

Таблица 2.1

Число остановок доменной печи по месяцам (общее число наблюдений n = 18)

Месяц

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Число

3

4

3

5

5

5

6

4

6

5

5

2

4

6

7

5

6

7

остано-

вок

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В данном примере число остановок доменной печи в течение месяца - это дискретная случайная величина. В первом из n = 18 измерений этой вели-

чины было получено значение x1 = 3, во втором – x2 = 4 и т.д., до x18 = 7. Приведенные в табл. 2.1 значения – это реализация такой случайной величины, как число остановок доменной печи в течение месяца.

Каждому значению дискретной случайной величины X (любому из событий А, когда случайная величина X принимает какое-либо строго определенное значение x), можно поставить в соответствие следующее отношение:

W = m

,

(2.2)

n

 

 

где m – число наблюдений, в которых дискретная случайная величина X оказалась равна x; n – общее количество наблюдений. Величину W называют частотой реализации события А.

В примере 2.1, в шести наблюдениях: i = 4, 5, 6, 10, 11 и 16, количество остановок доменной печи в течение месяца X оказалось равным пяти (X = 5),

следовательно, частота реализации такого события, как пять остановок, равна 6/18 = 0,33. Частоты реализаций для других событий (две, три, четыре и т.д.

остановки) приведены в табл. 2.2.

18

2. КРАТКИЕ СВЕДЕНИЯ ИЗ ТЕОРИИ ВЕРОЯТНОСТЕЙ …

Таблица 2.2

Частота остановок доменной печи

Число остановок x

2

3

4

5

6

7

 

 

 

 

 

 

 

Количество наблюдений m, в кото-

1

2

3

6

4

2

рых реализовалось событие X = x

 

 

 

 

 

 

Частота реализации, W = m / n

0,06

0,11

0,17

0,33

0,22

0,11

 

 

 

 

 

 

 

Если продолжить наблюдения за работой доменной печи в течение еще полутора лет, то, конечно же, совершенно не обязательно, что на протяжении следующих восемнадцати месяцев пять остановок будет снова зарегистрировано ровно в 6 случаях из 18 наблюдений, а частота реализации этого события опять окажется равной 0,33. Однако при возрастании числа повторений одного и того же комплекса условий частота реализации такого события, как, напри-

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

Предел, к которому стремится отношение m/n при неограниченном возрастании числа опытов n, называется вероятностью случайного события.

Вероятность P(А) события А – число от нуля до единицы, которое представляет собой предел частоты реализации события А при неограниченном числе повторений одного и того же комплекса условий.

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

определения (поскольку вероятность того, что она примет какое-либо конкрет-

ное свое значение, стремится к нулю).

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

чайной величины и соответствующими им вероятностями.

19

2. КРАТКИЕ СВЕДЕНИЯ ИЗ ТЕОРИИ ВЕРОЯТНОСТЕЙ …

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

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

Функция распределения F(x)– функция, определяющая для всех действительных х вероятность того, что случайная величина Х принимает значение не больше, чем х.

F(x)= P(X x).

(2.3)

Функция распределения F(x) имеет следующие свойства (рис.2.1, а):

1.Ее ордината, соответствующая произвольной точке х1, представляет собой вероятность того, что случайная величина X будет меньше, чем

х1, т.е. F(x1) = P(Х ≤ x1).

2.Функция распределения принимает значение, заключенное между нулем и единицей:

0 F (x)1.

(2.4)

3.Функция распределения стремится к нулю при неограниченном уменьшении х и стремится к единице при неограниченном возрастании х, то есть

lim F(x)=0,

lim F(x)=1.

(2.5)

x→−∞

x→+∞

 

4.Функция распределения представляет собой монотонно возрастающую кривую, то есть

F(x2)>F(x1), если х2>х1.

(2.5а)

5.Ее приращение на произвольном отрезке (х1; х2) равно вероятности

того, что случайная величина X попадет в данный интервал:

F(x2)F(x1)=P(Xx2)P(Xx1)=P(x1 <Xx2).

(2.6)

а

20

б

 

 

F(x)

 

F(x)

2. КРАТКИЕ СВЕДЕНИЯ ИЗ ТЕОРИИ ВЕРОЯТНОСТЕЙ …

Рассмотрим, какие особенности имеют функции распределения дискретных случайных величин. Пусть Х – дискретная случайная величина, принимающая возможные значения х1, х2,…, хn с вероятностями p1, p2, …, pn . Функция распределения вероятностей этой случайной величины Х равна

F (x)= P( X x) = pk ,

xk

где производится суммирование вероятностей всех возможных значений случайной величины Х, меньших чем х. Такая функция всегда разрывная, ступенчатая (рис.2.1, б): от −∞ до х1 включительно функция равна нулю, в точке х1 происходит скачок на величину p1, и функция остается постоянной до х2 включительно и т.д., то есть возможным значениям случайной величины соответствуют скачки функции, равные вероятностям этих значений. Последний скачок на pn происходит в точке хn, и функция равна единице от хn до +∞. Таким образом, сумма всех скачков равна единице.

Плотность распределения f(x) – первая производная (если она существу-

ет) функции распределения.

f (x )=

dF (x )

.

(2.7)

 

 

dx

 

21

2. КРАТКИЕ СВЕДЕНИЯ ИЗ ТЕОРИИ ВЕРОЯТНОСТЕЙ …

Плотность функции распределения f(x) имеет следующие свойства

(рис.2.2):

f(x)

f(x)

x1 dx

Mx M0 Me

x2

x

Рис.2.2. Дифференциальный закон распределения – плотность распределения f(x)

1.Плотность распределения вероятностей является неотрицательной функцией, т.е.

f (x)0.

(2.8)

Это свойство справедливо, так как F(x)

есть неубывающая функция.

2.Функция распределения случайной величины Х равна определенному интегралу от плотности распределения вероятностей в пределах (−∞,

х):

x

 

F (x )= f ( x ) dx .

(2.9)

3.Вероятность события, состоящая в том, что случайная величина Х

примет значение, заключенное в полуинтервале [x1 ,x2 ], равна опре-

деленному интегралу от плотности распределения вероятностей на

этом полуинтервале:

x 2

 

P (x 1 < X x 2 )= F (x 2 )F (x 1 )= f ( x ) dx .

(2.10)

x 1

 

22

2.КРАТКИЕ СВЕДЕНИЯ ИЗ ТЕОРИИ ВЕРОЯТНОСТЕЙ …

4.Интеграл плотности распределения в бесконечно большом интервале (-∞, + ∞) равен единице:

+∞f (x )dx = P(− ∞ < X ≤ +∞ )= 1,

(2.11)

−∞

 

так как попадание случайной величины в интервал −∞ < Х< + ∞ есть достоверное событие.

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

Например, для случайной величины, которая удовлетворяет так называемому нормальному закону распределения (закону распределения Гаусса), функцию распределения можно записать в виде

F ( x ) =

1

πσ

2

 

x

(x

M

x

)2

 

 

 

 

σ x2

 

 

 

 

e

 

2

 

dx ,

(2.12)

2

 

 

 

 

 

x

− ∞

 

 

 

 

 

 

 

а для случайной величины, имеющей, например, распределение Вейбу-

ла-Гнеденко (используемое для описания результатов экспериментов в случае хрупкого разрушения металла, а также в испытаниях на многоцикловую усталость), функция распределения определяется следующим выражением:

x x

H

b

 

 

 

 

 

 

 

c

 

 

 

 

F(x ) = 1 e

 

 

,

при Х>хн,

 

F(x) = 0,

 

 

 

 

при Х хн.

(2.13)

В функциях (2.12) и (2.13) константы Mx, σx2 и с, b, хн являются параметрами распределений, причем первое из этих двух выражений относится к двухпараметрическому виду закона распределения, а второе, соответственно, – к трехпараметрическому.

Параметр распределения – постоянная, от которой зависит функция распределения.

23

2. КРАТКИЕ СВЕДЕНИЯ ИЗ ТЕОРИИ ВЕРОЯТНОСТЕЙ …

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

Важнейшими параметрами распределения, задающими случайную величину Х, являются ее математическое ожидание Mx (характеризует центр рассеивания) и дисперсия σx2 (характеризует степень рассеивания).

Математическое ожидание Mx – среднее взвешенное по вероятностям значение случайной величины.

Для дискретной случайной величины математическое ожидание определяется выражением

M x =xi pi,

(2.14)

i

 

где хi – значения дискретной случайной величины, а pi = P(X= хi).

Если в условиях примера 2.1 предположить, что pi Wi (см. табл. 2.2), то для математического ожидания такой дискретной случайной величины, как число остановок доменной печи в течение месяца, можно получить следующее значение:

Mx = 2·0,06 + 3·0,11 + 4·0,17 + 5·0,33 + 6·0,22 + 7·0,11 = 4,87.

Для непрерывной случайной величины математическое ожидание определяется интегралом

M x = +∞xf (x )dx ,

(2.15)

−∞

 

где f(x) – плотность распределения непрерывной случайной величины.

Можно отметить, что геометрический смысл математического ожидания

непрерывной случайной величины – это абсцисса центра тяжести фигуры под кривой плотности распределения f(x). Сказанное проиллюстрируем на рис. 2.2,

где видно, что произведение f(x)dx есть площадь элементарного участка под кривой f(x), а x – абсцисса этого участка, т.е. расстояние от начала координат.

24

2. КРАТКИЕ СВЕДЕНИЯ ИЗ ТЕОРИИ ВЕРОЯТНОСТЕЙ …

Следовательно, интеграл (2.15) дает абсциссу центра тяжести всей площади фигуры под кривой f(x).

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

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

Для примера 2.1 (см. табл. 2.2), при условии, что pi Wi, мода Мо числа остановок доменной печи равна 5, поскольку именно этому значению данной дискретной случайной величины соответствует локальный максимум вероятности, равный 0,33.

Медиана Ме – значение случайной величины, для которого функция рас-

пределения принимает значение ½ , или имеет место «скачок» со значения, меньшего чем ½, до значения, большего чем ½.

Таким образом, для дифференциального закона распределения медиана есть такое значение непрерывной случайной величины Х, которое делит пополам площадь под кривой плотности распределения f(x).

В примере 2.1, если предположить, что функция распределения от четырех остановок F(4) (вероятность того, что число остановок доменной печи в течение месяца будет не более четырех) равна 0,06 + 0,11 + 0,17 = 0,34 , а функция распределения F(5) = 0,34 + 0,33 = 0,67, то медианой Ме такой дискретной случайной величины, как число остановок доменной печи в течение месяца, будет значение Ме = 5.

Дисперсия случайной величины σx2 – математическое ожидание случай-

ной величины

- Mx)2.

 

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

щим математическим выражением:

 

n

(xi M x )2 p(xi ).

 

σ x2 =

(2.16)

i=1

Впримере 2.1 (опять же, если предположить, что pi Wi) значение дис-

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

σx2 = (2 – 4,87)2·0,06 + (3 – 4,87)2·0,11 + (4 – 4,87)2·0,17 +

25

2.КРАТКИЕ СВЕДЕНИЯ ИЗ ТЕОРИИ ВЕРОЯТНОСТЕЙ …

+(5 – 4,87)2·0,33 + (6 – 4,87)2·0,22 + (7 – 4,87)2·0,11 = 1,7931.

Для непрерывной случайной величины дисперсия определяется выраже-

нием

σx2 = +∞(x Mx )2 f (x)dx,

(2.17)

−∞

 

где х – значения непрерывной случайной величины Х; f(х) – плотность распределения; Mx – математическое ожидание.

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

Среднее квадратичное отклонение σx – неотрицательный квадратный ко-

рень из дисперсии.

σx = + σx2 .

(2.18)

Для примера 2.1 среднее квадратичное отклонение числа остановок до-

менной печи в течение месяца равно σ x = + 1,7931 =1,34.

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

Квантиль порядка P, хр – значение случайной величины, для которого

функция распределения принимает значение P или имеет место «скачок» со

значения, меньшего чем P, до значения, большего чем P:

 

F(xp) = P.

(2.19)

Из этого определения квантиля следует, что медиана Ме – это квантиль

порядка ½, т.е. Ме = х0,5.

 

Вероятность попадания случайной величины Х в интервал [ хP1, хP2 ] рав-

на

 

P(x P1 < X x P2 ) = P(X x P2 ) P(X x P1 ) = F(x P2 )F(x P1 )= P2 P1.

(2.20)

В примере 2.1 квантиль порядка 0,95 числа остановок доменной печи

скорее всего равен семи х 0,95 = 7, поскольку F(6) ≈ 0,06 + 0,11 + 0,17 + 0,33 + 0,22 = 0,89, а F(7) ≈ 0,89 + 0,11 = 1,00.

26