Архитектура Аудит Военная наука Иностранные языки Медицина Металлургия Метрология
Образование Политология Производство Психология Стандартизация Технологии


СПОСОБ ОБРАБОТКИ ЭМПИРИЧЕСКИХ ДАННЫХ



Общие положения

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

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

Впервые этот метод предложил Гаусс. Рассмотрим применение этого метода для эмпирических формул, описывающих различные формы зависимостей.

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

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

Исходные данные для расчётов приведены в приложении 7.

Линейная функция

Пусть дано n – точек с координатами , представляющих данные эксперимента – таблица 1.

 

Таблица 1 – Экспериментальные данные в общем виде

i
n

 

Очевидно, что эмпирическую формулу нужно искать в виде линейной функции

(1)

Причём коэффициенты a и b нужно подобрать так, чтобы суммарное отклонение принимало минимальное значение.

Задача нахождения минимального значения функции – есть классическая задача дифференциального исчисления. Ординате i -ой точки на прямой соответствует точка - рисунок 1.

Рисунок 1 – График линейной функции

Отклонение находится по выражению

.

Сумма квадратов отклонений

(2)

 

Для нахождения экстремума функции двух переменных необходимо осуществить следующую последовательность операций.

Найти частные производные функции S по переменным a и b для чего продифференцировать S по a и b

,

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

.

или

(3)

.

Выражения (3) представляют собой систему двух линейных уравнений с двумя неизвестными a и b.

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

Прямая y = ax+b при найденных из системы (3) числовых значений и есть прямая наилучшего приближения к точкам по критерию метода наименьших квадратов. Эту прямую так же называют линией регрессии на .

Пример. Определить зависимость износа коренных шеек коленчатого вала двигателя Д-240 от наработки. Исходные данные приведены в таблице 2.

Нормальная система (3) согласно данным таблицы 2 имеет вид

(4)

В результате решения имеем

a=0, 09623; b=0, 02492

Значение параметров и можно определить по формулам:

,

Таблица 2 – Износ коренных шеек коленчатого вала двигателя Д-240 от наработки

Наработка , тыс. мотто-ч Износ , мм
0, 2 0, 05 0, 04 0, 010
0, 8 0, 08 0, 64 0, 064
1, 4 0, 17 1, 96 0, 238
2, 0 0, 21 4, 0 0, 420
2, 6 0, 27 6, 76 0, 702
3, 2 0, 37 10, 24 1, 184
3, 8 0, 40 14, 44 1, 520
4, 4 0, 42 19, 36 1, 848
18, 4 1, 97 57, 44 5, 988

 

Эмпирическое уравнение линейной функции имеет вид:

(5)

Квадратичная функция

Эмпирическую формулу нужно искать в виде квадратичной функции

(6)

Сумма квадратов отклонений имеет вид

Частные производные по , и

,

,

.

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

,

, (7)

.

При составлении нормальной системы данные опыта удобно свести в таблицу 3.

 

Таблица 3 – Экспериментальные данные в общем виде

 

Решая нормальную систему, найдём значения a, b, c для эмпирической формулы.

Параметры уравнения параболы второго порядка можно определить по следующим формула:

, , ,

где ; ; ;

,

число членов эмпирической зависимости; значение зависимой переменной y, корреляционно связанной с независимой переменной X; отклонение членов ряда независимой переменой X от их средней x, отнесённые к значению промежутка между членами ряда x, то есть . При нечётном числе членов ряда отклонения x выражается числами 0, 1, 2, 3…, а при чётном числе членов – 0, 54 1; 1, 5; 2, 5; 3, 5…

Рассмотренный приём можно использовать для аппроксимации опытных данных многочленами более высоких степеней. Коэффициенты многочлена в этом случае так же определяют из системы уравнений, получаемых приравниванием к нулю частных производных от суммы квадратов отклонений по соответствующим параметрам. Этот метод можно применять и в том случае, когда аппроксимация неполиноминальная.

Гиперболическая функция

Уравнение гиперболической функции

(8)

Сумма квадратов отклонений имеет вид

Приравнивая их к нулю, получим

(9)

При вычислении коэффициентов данные удобно свести в таблицу 4.

Таблица 4 – Экспериментальные данные в общем виде

k
n

Числовые значения параметров a и b найдём путем решения системы уравнений (9).

Параметры a и b уравнения (9) можно определить по формулам:

; ;

где

Показательная функция

Данные опыта могут быть аппроксимированы показательной кривой

(10)

Для получения параметров и прологарифмируем обе части функции. При этом учтем, что логарифм произведения равен сумме логарифмов сомножителей, а логарифм степени равен произведению показателя степени на логарифм ее основания. Затем следует найти величины и .

Сумма квадратов отклонений

(11)

Частные производные по и

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

(12)

Для решения системы удобно составить таблицу 5.

Таблица 5 – Экспериментальные данные в общем виде

i
k
n

Найденные из нормальной системы (12) значения и с помощью таблицы антилогарифмов позволяют легко определить a и b.

Показательная функция

Экспериментальные данные могут быть аппроксимированы показательной кривой в виде

( 13)

Прологарифмируем это выражение по основанию e

(14)

Сумма квадратов отклонений

.

Частные производные по основаниям b и ln a

(15)

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

Для решения системы уравнений удобно свести данные в таблицу 6.

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

 

Таблица 6 – Экспериментальные данные в общем виде

i
...
n

СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ

 

1. Тимирязев К.А. Наука и демократия. М.: Сельхозгиз, 1963. – 346с.

2. Павлов И.П. Избранные труды. М.: Учждгиз, 1964. – с.37.

3. Капустин В.П. Основы научных исследований и патентоведения: лекции к курсу. Тамбов. – 1996. – 84с.

4. Бронштейн И.Н., Семендяев К.А. Справочник по математике для инженеров и учащихся вузов. – М.: Наука 1981. – 719 с.

5. Доспехов Б.А. Методика полевого опыта. – М.: Агропромиздат, 198. – 31 с.

6. Кравченко В.С., Трубилин Е.И., Курасов В.С., Куцеев В.В. Основы научных исследований: Учебное пособие / КГАУ, Краснодар. – 2002. - 126 с.

7. Королюк В.С., Портенко Н.И., скороход А.В., Турбин А.Ф. Справочник по теории вероятностей и математической статистике. М.: Наука, 1985. – 640 с.

8. Рябушкин Т.В., Ефимова М.Р., Ипатова И.М., Яковлева Н.И. Общая теория статистики. - М.: Финансы и статистика, 1981.-279 с.

 

 

ПРИЛОЖЕНИЕ

 

Таблица 1 – Значения t - критерия (Стьюдента)

Число степеней свободы   Уровень значимости, %
0, 5
12, 71 63, 66 -
4, 30 9, 93 31, 60
3, 18 5, 84 12, 94
2, 78 4, 60 8, 61
2, 57 4, 03 6, 86
2, 45 3, 71 5, 96
2, 37 3, 50 5, 41
2, 31 3, 36 5, 04
2, 26 3, 25 4, 78
2, 23 3, 17 4, 59
2, 20 3, 11 4, 44
2, 18 3, 06 4, 32
2, 16 3, 01 4, 22
2, 15 2, 98 4, 14
2, 13 2, 95 4, 07
2, 12 2, 92 4, 02
2, 11 2, 90 3, 97
2, 10 2, 88 3, 92
2, 09 2, 86 3, 88
2, 09 2, 85 3, 85
2, 08 2, 83 3, 82
2, 07 2, 82 3, 79
2, 07 2, 81 3, 77
2, 06 1, 80 3, 75
2, 06 2, 79 3, 73
2, 06 2, 78 3, 71
2, 05 2, 77 3, 69
2, 05 2, 76 3, 67
2, 05 2, 76 3, 66
2, 04 2, 75 3, 65
50. 2, 01 2, 68 3, 50
100. 1, 98 2, 63 3, 39
1, 96 2, 58 3, 29

Таблица 2. Значения F – критерия (Фишера) на 5 %-ном уровне значимости

Степени свободы для меньшей дисперсии (знаменателя) Степени свободы для большей дисперсии (числителя)
1.
2. 18, 51 19.00 19, 16 19, 25 19, 30 19, 33 19, 36 19, 37 19, 38 19, 39 19, 41 19, 45 19, 47 19, 49
3. 10, 13 9, 55 9, 28 9, 12 9, 01 8, 94 8, 88 8, 84 8, 81 8, 78 8, 74 8, 54 8, 58 8, 56
4. 7, 71 6, 94 6, 59 6, 39 6, 26 6, 16 6, 09 6, 04 6, 00 5, 96 5, 91 5, 77 5, 70 5, 66
5. 6, 61 5, 79 5, 41 5, 19 5, 05 4, 95 4, 88 4, 82 4, 78 4, 74 4, 68 4, 53 4, 44 4, 40
6. 5, 99 5, 14 4, 36 4, 53 4, 39 4, 27 4, 21 4, 15 4, 10 4, 06 4, 00 3, 84 3, 75 3, 71
7. 5, 59 4, 74 4, 35 4, 12 3, 97 3, 87 3, 79 3, 73 3, 68 3, 63 3, 57 3, 41 3, 32 3, 28
8. 5, 32 4, 46 4, 07 3, 84 3, 69 3, 58 3, 50 3, 44 3, 39 3, 34 3, 28 3, 12 3, 03 2, 98
9. 5, 12 4, 26 3, 86 3, 63 3, 48 3, 37 3, 29 3, 23 3, 18 3, 13 3, 07 2, 90 2, 80 2, 76
10. 4, 96 4, 10 3, 71 3, 48 3, 33 3, 22 3, 14 3, 07 3, 02 2, 97 2, 91 2, 74 2, 64 2, 59
11. 4, 84 3, 98 3, 59 3, 36 3, 20 3, 09 3, 01 2, 95 2, 90 2, 86 2, 79 2, 61 2, 50 2, 45
12. 4, 75 3, 88 3, 49 3, 26 3, 11 3, 02 2, 92 2, 85 2, 80 2, 76 2, 69 2, 50 2, 40 2, 35
13. 4, 64 3, 80 3, 41 3, 18 3, 02 2, 92 2, 84 2, 77 2, 72 2, 67 2, 60 2, 42 2, 32 2, 26
14. 4, 60 3, 74 3, 34 3, 11 2, 96 2, 85 2, 77 2, 70 2, 65 2, 60 2, 53 2, 35 2, 24 2, 19
15. 4, 54 3, 60 3, 29 3, 06 2, 90 2, 79 2, 70 2, 64 2, 59 2, 55 2, 48 2, 29 2, 18 2, 12
16. 4, 49 3, 63 3, 24 3, 01 2, 85 2, 74 2, 66 2, 59 2, 54 2, 49 2, 42 2, 24 2, 13 2, 07
17. 4, 45 3, 59 3, 20 2, 96 2, 81 2, 70 2, 62 2, 55 2, 50 2, 45 2, 38 2, 19 2, 08 2, 02
18. 4, 41 3, 55 3, 16 2/93 2, 77 2, 66 2, 58 2, 51 2, 46 2, 41 2, 34 2, 15 2, 04 1, 98
19. 4, 38 3, 52 3, 13 2, 90 2, 74 2, 63 2, 55 2, 48 2, 43 2, 38 2, 31 2, 11 2, 00 1, 94
20. 4, 35 3, 49 3, 10 2, 87 2, 71 2, 60 2, 52 2, 45 2, 40 2, 35 2, 28 2, 08 1, 96 1, 90
21. 4, 32 3, 47 3, 07 2, 84 2, 68 2, 57 2, 49 2, 42 2, 37 2, 32 2, 25 2, 05 1, 93 1, 87
22. 4, 30 3, 44 3, 05 2, 82 2, 66 2, 55 2, 47 2, 40 2, 35 2, 30 2, 23 2, 03 1, 91 1, 84
23. 4, 28 3, 42 3, 03 2, 80 2, 64 2, 53 2, 45 2, 38 2, 32 2, 28 2, 20 2, 00 1, 88 1, 82
24. 4, 26 3, 40 3, 01 2, 78 2, 62 2, 51 2, 43 2, 36 2, 30 2, 26 2, 18 1, 98 1, 86 1, 80
25. 4, 24 3, 38 2, 99 2, 76 2, 60 2, 49 2, 41 2, 34 2, 25 2, 24 2, 16 1, 96 1, 84 1, 77
26. 4, 22 3, 37 2, 98 2, 74 2, 59 2, 47 2, 39 2, 32 2, 27 2, 22 2, 15 1, 95 1, 82 1, 76
28. 4, 20 3, 34 2, 95 2, 71 2, 56 2, 44 2, 36 2, 29 2, 24 2, 19 2, 13 1, 91 1, 78 1, 72
30. 4, 17 3, 32 2, 92 2, 69 2, 53 2, 42 2, 34 2, 27 2, 21 2, 12 2, 09 1, 89 1, 76 1, 69
40. 4, 08 3, 23 2, 84 2, 61 2, 45 2, 34 2, 25 2, 18 2, 12 2, 07 2, 00 1, 79 1, 66 1, 59
4, 03 3, 18 2, 79 2, 56 2, 40 2, 29 2, 20 2, 13 2, 07 2, 02 1, 95 1, 74 1, 60 1, 52
100. 3, 94 3, 09 2, 70 2, 46 2, 30 2, 19 2, 10 2, 03 1, 97 1, 92 1, 85 1, 63 1, 48 1, 39

Таблица 3. Значения - критерия (Пирсона)

Число степеней свободы Уровень значимости
0, 95 0, 50 0, 10 0, 05
1. 0, 45 2, 71 3, 84
2. 0, 10 1, 39 4, 61 5, 99
3. 0, 35 2, 37 6, 25 7, 81
4. 0, 71 3, 36 7, 78 9, 49
5. 1, 15 4, 35 9, 24 11, 07
6. 1, 64 5, 35 10, 64 12, 59
7. 2, 17 6, 35 12, 02 14, 07
8. 2, 73 7, 34 13, 36 15, 51
9. 3, 33 8, 34 14, 68 16, 92
10. 3, 94 9, 34 15, 99 18, 31
11. 4, 57 10, 34 17, 28 19, 68
12. 5, 23 11, 34 18, 55 21, 03
13. 5, 89 12, 34 19, 81 22, 36
14. 6, 57 13, 34 21, 06 23, 68
15. 7, 26 14, 34 22, 31 25, 00
16. 7, 96 15, 34 23, 54 26, 30
17. 8, 67 16, 34 24, 77 27, 59
18. 9, 39 17, 34 25, 99 28, 87
19. 10, 12 18, 34 27, 20 30, 14
20. 10, 85 19, 34 28, 41 31, 41
21. 11, 59 20, 34 29, 62 32, 67
22. 12, 34 21, 34 30, 81 33, 92
23. 13, 09 22, 34 32, 01 35, 17
24. 13, 85 23, 34 33, 20 36, 42
25. 14, 61 24, 34 34, 38 37, 65
26. 15, 38 25, 34 35, 56 38, 89
27. 16, 15 26, 34 36, 74 40, 11
28. 16, 93 27, 34 37, 92 41, 34
29. 17, 71 28, 34 39, 09 42, 56
30. 18, 49 29, 34 40, 26 43, 77
40. 26, 51 39, 34 51, 80 55, 76
50. 34, 76 49, 33 63, 17 67, 50
60. 43, 19 59, 33 74, 40 79, 08
70. 51, 74 69, 33 85, 53 90, 53
80. 60, 39 79, 33 96, 58 101, 88
90. 69, 13 89, 33 107, 56 113, 14
100. 77, 93 99, 33 118, 50 124, 34

Таблица 4. Значения функции Лапласа


Поделиться:



Популярное:

Последнее изменение этой страницы: 2016-04-10; Просмотров: 765; Нарушение авторского права страницы


lektsia.com 2007 - 2024 год. Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав! (0.06 с.)
Главная | Случайная страница | Обратная связь