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

Математическая статистика. Проверка гипотезы о виде закона распределения. Практикум

.pdf
Скачиваний:
0
Добавлен:
11.08.2026
Размер:
569 Кб
Скачать

 

 

 

 

 

 

 

nl n

 

1

 

b* xl 1 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*

a

*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

n5 n

 

 

 

 

b* x4 50 0,0343 38,984

34 8,5 .

 

 

 

b

*

a

*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для вычисления наблюдаемого значения критерия Пирсона набл2

составим расчётную табл. 2.10.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2.10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Номер

Интер-

 

 

Час-

Гипотетиче-

 

Теорети-

n nP

 

n nP

 

2

n nP

2

интер-

вал

 

 

 

тота

ская вероят-

 

ческая

i

i

 

 

 

 

 

 

 

 

 

 

i

i

 

 

i

i

 

вала i

 

 

 

 

ni

 

ность Pi

 

частота nPi

 

 

 

 

 

 

 

 

nPi

 

1

[10; 16)

 

 

11

 

0,212

 

 

 

10,6

0,4

 

 

 

0,16

 

 

0,0151

 

2

[16; 22)

 

 

10

 

0,206

 

 

 

10,3

– 0,3

 

 

 

0,09

 

 

0,0087

 

3

[22; 28)

 

 

10

 

0,206

 

 

 

10,3

–0,3

 

 

 

0,09

 

 

0,0087

 

4

[28; 34)

 

 

11

 

0,206

 

 

 

10,3

0,7

 

 

 

0,49

 

 

0,0445

 

5

[34; 40]

 

 

8

 

0,170

 

 

 

8,5

–0,5

 

 

 

0,25

 

 

0,0294

 

6

[10; 16)

 

 

11

 

0,212

 

 

 

10,6

0,4

 

 

 

0,16

 

 

0,0151

 

 

 

 

 

50

 

1

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

набл 0,1066

4) сравним эмпирические и теоретические частоты с помощью критерия Пирсона, приняв число степеней свободы k s 1 r 5 1 2 2 ,

так какравномерноераспределениезависитотдвух r 2 параметров. По табл. П3 (прил. 3) найдем критическую точку правосторонней

критической области кр2 0,952 2 5,99 .

Так как набл2 0,1066 кр2 5,99 , то наблюдаемое значение набл2 меньше критического кр2 , следовательно, нет оснований отвергнуть

нулевую гипотезу о равномерном распределении генеральной совокупности.

Рассмотрим конкурирующую гипотезу H1.

Предположим, что полученное эмпирическое распределение случайной величины похоже на показательное, сдвинутое на 10 единиц вправо по оси 0x. Получим несдвинутое показательное распределение, вычитая из всех членов вариационного ряда число 10, т.е. переходя от случайной величины xi к новой величине ui (табл. 2.11):

ui xi 10 .

wi 1.

31

Таблица 2.11

Номер ин-

[ui ; ui+1)

 

ni

n j

 

Относительные частоты

wi ni

тервала i

 

 

 

 

[ui ;ui 1)

 

 

n

 

 

 

 

 

 

 

 

1

[0;

6)

 

 

11

 

0,22

 

2

[6; 12)

 

 

10

 

0,20

 

3

[12;

18)

 

 

10

 

0,20

 

4

[18;

24)

 

 

11

 

0,22

 

5

[24;

30)

 

 

8*

 

0.16

 

___________

 

 

 

 

 

 

 

 

* Интервалы 24; 30

и 30;

объединены.

 

 

Среднее выборочное несдвинутого распределения u x 30 24,4 10 14,4 .

Если случайная величина распределена по показательному закону с параметром , то её математическое ожидание M u 1 .

Так как M u u 14,4 , то M1u 14,1 4 0,0694 .

Проверим гипотезу о показательном распределении с помощью критерия согласия Пирсона (критерий 2 – «хи квадрат»).

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

P U b e e b ,

где и b концы интервала.

P1 P 0 U 6 e 0 e 0,0694 6 1 0,6592 0,3408 ;

P2 P 6 U 12 e 0,0694 6 e 0,0694 12 0,6592 0,4346 0,2246 ;

P3 P 12 U 18 e 0,0694 12 e 0,0694 18 0,4346 0,2865 0,1481 ;

P4 P 18 U 24 e 0,0694 18 e 0,0694 24 0,2865 0,1889 0,0976 ;

P5 P 24 U e 0,0694 24 e 0,1889 0 0,1889 .

Интервал 24; 30 объединён с интервалом 30; .

Для вычисления наблюдаемого значения критерия Пирсона набл2 составим расчётную табл. 2.12.

32

 

 

 

 

 

 

 

 

 

 

 

Таблица 2.12

 

 

 

 

 

n nP

 

 

 

 

 

 

 

 

Номер

Интер-

Часто-

Гипотети-

Теоретиче-

n nP

 

2

n nP

 

2

интер-

вал

та ni

ческая ве-

ская частота

i

i

 

 

 

 

i

i

 

 

i

i

 

 

вала i

 

 

роятность Pi

nPi

– 6,04

 

 

 

 

 

nPi

 

 

1

[0 6)

11

0,3408

17,04

36,482

 

 

2,141

 

 

2

[6;12)

10

0,2246

11,23

– 1,23

1,513

 

 

 

0,151

 

 

3

[12;18)

10

0,1481

7,41

2,59

 

6,708

 

 

 

0,905

 

 

4

[18;24)

11

0,0976

4,88

6,12

 

37,454

 

 

7,675

 

 

5

[24;+∞)

8

0,1889

9,45

–1,45

 

2,103

 

 

 

0,222

 

 

 

50

1

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

набл 11,09

Сравним эмпирические и теоретические частоты с помощью критерия Пирсона.

Критическую величину

2

2

k найдём по таблице критиче-

 

кр

1

 

ских точек распределения .2 (прил. 3). Уровень значимости задан:0,05 . Подсчитаем число степеней свободы

k l 1 r ,

где l – число групп (частичных интервалов) выборки;

r – числопараметровпредполагаемогопоказательногораспределения.

Так как l = 5, r = 1 (подсчитан всего 1 параметр ), то k 5 1 1 3 . По табл. П3 (прил. 3) найдем критическую точку пра-

восторонней критической области кр2 0,952 3 7,8 .

Так как набл2 11,094 кр2 7,8 , то наблюдаемое значение набл2 превышает критическое значение кр2 . Это означает, что результаты

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

Пример 4

Имеютсярезультаты50 измерений некоторойслучайнойвеличины: 2,2; 2,9; 2,5; 2,7; 3,2; 2,8; 3,1; 3,5; 2,0; 3,6; 2,4; 2,6; 2,8; 2,9; 3,0; 2,9; 3,1; 2,7; 3,9; 2,5; 3,8; 3,3; 2,3; 4,0; 2,4; 2,8; 3,0; 3,5; 3,1; 2,1; 3,6; 3,0; 2,0; 3,1;

2,6; 2,9; 3,3; 3,0; 2,2; 3,8; 2,7; 3,0; 2,1; 2,3; 2,0; 3,3; 2,9; 3,5; 2,8; 3,0.

Принять l = 5.

1. Составим вариационный ряд:

j

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

xj

2,0

2,1

2,2

2,3

2,4

2,5

2,6

2,7

2,8

2,9

3,0

3,1

3,2

3,3

3,5

3,6

3,8

3,9

4,0

nj

1

1

2

2

3

2

2

3

4

5

6

5

1

3

3

2

2

2

1

33

 

 

 

 

nj 50 n .

 

 

 

Составим интервальный вариационный ряд.

 

 

Длина частичного интервала

h

xN x1

 

4 2 0, 4

(табл. 2.13).

 

 

 

 

 

 

 

 

 

 

 

5

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2.13

 

 

Частичный

 

 

 

 

 

Номер

 

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

Сумма частот вари-

 

Плотность

интервала

 

интервал

 

го интервала

ант интервала

 

частоты ni

i

 

[xi ; xi+1)

 

 

x x

 

 

 

 

 

 

x*

 

 

 

 

h

 

 

 

 

i

i 1

ni

n j

 

 

 

 

 

i

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[xi ;xi 1)

 

 

1

 

[2,0; 2,4)

 

 

2,2

 

 

 

 

 

9

 

2,25

2

 

[2,4; 2,8)

 

 

2,6

 

 

 

 

 

9

 

2,25

3

 

[2,8; 3,2)

 

 

3,0

 

 

 

 

 

19

 

4,75

4

 

[3,2; 3,6)

 

 

3,4

 

 

 

 

 

7

 

1,75

5

 

[3,6; 4,0)

 

 

3,8

 

 

 

 

 

6*

 

1,50

___________

 

 

 

 

 

 

 

 

 

 

 

 

 

* Интервалы 3,6; 4,0

и 4,0;

объединены.

 

 

 

2. Построим эмпирическую функцию распределения.

Для того чтобы построить эмпирическую функцию распределения

F xi nni ,

можно считать, что все значения непрерывной случайной величины, попавшие в интервал от xi до xi+1, совпадают со значением xi.

F x0 P X 2,0 0 ;

F x

P X 2, 4

9

0,18 ;

 

 

 

 

 

1

50

 

 

 

 

 

 

 

 

 

F x

P X 2,8 9 9 18 0,36 ;

 

2

50

50

 

 

 

 

 

F x

P X 3,2 18 19

37 0,74 ;

 

3

50

50

 

 

 

 

 

F x

P X 3,6 37 7

44 0,88 ;

 

4

50

50

 

 

 

 

 

F x5

P X 4,0 n nx 4 50 1

0,98

;

 

 

 

n

50

 

 

F x P X 4,0 5050 1.

34

Таким образом:

 

 

 

 

0

при

x 2,0;

 

0,18

при

2,0 x 2,4;

 

 

 

при

2,4 x 2,8;

0,36

F x 0,74

при

2,8 x 3,2;

0,88

при

3,2 x 3,6;

 

 

при

3,6 x 4,0;

0,98

 

 

при

x 4,0.

1,00

Построим график эмпирической функции распределения (рис. 2.7) и гистограмму (рис. 2.8).

F(x)

1,00

0,98

0,88

0,74

0,36

0,18

0

2,0

2,4

2,8

3,2

3,6

4,0

x

Рис. 2.7. График эмпирической функции распределения к примеру 4

3. Точечные оценки для математического ожидания и дисперсии. Найдём точечные оценки для интервального вариационного ряда:

 

 

 

x

l

ni xi*

 

 

 

 

i 1

 

,

 

 

 

 

n

 

 

 

 

 

 

где

x*

xi xi 1

– серединные значения вариант в интервалах.

 

 

i

2

 

 

 

 

 

 

 

 

 

 

35

ni h

 

 

4,75

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,25

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,75

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,50

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4,0

x

 

 

 

 

0

 

 

 

2,0

2,4

2,8

3,2

3,6

 

 

 

 

 

 

 

 

 

 

Рис. 2.8. Гистограмма к примеру 4

 

 

 

x

1

 

 

9

2, 2

9 2,6 19 3,0 7

3, 4 6

3,8)

146,8

2,936 .

50

50

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D*

ni xi2

x 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

D*

1

 

9 2, 22

9 2,62

19 3,02

7 3, 42

6 3,82 2,9362

 

50

 

 

 

 

 

 

 

 

 

 

442,96 8,6201 0, 2391 .

 

 

 

 

 

 

 

 

 

 

 

50

 

 

 

 

 

 

 

 

 

 

 

Найдём исправленную выборочную дисперсию, которая служит несмещённой оценкой генеральной дисперсии:

S2

n

 

D*

50

0, 2391 0,244 ;

S

S2

0,244 0,494 .

n 1

49

 

 

 

 

 

 

4. Выберем закон распределения и по критерию Пирсона проверим гипотезу о его согласованности с результатами наблюдений при уровне значимости 0,05 .

По виду гистограммы – распределение нормальное. Проверим гипотезу о нормальном распределении с математическим ожиданием

M X x = 2,936 исреднимквадратическимотклонением = S = 0,494.

36

Вычислим гипотетические вероятности Pi попадания случайной величины в интервалы ai ; bi по формуле

P a x b b Sx a Sx ,

где Ф x функция Лапласа.

 

P P X 2,4 2,4 2,936

 

 

 

 

 

1

 

 

 

 

 

0, 494

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,09 0,3621 0,5 0,1379;

 

 

 

 

 

 

P P 2,4

X

2,8

2,8 2,936

 

 

2, 4 2,936

 

 

 

 

2

 

 

 

 

 

0, 494

 

 

 

 

0,494

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,28 1,09 0,1103 0,3621 0,2519;

 

 

 

 

P P 2,8

X 3, 2

3, 2 2,936

 

 

2,8 2,936

 

 

 

 

3

 

 

 

 

 

0, 494

 

 

 

 

0, 494

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,53 0,28 0,2019 0,1103 0,3122;

 

 

 

 

 

P P 3, 2

X

3,6

3,6 2,936

 

 

3, 2 2,936

 

 

 

 

4

 

 

 

 

 

0, 494

 

 

 

 

0,494

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,34 0,53 0,4099 0,2019 0,2080;

 

 

 

 

 

 

P P 3,6 X 3,6 2,936

 

 

 

 

 

 

5

 

 

 

 

 

 

 

0, 494

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,5 1,34 0,5 0,4099 0,0901.

 

 

 

 

 

 

 

Для определения набл2

составим табл. 2.14.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2.14

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Номер

Интервал

Час-

Гипотети-

Теорети-

n nP

n nP

 

2

n nP

 

2

интер-

 

 

тота

ческая ве-

ческая

 

 

i

i

 

 

 

 

 

 

 

 

i

i

 

 

 

i i

 

 

вала i

 

 

ni

роятность Pi

частота nPi

 

 

 

 

 

 

 

 

 

nPi

 

 

1

[2,0; 2,4)

 

9

0,1379

 

6,895

 

2,105

4,431

 

 

 

0,64

 

 

2

[2,4; 2,8)

 

9

0,2519

12,595

–3,595

12,92

 

 

 

1,03

 

 

3

[2,8; 3,2)

 

19

0,3122

15,610

 

3,390

11,49

 

 

 

0,74

 

 

4

[3,2; 3,6)

 

7

0,2080

10,400

–3,400

11,56

 

 

 

1,11

 

 

5

[3,6; )

 

6

0,0901

 

4,505

 

1,495

2,235

 

 

 

0,50

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

50

 

1

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

набл 4,02

37

Сравним эмпирические и теоретические частоты с помощью критерия Пирсона, приняв число степеней свободы k l 1 r 5 1 2 2 , так какнормальноераспределениезависитотдвух(r = 2) параметровM(X) и .

По табл. П3 (прил. 3) найдем критическую точку правосторонней критической области кр2 0,952 2 5,99 .

Так как набл2 4,02 кр2 5,99 , то наблюдаемое значение набл2 меньше критического значения кр2 . Это означает, что результаты

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

Рассмотрим конкурирующую гипотезу H1.

Предположим, что полученное эмпирическое распределение случайной величины похоже на показательное, сдвинутое на 2 единицы вправо по оси 0x.

Получим несдвинутое показательное распределение, вычитая из всех членов вариационного ряда число 2, т.е. переходя от случайной величины xi к новой величине ui (табл. 2.15).

 

 

ui

xi

2,0 .

 

 

 

 

 

 

 

Таблица 2.15

 

 

 

 

 

 

Номер ин-

[ui; ui 1)

ni

n j

Относительные частоты

wi ni

тервала i

 

 

[ui ;ui 1)

 

n

 

 

 

 

 

1

[0,0; 0,4)

 

9

 

0,18

 

2

[0,4; 0,8)

 

9

 

0,18

 

3

[0,8; 1,2)

 

19

 

0,38

 

4

[1,2; 1,6)

 

7

 

0,14

 

5

[1,6; 2,0]

 

6

 

0.12

 

wi 1 .

Среднее выборочное несдвинутого распределения u x 2,0 2,936 2,0 0,936 .

Если случайная величина U распределена по показательному за-

кону с параметром , то её математическое ожидание

M U

1

.

 

 

 

 

 

 

 

 

 

 

Так как M U

u

0,936 , то

1

 

1

1,068 .

 

 

 

M U

0,936

 

 

 

 

 

 

 

 

 

 

 

38

Проверим гипотезу о показательном распределении с помощью критерия согласия Пирсона (критерий 2 – «хи квадрат»):

 

l

n

nP

2

набл2

 

i

 

i

 

,

 

nP

 

 

 

i 1

 

i

 

 

 

где ni эмпирические частоты;

ni nPi – теоретические частоты.

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

P U b e e b ,

где и b концы интервала.

P1 P 0,0 U 0,4 e 0 e 1,068 0,4 1 0,6523 0,3477 ;

P2 P 0,4 U 0,8 e 1,068 0,4 e 1,068 0,8 0,6523 0,4255 0,2268 ;

P3 P 0,8 U 1,2 e 1,068 0,8 e 1,068 1,2 0,4255 0,2776 0,1479 ; P4 P 1,2 U 1,6 e 1,068 1,2 e 1,068 1,6 0,2776 0,1811 0,0965 ;

P5 P 1,6 U e 1,068 1,6 e 0,1811 0 0,1811.

Интервал 1,6 U 2,0 объединён с интервалом 2,0 U .

Для вычисления наблюдаемого значения критерия Пирсона набл2 составим расчётную табл. 2.16.

 

 

 

 

 

 

 

 

 

 

 

Таблица 2.16

 

 

 

 

 

n nP

 

 

 

 

 

 

 

 

Номер

Интервал

Час-

Гипотети-

Теоретиче-

n nP

 

2

n nP

 

2

интер-

 

тота

ческая ве-

ская часто-

i

i

 

 

 

 

 

i

i

 

 

i

i

 

 

вала i

 

ni

роятность Pi

та nPi

–8,39

 

 

 

 

 

 

nPi

 

 

1

[0,0 ; 0,4)

9

0,3477

17,39

 

70,39

 

 

4,048

 

 

2

[0,4 ; 0,8)

9

0,2268

11,34

–2,34

 

5,476

 

 

0,483

 

 

3

[0,8 ; 1,2)

19

0,1479

7,40

11,60

 

134,56

 

 

18,183

 

4

[1,2 ; 1,6)

7

0,0965

4,83

2,17

 

4,708

 

 

0,975

 

 

5

[1,6; )

6

0,1811

9,06

–3,06

 

9,364

 

 

1,034

 

 

 

 

 

 

 

2,105

 

 

 

 

 

 

 

 

6

[2,0; 2,4)

9

0,1379

6,895

 

4,431

 

 

 

0,64

 

 

 

50

1

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

набл 24,723

Сравним эмпирические и теоретические частоты с помощью критерия Пирсона.

39

Критическую величину 2 2 k найдём по таблице критиче-

кр 1

ских точек распределения .2 (прил. 3). Уровень значимости задан:

0,05 .

Подсчитаем число степеней свободы k l 1 r , где l число групп (частичных интервалов) выборки; r число параметров предполагаемого показательного распределения. Так как l = 5, r = 1 (подсчитан всего 1 параметр ), то k 5 1 1 3 .

По табл. П3 (прил. 3) найдем критическую точку правосторонней критической области кр2 0,952 3 7,8 .

Так как набл2 24,723 кр2 7,8 , то наблюдаемое значение набл2 превышает критическое значение кр2 . Это означает, что результаты

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

40

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]