Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
ВКР-2010(Поляков).pdf
Скачиваний:
971
Добавлен:
15.05.2015
Размер:
4.39 Mб
Скачать

181

*

 

 

 

1

 

 

 

x mx

f (xi

 

f (x)

 

 

 

 

 

 

 

 

exp

 

 

 

 

 

 

2 2

 

2

 

 

 

 

 

 

 

x1

x2

xk

xk+1

x

 

а) распределение, близкое к нормальному (Гаусса)

 

 

f*(xi)

 

 

 

 

 

 

f(x) = exp ( - x)

 

 

x1

x2

 

xk

xk+1

x

 

б) распределение, близкое к экспоненциальному

 

 

 

f*(xi)

f(x) = 1 / ( - )

 

 

 

 

 

 

 

 

 

 

x1

x2

xk

xk+1

x

в) распределение, близкое к равномерному распределению

Рисунок 6.37 - Гистограммы статистических распределений

6.5.2 Расчет параметров распределения случайной величины

Каждой числовой характеристике случайной величины соответствует статистическая аналогия: математическому ожиданию – статистическое среднее; дисперсии - статистическая дисперсия; среднему квадратическому отклонению – эмпирический стандарт и т.д.

Если число опытов сравнительно невелико (n < 100), то статистические характеристики можно определить из простого статистического ряда по формулам:

182

- статистическое среднее

 

 

n

 

 

mx 1

n

xi

,

(6.94)

 

i 1

 

 

 

 

 

 

- статистическая дисперсия

 

 

 

 

n

 

 

 

x 1

 

(xi mx )2 .

(6.95)

D

n

 

 

 

i 1

 

 

 

 

 

 

При большом количестве опытов (n > 100) для определения статистических характеристик целесообразно воспользоваться статистическим рядом. Тогда для числовых характеристик статистического распределения имеем:

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

mx x j p*j

;

(6.96)

 

 

 

 

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

(x j mx )2 p*j ,

(6.97)

 

 

 

 

 

 

 

Dx

 

 

 

 

 

 

 

 

 

 

 

j 1

 

 

где

x

 

1

(x

 

x

j 1

)

- «представитель» j-го разряда; сред-

 

 

j

2

 

j

 

 

 

 

 

 

нее значение j-го разряда.

6.5.3 Выравнивание статистического ряда

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

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

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

- нормальное распределение (распределение Гаусса)

183

 

 

1

 

 

 

(x mx )2

 

;

(6.98)

f (x)

 

 

 

exp

 

 

 

 

 

 

 

 

 

2 x

 

 

2

 

 

 

 

 

 

 

2 x

 

 

 

- экспоненциальное распределение

0

при x 0;

(6.99)

f (x)

 

e x при x 0;

 

 

 

 

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

 

 

1

при x ;

(6.100)

 

 

 

 

f (x)

 

 

 

0

при x

и x ,

 

 

 

 

 

 

 

где σx – среднее квадратическое отклонение; mx – математическое ожидание;

λ – интенсивность потока событий; α, β – границы интервала распределения СВ (для распреде-

ления равномерной плотности).

Пределы изменения случайной величины для распределения равномерной плотности определяются из уравнений:

mx

 

и Dx

 

( )2

.

12

 

 

2

 

 

Решая эти уравнения, получим

 

 

 

 

 

 

mx

 

и mx 3Dx .

(6.101)

3Dx

Нормальный закон распределения (закон Гаусса) используется для описания большинства случайных величин в технике, природе; экспоненциальное распределение – для описания времени наработки до отказа; закон равномерной плотности – для описания случайных величин, вероятность появления которых примерно постоянна на данном интервале.

Следует отметить, что в науке и технике используется более десяти законов распределения [3.10].

Отметим, что гистограмма распределения может быть описана не только известными законами распределения, но и специфическими для данного распределения зависимостями. Подбор аналитических функций для гистограмм плотности может быть произведен методом наименьших квадратов. Но любая аналитическая функция f(x), с помощью которой выравнивается статистическое распределение, должно обладать основными свойствами плотности распределения:

 

184

 

f (x) 0;

 

(6.102)

 

f (x) dx 1.

 

 

Предположим, что вид аналитической функции f(x) вы-

бран. Известно, что в выражение этой функции входят несколько параметров a, b, c и т.д. Требуется подобрать эти параметры так, чтобы функция f(x) наилучшим образом описывала данный статистический материал. Согласно одному из методов подбора – методу моментов – параметры a, b, c выбирают с таким расчетом, чтобы несколько важнейших числовых характеристик (моментов) теоретического распределения были равны соответствующим статистическим характери-

стикам, например: mx

 

mx ; Dx

 

 

x

 

и т.д.

 

D

 

 

Для экспоненциального закона распределения, который

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

 

 

1

 

или

 

1

.

(6.103)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

Dx

 

 

 

 

 

 

 

Для нормального закона и закона равномерной плотно-

сти, которые определяются двумя параметрами, -

 

 

mx mx ; Dx

 

 

x .

 

(6.104)

 

D

 

6.5.4 Проверка правдоподобия гипотезы о выборе закона распределения случайной величины

Критерии согласия отвечают на вопрос о степени согласованности (соответствия) теоретического и статистического распределений. Одним из наиболее распространенных крите-

риев является критерий К. Пирсона, или критерий «хи - квад-

рат», или χ2 - квадрат.

К. Пирсон показал, что мера расхождения

 

k

n p j )

2

 

 

 

 

2

 

(m j

 

(6.105)

 

n p j

 

 

j 1

 

 

зависит только от числа «степеней свободы»

r = k s, (6.106)

где k - число разрядов;

s - число независимых условий («связей»), наложенных на выбор теоретического распределения;

pj - вероятность попадания СВ в j - ый интервал.

185

Для экспоненциального закона распределения число связей s = 2, для нормального закона и закона равномерной плотности s = 3, а именно:

 

k

 

 

1 я связь : p*j 1;

(6.107)

j 1

 

 

 

2 я связь :mx

mx ;

 

 

 

 

3 я связь :

D

х

Dx .

 

Вероятности попадания СВ в j -ый интервал определяются

p j F(x j 1 ) F(x j ) .

(6.108)

Функция распределения определяется зависимостями: - для экспоненциального закона

0

при x 0;

(6.109)

F (x) 1 e x

 

при x 0;

 

 

 

 

- для закона равномерной плотности

 

0

 

при x ;

(6.110)

x

 

F (x)

 

 

при x ;

 

 

 

 

 

 

 

1

при x .

 

 

 

 

 

 

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

вы и равны pj = 1 / k.

Для нормального закона распределения вероятности попадания в j -ый интервал определяются по формуле

 

 

 

1

 

x

j 1

m

x

 

x

j

m

x

 

 

p

 

 

 

Ф

 

 

 

Ф

 

 

 

.

j

 

 

 

 

 

 

 

 

 

 

 

2

 

 

x

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Значения функции Лапласа

 

2

 

 

Ф(x)

 

e t 2 dt ,

 

 

 

 

 

 

 

0

где t x mx .x

(6.111)

(6.112)

Значения функций Лапласа (Приложение Б) и распределения χ2 (Приложение Д) табулированы и приведены в руководствах по теории вероятностей, а также могут быть определены по встроенным функциям MS Excel ВЕРОЯТНОСТЬ(…) и ХИ2ОБР(Вероятность, степени свободы).

Фрагменты таблицы «Значения χ2 в зависимости от r и p» приведены в таблицах 6.12 и 6.13.

186

Используя таблицу «Значения χ2 в зависимости от r и p», для каждого значения χ2 и числа «степеней свободы» r можно найти вероятность того, что величина, распределенная по закону χ2, превзойдет фактическое значение этой величины.

Если вероятность р меньше вероятности нижней границы рн.г. = 0,1, т.е. р < рн.г., то выбранное теоретическое распределение не может быть принято. Если вероятность р больше вероятности верхней границы рв.г. = 0,9, т.е. р > рв.г., то следует исключить «подчистку данных», когда некоторые результаты изменяются или исключаются.

Пользоваться таблицей возможно двумя способами:

1 По числу «степеней свободы» r и значению χ2, полученному по формуле (6.105), определяют точное значение вероятности р. Если выполняется условие рн.г. < р < рв.г. , то выбранное теоретическое распределение может быть принято в качестве рабочей гипотезы. Для r = 8 и χ2 = 5,53 последовательность действий показана на номограмме (таблица 6.12).

Таблица 6.12 - Значения χ2 в зависимости от r и p (Фрагмент)

p

0,99

рв.г. = 0,9

0,7

рн.г. = 0,1

0,001

r

 

 

 

 

 

 

 

 

 

1

0,000

0,016

 

2,71

10,83

2

0,020

0,211

 

4,60

13,82

8

1,65

3,49

2,82

5,53

13,36

21,6

2 Для верхней рв.г. = 0,9 и нижней рн.г. = 0,1 границ вероятностей и числа «степеней свободы» r определяют диапазон значений χ2. Если значение χ2, полученное по формуле (6.105) попадает в диапазон, то выбранное теоретическое распределение может быть принято в качестве рабочей гипотезы. Для r = 8 последовательность действий показана на номограмме (таблица 6.13).

Таблица 6.13 - Значения χ2 в зависимости от r и p (Фрагмент)

p

0,99

рв.г. = 0,9

рн.г. = 0,1

0,001

r

 

 

 

 

 

 

 

 

 

1

0,000

0,016

2,71

10,83

2

0,020

0,211

4,60

13,82

8

1,65

3,49

 

 

13,36

21,6