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

Министерство образования и науки Российской Федерации

Федеральное государственное бюджетное образовательное учреждение

высшего профессионального образования «Оренбургский государственный университет»

Кафедра сопротивления материалов и теоретической механики

Н. А. Морозов

СТАТИСТИЧЕСКИЙ АНАЛИЗ РЕЗУЛЬТАТОВ ОДНОФАКТОРНОГО ЭКСПЕРИМЕНТА

Рекомендовано к изданию Редакционно-издательским советом федерального государственного бюджетного образовательного учреждения высшего профессионального образования «Оренбургский государственный университет» в качестве методических указаний для студентов, обучающихся по программам высшего профессионального образования по направлению подготовки 151600.62 Прикладная механика

Оренбург

2014

УДК 001.8 ББК 72.4 (2)

М 80

Рецензент – профессор, доктор технических наук В. М. Кушнаренко

Морозов, Н.А.

М80 Статистический анализ результатов однофакторного эксперимента: методические указания к лабораторной работе по дисциплине «Основы научных исследований» / Н.А. Морозов; Оренбургский гос. ун-т. – Оренбург: ОГУ, 2014 – 18 с.

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

УДК 001.8 ББК 72.4 (2)

©Морозов Н.А., 2014 ©ОГУ, 2014

2

 

Содержание

 

Введение.......................................................................................................................

4

1

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

5

2

Пример выполнения лабораторной работы...........................................................

12

Список использованных источников........................................................................

18

3

Введение

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

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

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

4

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

Цели работы:

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

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

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

следующей последовательности.

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

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

порядковый номер.

 

 

 

2 Вычисление среднего арифметического

х

результатов эксперимента,

принимаемого за результат эксперимента:

 

 

 

 

1

n

 

 

х

 

xi ,

(1)

 

n

 

 

 

i 1

 

 

 

 

 

 

где n - количество результатов эксперимента, xi - i-й результат эксперимента.

Определить х можно с помощью программы Microsoft Excel,

воспользовавшись функцией «СРЗНАЧ».

3 Вычисление оценки среднего квадратичного отклонения результата эксперимента (дисперсии):

 

1

 

n

 

 

 

(xi

x

)2 .

(2)

n 1

 

i 1

 

 

 

 

 

5

Определить можно с помощью программы Microsoft Excel,

воспользовавшись функцией «СТАНДОТКЛОН». 3 Избавление от промахов.

Для определения промахов будем пользоваться двумя критериями:

критерием 3 и критерием Романовского [1].

Критерий 3 основан на том, что при нормальном законе распределения

модуль ни одной случайной погрешности xi x , с вероятностью, практически

равной 1, не может превышать значения 3 . Таким образом, результат измерения считается промахом, если выполняется неравенство:

xi

x

3 .

(3)

Промах исключают из вариационного ряда и заново повторяют вычисления среднего арифметического и дисперсии, затем делают проверку на промахи еще раз. Данную процедуру повторяют до тех пор, пока не останется ни одного промаха.

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

как правило, выбирают из ряда 0.05, 0,02, 0.01, 0.005. Чем больше количество полученных результатов эксперимента, тем меньше должна быть вероятность.

По таблице 1 для заданной вероятности определяют квантиль распределения Стьюдента t. Наличие промаха по критерию Романовского подтверждается, если выполняется следующее неравенство:

xi

x

t .

(4)

Промах исключается из вариационного ряда, проверка повторяется.

4 Определение оценки среднего квадратичного отклонения среднего арифметического:

6

 

 

S

x

 

 

 

.

 

 

(5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

Таблица 1 – Квантили распределения Стьюдента

 

 

 

 

 

 

 

 

 

 

Число

 

 

 

 

Вероятность

 

измерений

0.05

0.02

 

 

 

 

0.01

0.001

10

2,372

2,959

 

 

 

 

3,409

5,014

11

3,327

2,887

 

 

 

 

3,310

4,791

12

2,291

2,829

 

 

 

 

3,233

4,618

13

2,261

2,782

 

 

 

 

3,170

4,481

14

2,236

2,743

 

 

 

 

3,118

4,369

15

2,215

2,710

 

 

 

 

3,075

4,276

16

2,197

2,683

 

 

 

 

3,038

4,198

17

2,181

2,658

 

 

 

 

3,006

4,131

18

2,168

2,637

 

 

 

 

2,997

4,074

19

2,156

2,618

 

 

 

 

2,953

4,024

20

2,145

2,602

 

 

 

 

2,932

3,979

21

2,135

2,587

 

 

 

 

2,912

3,941

22

2,127

2,575

 

 

 

 

2,895

3,905

23

2,119

2,562

 

 

 

 

2,880

3,874

24

2,112

2,552

 

 

 

 

2,865

3,845

25

2,105

2,541

 

 

 

 

2,852

3,819

26

2,099

2,532

 

 

 

 

2,840

3,796

27

2,094

2,524

 

 

 

 

2,830

3,775

28

2,088

2,517

 

 

 

 

2,820

3,755

29

2,083

2,509

 

 

 

 

2,810

3,737

30

2,079

2,503

 

 

 

 

2,802

3,719

40

2,048

2,456

 

 

 

 

2,742

3,602

60

2,018

2,411

 

 

 

 

2,683

3,402

120

1,988

2,368

 

 

 

 

2,628

3,388

 

1,960

2,326

 

 

 

 

2,576

3,291

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

нормальному распределению.

Разобьем вариационный ряд на интервалы. Определим число интервалов

разбиения, используя формулу Старджесса [2]:

m 1 3,31 lgn.

(6)

7

Полученное значение m рекомендуется округлить до целого нечетного числа.

Определим ширину интервала по формуле:

l

xmax xmin

,

(7)

 

 

m

 

где xmax ,xmin - максимальное и минимальное значение вариационного ряда.

Определим границы каждого интервала. Верхняя граница k -го интервала может быть определена из формулы:

xmaxk xmin k l.

(8)

Найдем середины интервалов:

x

 

xmink xmaxk

.

(9)

 

k0

2

 

 

 

 

 

где xmink - нижняя граница k -го интервала (верхняя граница (k 1)-го интервала).

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

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

 

v

 

P

k

.

(10)

 

k

n

 

Результаты вычислений представим в виде таблицы 2.

8

Таблица 2 – Промежуточные значения интервального ряда

Номер

Границы

Среднее

Частота

Статистическая

интервала

интервала

значение

попадания в

вероятность

 

xmin k xmaxk

интервала

интервал

Pk

 

 

xk0

vk

 

...

 

 

 

 

 

По значениям таблицы построим гистограмму. Если гистограмма имеет колоколообразный вид, это свидетельствует о возможности принадлежности распределения экспериментальных данных нормальному закону [3, 4, 5]. Для подтверждения данного результата воспользуемся проверкой нормальности распределения по составному критерию d [2].

Зададимся уровнем значимости по критериям I и II и определим уровень значимости по составному критерию:

 

 

 

I II .

(11)

Вычислим смещенную оценку среднего квадратичного отклонения

результата эксперимента:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

n

 

 

*

 

(xi

x

)2 .

 

 

(12)

n

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

:

 

Проверяем выполнение критерия I. Вычислим значение d

 

 

 

 

1

 

 

m

 

 

 

 

 

 

 

xk0

 

 

 

vk .

(13)

d

 

 

x

 

n

*

 

 

 

 

 

 

 

 

 

 

 

 

k 1

 

9

Гипотеза о принадлежности распределения экспериментальных данных нормальному закону справедлива в случае, если:

 

 

 

 

 

 

 

 

d1 0,5 I

 

d d0,5 I

,

 

 

 

 

 

 

 

 

 

(14)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где d1 0,5 I

, d0,5 I - квантили распределения d , определяемые по таблице 3.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 3 – Квантили распределения d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

d

0,01

 

 

d0,05

 

 

 

 

 

d0,95

 

 

 

 

 

 

d0,99

11

 

 

 

0,9359

 

0,9073

 

 

 

 

0,7153

 

 

 

 

 

 

0,6675

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16

 

 

 

0,9137

 

0,8884

 

 

 

 

0,7236

 

 

 

 

 

 

0,6829

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

21

 

 

 

0,9001

 

0,8768

 

 

 

 

0,7304

 

 

 

 

 

 

0,6950

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

26

 

 

 

0,8901

 

0,8685

 

 

 

 

0,7360

 

 

 

 

 

 

0,7040

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

31

 

 

 

0,8827

 

0,8625

 

 

 

 

0,7404

 

 

 

 

 

 

0,7110

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

36

 

 

 

0,8769

 

0,8578

 

 

 

 

0,7440

 

 

 

 

 

 

0,7167

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

41

 

 

 

0,8722

 

0,8540

 

 

 

 

0,7470

 

 

 

 

 

 

0,7216

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

46

 

 

 

0,8682

 

0,8508

 

 

 

 

0,7496

 

 

 

 

 

 

0,7256

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

51

 

 

 

0,8648

 

0,8481

 

 

 

 

0,7518

 

 

 

 

 

 

0,7291

 

 

 

 

 

 

 

 

 

 

Значения квантилей для

n, не представленных в таблице 3, находятся с

помощью линейной интерполяции.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Проверяем выполнение критерия II. Гипотеза о нормальности

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

 

если не более

 

 

разностей

 

xk0

x

 

 

превзошли

 

 

 

 

 

значение z0,5P , где z0,5P

 

- верхняя

квантиль

интегральной

 

функции

нормированного распределения Лапласа для вероятности 0,5P. Значения и P

определяются в

зависимости от II

 

и n

по

таблице

4.

Затем по

таблице 5

определяется z0,5P .

6 Определение доверительного интервала и запись результата эксперимента.

10

Соседние файлы в папке новая папка 1