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

2018

.pdf
Скачиваний:
1
Добавлен:
08.01.2021
Размер:
345.11 Кб
Скачать

21

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 1

 

 

 

 

 

Результаты расчета статистических характеристик

 

 

 

 

 

 

 

Номера интервалов

 

 

 

 

Параметры

 

 

 

 

 

 

 

п/п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

3

4

5

6

7

8

 

9

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

Границы

 

раз-

20-24

24-28

28-32

32-36

36-40

40-44

44-48

48-52

 

52-56

56-58

 

рядов

 

 

 

 

 

 

 

 

 

 

 

 

 

2

Опытные

 

час-

20

31

74

121

153

114

80

27

 

15

6

 

тоты

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

Опытные

 

час-

.031

.048

.115

.189

.239

.178

.125

.042

 

.023

.009

 

тоты

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

Границы

 

раз-

 

 

 

 

 

 

 

 

 

 

 

 

рядов

в

откло-

-20

-16

-12

-8

-4

0

4

8

 

12

16

 

нениях

от

лож-

-16

-12

-8

-4

0

4

8

12

 

16

20

 

ного нуля

 

 

 

 

 

 

 

 

 

 

 

 

 

5

Середина

раз-

-18

-14

-10

-6

-2

2

6

10

 

14

18

 

ряда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

Накопленные

 

 

 

 

 

 

 

 

 

 

 

 

опытные

часто-

0.031

0.079

0.194

0.383

0.622

0.8

0.925

0.967

 

0.991

1

 

сти

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7

Центрир.

 

и

 

 

 

 

 

 

 

 

 

 

 

 

нормиров.

 

От-

 

 

 

 

 

 

 

 

 

 

 

 

клонения

 

для

-2.12

-1.59

-1.05

-0.53

-.003

0.52

1.05

1.58

 

2.11

2.64

 

середины

раз-

 

 

 

 

 

 

 

 

 

 

 

 

рядов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

Таблич.

плот-

 

 

 

 

 

 

 

 

 

 

 

 

ности

вероят-

.042

0.11

0.23

0.35

0.4

0.35

0.23

0.11

 

0.043

0.012

 

ностей

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9

Теоретич.

час-

 

 

 

 

 

 

 

 

 

 

 

 

тоты попадания

14.4

38.3

77

117

135

118

77

38

 

15

4

 

в разряды

 

 

 

 

 

 

 

 

 

 

 

 

 

10

Правые

преде-

-16

-12

-8

-4

0

4

8

12

 

16

20

 

лы разрядов

 

 

 

 

 

 

 

 

 

 

 

 

11

Квадраты

от-

2.18

1.39

0.123

0.11

2.36

0.11

0.08

3.49

 

0.013

0.82

 

клонения

 

 

 

 

 

 

 

 

 

 

 

 

 

12

Правые

преде-

-16

-12

-8

-4

0

4

8

12

 

16

20

 

лы разрядов

 

 

 

 

 

 

 

 

 

 

 

 

13

Центрирован-

 

 

 

 

 

 

 

 

 

 

 

 

ные и нормиро-

 

 

 

 

 

 

 

 

 

 

 

 

ванные

 

откло-

-1.85

-1.32

-0.79

-0.27

0.26

0.79

1.32

1.85

 

2.37

2.9

 

нения

правых

 

 

 

 

 

 

 

 

 

 

 

 

пределов

разря-

 

 

 

 

 

 

 

 

 

 

 

 

дов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14

Интегральная

.033

.094

0.21

0.37

0.59

0.79

0.9

0.97

 

0.991

0.998

 

функция

 

 

 

 

 

 

 

 

 

 

 

 

 

Алгоритм решения

1. Определяют опытные частости (см. лабораторную работу № 1) и за-

22

писывают их в строке 3 табл.1.

2. По полученным значениям строят гистограмму (рис.1). Находят границы разрядов в отклонениях от ложного нуля.

Если a1, b1 – левая и правая границы разрядов соответственно, то:

для первого разряда а1= 20 - 40 = -20.

b1= 24 - 40

= -16;

для второго

а1= 24 - 40 = -16,

b1= 28 - 40

= -12

и т.д.

 

 

 

Результаты записывают в строку 4 табл. 1.

Построение гистограммы

выполняют с использованием середины разрядов, которые определяют как среднее арифметическое между левой и правой границами разряда и записывают в строку 5 табл.1. Тогда для первого интервала

(-20 + (-16))/2 - - 18 (-16 + (-12))/2 = -19 и т.д.

Рис.1. Гистограмма и выравнивающая ее теоретическая кривая распределения опытных частостей

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

23

ния случайной величины в разряды и теоретические частоты (см. лабораторную работу № 1).

Для нашего примера математическое ожидание будет равно М(х) = -2 статистическая дисперсия

D(x) = 52,

а несмещенная оценка для среднеквадратичного отклонения:

б(х) = 7.6.

Плотность вероятности нормального закона выражается с помощью за-

висимости

 

f (x) =1/(σ 2 π ) E ( xM ( x))2 /(2 σ 2 ) .

(3)

4. Вычисляют вероятность попадания случайной величины в разряды. Для этого необходимо вычислить опытные частости (строка 6 табл.1) как сумму опытных частот предыдущих разрядов и определяемого.

Для первого разряда F(x) = Р1 =.031;

Для второго разряда F(x) = Р1 + Р2=.031 +.048 = .079; Для третьего разряда

F(x) = Р1 + Р2 3 =.031 +.048 +.115 =.194.

Далее, применяя метод моментов, необходимо выровнять экспериментальное распределение нормальным законом.

Для построения графика нормального закона производят центрирование и нормирование функции плотности вероятности, для чего предполагают X = 0 и σ = 1. Полученная в результате этого центрирования и нормирования

Функция называется табличной и имеет вид

 

f (t) =1/(σ 2 π Eti 2 / 2 ,

(4)

где ti - центрированные и нормированные отклонения для середины разрядов

ti = (xiсс x) /σ .

(5)

Тогда для первого разряда t1 = (-18 - (-2)) /7.6 = -2.12;

для второго разряда t2 = (-14 - (-2)) / 7.6 = -1.59.

Результаты записываем в строку 7 табл.1.

24

Определив значения центрированных и нормированных отклонений для середины разрядов, вычисляют по формуле (5) табличные плотности вероятностей f(t) и теоретическую вероятность Pi попадания в разряды.

f (t ) = 1 /

2 * 3.14 2.73 ( 2.12 ) 2 / 2 = .042,

 

Pi = f (t)

x /σ = 0.042 * 4 / 7.6 = 0.022 .

 

5. Вычисляют теоретические частоты попадания случайной величины в

разряды по формуле (6)

 

 

mi = Pi N

(6)

mi = 0.022 * 641 = 14.4 (строка 10 табл.1)

6. Для выяснения принадлежности опытных данных к нормальному закону распределения вычисляют значение квадрата Пирсона. Для этого по каждому разряду вычисляют отношение квадратов отклонений теоретических частот от опытных к теоретической частоте попаданий (строка 11 табл.1)

(mi* mi ) / mi = (20 14)2 /14 = 2.18 .

Значения интегральной функции P(t) для нашего примера записаны в строку 14 табл.1 на основании данных табл. П.1. Сумма величин, записанных в строку 11 табл. 1 представляет собой значение квадрата Пирсона в рассматриваемом примере:

n

 

χ2 = (mi* mi )2 / m

(7)

i=1

χ2 = 2.18 +1.39 +... =10.7 .

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

n P * = 1;

i =1

 

 

25

M ( x ) = n

x i * Pi

(8)

i =1

 

 

D ( x ) = n

( x i M ( x )) 2 Pi

 

i =1

 

 

Тогда число степеней свободы будет равно Ч=10 - 3 = 7. Используя стандартную таблицу (табл. П.3), находят уровень значимости

P(χ2 ,Ч) = P(11,7) = 0.088 .

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

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

Pоп (χ2 ,Ч) с заданным значением уровня значимости, т.е.

> q - гипотезанеотвергается Pоп (χ2 ,Ч) q - гипотезаотвергается

или 0.088 > 0.055 < 0.1

Как видно, при заданном уровне значимости q =0.05 гипотеза не отвергается, а при q = 0.1 отвергается, следовательно, мера жесткости проверки правдоподобия принятой гипотезы существенно влияет на принятие окончательного решения.

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

(χ2 Ч). Ч / 2 <3.

 

(9)

Для нашего примера (10.7 7) /

7

=1.98 < 3. Отсюда следует, что со-

 

2

 

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

26

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

интегральной функцией и теоретической.

 

D

= (F * (x) F(x))

max

.

(10)

max

 

 

 

Опытная интегральная Функция строится на основании полученных ранее значений (строка 6 табл.1).

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

ка 12 табл.1) по формуле

 

ti = (bi M (x)) /σi ;

 

ti = (16 + 2) / 7.6 = −1.85 .

(11)

Интегральная функция P(t) (строка 14 табл.1) определяется по табл. П. 2.

На основе числовых значений вероятностных функций, представлен-

ных в строках 6 и 14 табл.1 строят график по формулам

 

n

n

 

F * (x) = Pi

* и P(t) = F(x)dx

(12)

i=1

0

 

Рис. 2. График опытной (1) и теоретической (2) интегральной функции

27

 

Используя полученный график, вычисляют , λ=Dmax

N и с помощью

специальной табл. П. 3 находят вероятность

 

P = P[Dmax N ].

(13)

0.6 - гипотезанеотвергается, При этом если P > 0.6 - гипотезаотвергается

Для нашего примера наибольшее расхождение между опытной и теоретической кривыми интегральной функции находится в сечении t = 38 и составляет

Dmax = 0.622 0.590 0.022 ,

а следовательно, λ = 0.022 641 = 0.56 .

По специальной таблице критических вероятностей для критерия Колмагорова Р(0.56) = 0.90. Поэтому гипотеза о принадлежности опытных данных нормальному закону не отвергается, так как 0.9 > 0.6.

10. Находят величину доверительного интервала разброса среднего результата, отвечающего доверительной вероятности Pg = 95% . Для этого вы-

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

σM*

( x) =σ(x) /

n

(14)

σM*

( x) = 7.57 / 3.3 = 2.39.

 

Половина доверительного интервала составит

 

σ =σm (x) a,

 

(15)

σ = 2.39 1.64 = 3.92 .

 

Значение Pg

= 90% определяется обратной интерполяцией по таблице

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

Jg =38.023.92X <38.02+3.06.

Это значит, что с вероятностью, равной Pg =95%, можно утверждать,

что математическое ожидание составляет не менее 39.51 и не более 45.63.

28

Таким образом, задача решена.

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

ем M (x) = 38.02 .

Порядок выполнения работы

1. Записать условие задачи и исходные данные.

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

2.Вычертить таблицу 1 и заполнить первые три строки, воспользовавшись результатами лабораторной работы № 1.

Провести расчет, используя табличный процессор Excel, и заполнить остальные строки таблицы.

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

3.Проверить принятую гипотезу с помощью критерия согласия.

4.Сделать интервальную оценку математического ожидания и дать заключение.

29

Библиографический список Основная литература

1.Ключников, В. И. Основы научных исследований [Текст] : учеб. пособие / В. И. Ключников ; ВГЛТА. – Воронеж, 2002. – 100 с.

2.Гмурман, В. Е. Теория вероятностей и математическая статистика [Текст] : учеб. пособие / В. Е. Гмурман. – М. : Высш. школа, 1997. – 150 с.

Дополнительная литература 4. Статистические методы для ЭВМ [Текст] / под ред. К. Энслейна. – М. :

Наука, 1986. – 460 с.

30

Приложение Таблица П.1 Вероятность интегральной функции нормального закона для нормиро-

ванной и центрированной случайной величины

t

Pt

t

Pt

t

Pt

 

 

 

 

 

 

-0.00

0.500

-2.70

0.003

1.40

0.919

-0.10

0.460

-2.80

0.002

1.50

0.933

-0.20

0.420

-2.90

0.001

1.60

0.945

-0.30

0.382

-3.00

0.001

1.70

0.955

-0.40

0.344

-3.10

0.001

1.80

0.964

-0.50

0.308

-3.20

0.000

1.90

0.971

-0.60

0.274

-3.30

0.000

2.00

0.977

-0.70

0.242

-3.40

0.000

2.10

0.982

-0.80

0.211

-3.50

0.000

2.20

0.986

-0.90

0.184

-3.60

0.000

2.30

0.989

-1.00

0.584

-3.70

0.000

2.40

0.991

-1.20

0.116

-3.90

0.000

2.60

0.995

-1.30

0.096

0.00

0.500

2.70

0.996

-1.40

0.080

0.10

0.539

2.80

0.997

-1.50

0.066

0.20

0.579

2.90

0.998

-1.60

0.054

0.30

0.617

3.00

0.998

-1.70

0.044

0.40

0.665

3.10

0.999

-1.80

0.035

0.50

0.691

3.20

0.999

-1.90

0.028

0.60

0.725

3.30

0.999

-2.00

0.022

0.70

0.758

3.40

0.999

-2.10

0.017

0.80

0.788

3.50

0.999

-2.30

0.010

1.00

0.841

3.70

0.999

-2.40

0.008

1.10

0.864

3.80

0.999

-2.50

0.006

1.20

0.864

3.90

1.000

 

 

 

 

 

 

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