Лабораторная работа №3

 

Изучение нормального распределения

случайной величины на доске Гальтона

 

ЦЕЛЬ РАБОТЫ

 

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

 

ТЕОРЕТИЧЕСКИЕ ОСНОВЫ РАБОТЫ

 

Данная работа посвящена определению по ограниченному коли-честву опытов неизвестных параметров, от которых зависит рас-пределение случайной величины. Надо отметить, что любое значе-ние искомого параметра, вычисленное на основании ограниченного числа опытов, всегда будет содержать элемент случайности. Такое приближенное, случайное значение называется оценкой параметра. Например, оценкой для математического ожидания может служить среднее арифметическое наблюденных значений случайной вели-чины в  независимых опытах. При очень большом числе опытов среднее арифметическое будет с большой вероятностью весьма близким к математическому ожиданию. Если же количество опытов  невелико, то замена математического ожидания средним значе-нием приведет к некоторой ошибке. Эта ошибка будет тем больше, чем меньше число опытов. Так же будет обстоять дело и с оцен-ками других параметров. Любая из таких оценок случайна, при пользовании ею неизбежны ошибки. Желательно выбрать такую оценку, чтобы ошибки были по возможности минимальными.

Возьмем случайную величину , закон распределения которой содержит неизвестный параметр . Требуется найти подходящую оценку для параметра  по результатам  независимых опытов, в которых величина  принимает определенное значение. Обозна-чим наблюденные значения случайной величины:

 

.                                (3.3.1)

 

Обозначим как  оценку для параметра . Любая оценка
должна представлять собой функцию величин :

 

.                         (3.3.2)

 

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

 

.                                  (3.3.3)

 

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

.                              (3.3.4)

 

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

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

Рассмотрим задачу о доверительном интервале для математи-ческого ожидания. Пусть произведено  независимых опытов над случайной величиной , характеристики которой (математическое ожидание  и дисперсия ) неизвестны. Для этих параметров по-лучаем оценки:

; .                     (3.3.5)

 

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

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

 

.                           (3.3.6)

 

Применяя распределение Гаусса, получим:

 

,                                    (3.3.7)

 

где       среднее квадратичное отклонение оценки ;

           функция, определяемая из распределения Гаусса.

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

 

.                        (3.3.8)

 

Таким образом, задача построения доверительного интервала решена:

.                          (3.3.9)

 

Равенство (3.3.9) означает, что неизвестное значение математи-ческого ожидания  попадет в интервал  с вероятностью .

При этом необходимо отметить, что ранее мы рассматривали ве-роятность попадания случайной величины в заданный интервал. Здесь дело обстоит иначе: значение математического ожидания величина не случайная, зато случаен интервал . Случайно его положение на оси абсцисс, определяемое его центром ; случай-на и длина интервала, равная , т.к. величина  вычисляет-ся, как правило, по опытным данным. Поэтому в данном случае лучше будет толковать величину  не как вероятность попадания точки  в интервал , а как вероятность того, что случайный ин-тервал  накроет точку .

Мы рассмотрели приближенный метод построения доверительно-го интервала для математического ожидания . Идея точных по-строений доверительных интервалов сводится к следующему. Лю-бой доверительный интервал находится из условия, выражающего вероятность выполнения некоторых неравенств, в которые входит интересующая нас оценка . Закон распределения оценки  зави-сит от самих неизвестных параметров величины . Однако иногда удается перейти от случайной величины  к какой-либо другой функции наблюдаемых значений , закон распределе-ния которой не зависит от неизвестных параметров, а зависит толь-ко от количества опытов  и от вида закона распределения .


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

 

,                              (3.3.10)

где

                      (3.3.11)

 

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

 

ОПИСАНИЕ ЭКСПЕРИМЕНТАЛЬНОЙ УСТАНОВКИ

 

Работа выполняется на так называемой доске Гальтона, изобра-женной на рис.3.3.1.

Установка состоит из вертикальной доски-1, на которой закреп-лены в шахматном порядке стержни-2, служащие для рассеивания шариков, поступающих из хранилища-3, расположенного вверху доски. Под стержнями расположены 15 одинаковых ячеек, разде-ленных перегородками одинаковой высоты-4. Шарики удержива-ются в хранилище стерженьком-5, закрывающим отверстие, через которое шарики высыпаются из хранилища. Лицевая часть доски закрыта стеклом. Выпускное отверстие хранилища шариков распо-ложено над 8-й ячейкой.


 

Если бы не было стержней, то шарики, выпущенные из хранили-ща, попали бы в 8-ю ячейку. В нашем же опыте шарик, соударяясь с рядом стержней, может попасть практически в любую ячейку. Иначе говоря, попадание шарика в ту или другую ячейку носит случайный характер. Если выпустить три шарика, то, скорее всего, они попадут в разные ячейки, номера которых будут отличаться от номера ячей-ки, над которой расположено выпускное отверстие. Чаще всего даже средние значения будут отличаться от истинного значения, т.е. . При повторении этого опыта несколько раз вероятнее всего полу-чатся иные результаты, нежели в первый раз. В каждой серии изме-рений, состоящих из трех опытов , можно найти доверитель-ный интервал по формуле Стьюдента.

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

                                (3.3.12)

 

Если выпустить все шарики,то на доске Гальтона они располо-жатся так, как показано на рис.3.3.2.

 


Итак, мы получаем экспериментальную кривую вероятности  в зависимости от номера ячеек . По форме она напоминает форму распределения шариков по ячейкам (рис.3.3.2). Если опыт прове-ден аккуратно, то полученная кривая должна совпасть с теорети-ческой кривой распределения случайной величины Гаусса.

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

ПОРЯДОК ВЫПОЛНЕНИЯ РАБОТЫ

 

Перед началом работы необходимо убедиться, что в ячейках нет шариков.

Часть 1

 

1. Постепенно выпустить из хранилища все шарики.

2. С помощью линейки измерить высоты получившихся столби-ков шариков в каждой ячейке  и занести результаты, округлив их до 0,5 см, в таблицу 1.

3. Перевернув доску Гальтона, пересыпать все шарики в храни-лище, расположенное в верхней части доски.

 

Таблица 1

п/п

 

 

 

 

 

 

1

2

3

15

 

 

 

 

 

 

 

 

 

Часть 2

 

1. Выпустить из хранилища поочередно 3 шарика и записать в таблицу 3 значения координаты , т.е. номера ячеек, куда попали шарики.

2. Аналогичные опыты проделать еще 4 раза и заполнить четыре таблицы, аналогичные таблице 2 (3,4,5,6).

 

Таблица 2

 

 

 

 

Доверительный интервал

1

2

3

 

 

 

 

 

Среднее

 

 

 

Обработка реЗУЛЬТАТОВ ИЗМЕРЕНИЙ

 

Часть 1

 

1. Подсчитать сумму высот  во всех ячейках .

2. По формуле (3.3.12) определить экспериментальное значение вероятности попадания шарика в каждую ячейку . Резуль-таты занести в таблицу 1.

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

4. Вычислить величину  и занести в таблицу 1.

5. Вычислить и занести в таблицу 2 значения

6. По формуле (1.23) определить среднее квадратичное откло-нение .

7. Для каждого  подсчитать по формуле (1.26) значение функ-ции , где  и .

8. В одних осях координат построить график вероятности  в зависимости от координаты  и график функции . Сравнить полученные кривые.

 

Часть 2

 

1. Найти оценку математического ожидания  координаты вы-пускного отверстия по формуле (3.3.3) для каждой из пяти серий измерений.

2. Найти погрешность :

,

 

а также  для каждого измерения и занести результаты в таблицы 2-6.

3. По формуле (3.3.7), заменив предварительно величину  на коэффициент Стьюдента , и с использованием формулы (3.3.8) вычислить  отдельно для  и , т.е. использовать формулу:

.

 

4. Записать в каждой таблице, соответственно, для двух надеж-ностей ответы в виде:

.

 

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

 

Контрольные вопросы

 

1. Какой процесс называется случайным?

2. Что такое случайная величина?

3. Какая величина называется частотой события?

4. Что называется законом распределения случайной величины?

5. Что такое функция распределения случайной величины?

6. Запишите основные свойства функции распределения случай-ной величины.

7. Укажите основные числовые параметры, характеризующие закон распределения случайной величины и объясните их смысл.

8. Что называется математическим ожиданием случайной вели-чины?

9. Дайте определение величины дисперсии.

10. Что называется средним квадратичным отклонением?

11. В чем удобство использования среднего квадратичного от-клонения по сравнению с использованием дисперсии?

12. Какова связь между функцией распределения и плотностью распределения?

13. Для каких случайных величин существует плотность распре-деления– дискретных или непрерывных?

14. Запишите основные свойства плотности распределения.

15. Что такое кривая распределения?

16. Запишите выражение для функции плотности распределения непрерывной случайной величины.

17. Какова вероятность принятия случайной величиной конкрет-ного значения при дискретном распределении? При непрерывном распределении?

18. Как влияет дисперсия случайной величины на форму кривой распределения?

19. Укажите аналог кривой распределения для дискретных слу-чайных величин.

20. Укажите оценку основных параметров распределения.

21. Какая оценка называется точечной?

22. Какая оценка называется несмещённой?

23. Что такое доверительная вероятность или надежность изме-рения?

24. Как измеряется доверительный интервал для среднего значе-ния измеряемой величины и что он обозначает?

25. В классической физике имеет место классическое распреде-ление Максвелла. Что это за распределение и чем оно отличается от нормального распределения?

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

27. Воспользовавшись условиями предыдущего вопроса, укажи-те механическую интерпретацию дисперсии.

28. Как составляется статистический ряд?

29. Как производится построение гистограммы?

30. Как по гистограмме получить характеристики распределения случайной величины?

 

ForStu / Практика / Физика / МЕТОДИЧКИ+ЛАБЫ(2 курс МАТИ, 5 факультет)

Copyright © 2004-2017, ForStu

Яндекс.Метрика