11
Система из уравнений общего вида в данном случае будет:
∑Y = a∑n + b∑ X + c∑ X 2
∑XY = a∑ X + b∑ X 2 + c∑ X 3
∑X 2Y = a∑ X 2 + b∑ X 3 + c∑ X 4
Для получения числовых значений данной системы составляется вспомогательная таблица.
Таблица 2.3 Вспомогательные расчеты для вычисления конкретного уравнения
параболы
|
|
|
|
|
|
|
X |
Y |
Х.Y |
Х2 |
Х2.Y |
Х3 |
Х4 |
12 |
16,00 |
192,00 |
144 |
2304,00 |
1728 |
20736 |
|
|
|
|
|
|
|
16 |
18,00 |
288,00 |
256 |
4608,00 |
4096 |
65536 |
|
|
|
|
|
|
|
20 |
20,15 |
403,00 |
400 |
8060,00 |
8000 |
160000 |
|
|
|
|
|
|
|
24 |
22,14 |
531,36 |
576 |
12752,64 |
13824 |
331776 |
|
|
|
|
|
|
|
28 |
23,48 |
657,44 |
784 |
18408,32 |
21954 |
614656 |
|
|
|
|
|
|
|
32 |
23,65 |
756,80 |
1024 |
24217,60 |
32768 |
1048576 |
|
|
|
|
|
|
|
36 |
24,62 |
886,32 |
1296 |
31907,52 |
46656 |
1679616 |
|
|
|
|
|
|
|
40 |
26,00 |
1040,00 |
1600 |
41600,00 |
64000 |
2560000 |
|
|
|
|
|
|
|
44 |
27,00 |
1188,00 |
1936 |
52272,00 |
85184 |
3747096 |
|
|
|
|
|
|
|
∑∑ 252 |
201,04 |
5942,92 |
8016 |
196130,08 |
278208 |
10228992 |
|
|
|
|
|
|
|
Полученные результаты подставляются в вышеприведенную систему из трех уравнений, после решения которых получаем следующие конкретное уравнение:
Y=8,85+0,68.Х-0,006.X2 . |
(II) |
2.3Оценка точности выравнивания
2.3.1Подготовительные работы. Проверка правильности вычислений
уравнений С целью последующего анализа результатов применения полученных
уравнений вычисляются вероятные (теоретические) значения зависимого признака по первому уравнению (HI) и второму уравнению (HII), которые сравниваются с исходными (опытными) данными (Y). Указанные сравнения
12
(∆ = Y-H) производятся по всем классам Х, а их результат возводится в квадрат (∆2) с целью устранения влияния знака. Данные вычисления показаны в табл.2.4 для прямой линии и табл. 2.5 – для параболы.
Таблица 2.4 Сравнение исходных и вероятных высот деревьев, полученных по
уравнению прямой линии
Исходные данные |
Вероятные высоты |
Отклонения, м |
|
||||
диаметр, |
высота, |
|
|
|
|
∆I2 |
∆II2 |
см |
м |
H I |
H II |
∆I |
∆ II |
||
X |
Y |
|
|
|
|
|
|
12 |
16,00 |
16,60 |
17,06 |
-0,60 |
-1,06 |
0,36 |
1,12 |
16 |
18,00 |
18,00 |
18,38 |
0,00 |
-0,38 |
0,00 |
0,14 |
20 |
20,15 |
19,40 |
19,70 |
+0,75 |
+0,45 |
0,56 |
0,20 |
24 |
22,14 |
20,80 |
21,02 |
+1,34 |
+1,12 |
1,80 |
1,26 |
28 |
23,48 |
22,20 |
22,34 |
+1,28 |
+1,14 |
1,64 |
1,30 |
32 |
23,65 |
23,60 |
23,66 |
+0,05 |
-0,01 |
0,00 |
0,00 |
36 |
24,62 |
25,00 |
24,98 |
-0,38 |
-0,36 |
0,14 |
0,13 |
40 |
26,00 |
26,40 |
26,40 |
-0,40 |
-0,40 |
0,16 |
0,16 |
44 |
27,00 |
27,80 |
27,62 |
-0,80 |
-0,62 |
0,64 |
0,39 |
|
|
|
|
|
∑∑ |
5,30 |
4,70 |
Приведенные в табл. 2.4 и 2.5 данные позволяют прежде всего проверить правильность вычислений, выполненных при получении конкретных уравнений, на предмет обнаружения грубых арифметических ошибок.
Правильность вычисления уравнений связи проверяется путем сравнения исходных значений Y с вероятными (Ŷ), полученными по уравнению I (HI) и уравнению II (HII). Критерием правильности вычислений уравнения I будет совпадение вероятных значений HI с исходными значениями Y для тех классов, в которых использованы координаты точек в качестве исходных для получения конкретного уравнения 1. В нашем примере для уравнения прямой линии значение HI равно 18,0, соответствует исходным данным Y во втором классе, то есть также 18,0. Аналогичное положение и в следующем, шестом классе: YI=23,6 практически не отличается от YI =23,65. Совпадение Y и H в остальных классах не обязательно и может наступить только случайно.
Некоторый контроль правильности уравнения II можно получить путем сопоставления Y и HI – во всех классах. В этом случае должно наблюдаться такое сочетание знаков (плюс и минус), которое отражает «срединное»
13
положение выравнивающей прямой между выравниваемыми исходными значениями Y.
Таблица 2.5 Сравнение исходных и вероятных высот деревьев, полученных по
уравнению параболы
Исходные данные |
Вероятные высоты |
Отклонения, м |
|
|||||
диаметр, |
высота, |
|
|
|
|
∆I2 |
∆II2 |
|
см |
м |
HI |
HII |
∆I |
∆ II |
|||
X |
Y |
|||||||
|
|
|
|
|
|
|||
12 |
16,00 |
16,00 |
16,15 |
0 |
-0,15 |
0 |
0,02 |
|
16 |
18,00 |
18,20 |
18,17 |
-0,20 |
-0,17 |
0,04 |
0,03 |
|
20 |
20,15 |
20,15 |
19,99 |
0 |
+0,16 |
0 |
0,03 |
|
24 |
22,14 |
21,84 |
21,60 |
+0,30 |
+0,54 |
0,09 |
0,29 |
|
28 |
23,48 |
23,27 |
23,02 |
+0,21 |
+0,46 |
0,04 |
0,21 |
|
32 |
23,65 |
24,44 |
24,23 |
-0,79 |
-0,58 |
0,62 |
0,34 |
|
36 |
24,62 |
25,35 |
25,23 |
-0,73 |
-0,61 |
0,53 |
0,37 |
|
40 |
26,00 |
26,00 |
26,02 |
0 |
-0,02 |
0 |
0,00 |
|
44 |
27,00 |
26,40 |
26,62 |
+0,60 |
+0,38 |
0,36 |
0,14 |
|
|
|
|
|
|
∑∑ |
1,68 |
1,43 |
|
О явной неправильности полученного уравнения будет свидетельствовать наличие во всех классах только +, равно как и знаков -, а также, если в нескольких начальных классах будут наблюдаться отклонения с одним и тем же знаком (+ или -), а во всех последующих классах с противоположным, а именно:
+ + + + + + +
- - - - - - -
+ + + + - - -
- - - - + + +
Заметим, что описанные критерии правильности вычислений I и II уравнений распространяются и на выравнивание по всем другим линиям связи, в том числе и на уравнение параболы 2-го порядка.
Убедившись в том, что полученные регрессионные уравнения вычислены правильно, можно приступить к оценке точности их «работы», которая может быть выполнена с учетом следующих критериев: мера выравнивая (r1); величина несовпадения исходных и вероятных значений (∆2); среднеквадратическая ошибка уравнения (mY); показатель эффективности уравнения (Fэф).
14
Приведенные критерии необходимы для выбора оптимального уравнения. Подробное описание критериев приводится в последующих разделах.
2.3.2 Мера выравнивания
Показатель меры выравнивания (r1) по Дворецкому [6] вычисляется по формуле
______________
r1 = √(∑α2 - ∑a2) / ∑α2 .
Из анализа формулы видно, что чем выше абсолютное значение r, тем точнее «работает» конкретное уравнение регрессии. Если r1 > 0,95, можно считать, что уравнение аппроксимирует опытные данные достаточно точно (надежно).
Для вычисления ri необходимо, прежде всего, получить среднеарифметическое (Ŷ) из имеющихся опытных значений Y, использовав вышеприведенную (см. раздел 1.1) формулу
N
H = i=1∑Yi / N ,
где N-число наблюдений; Yi-опытные данные отдельных наблюдений; Ŷ
– среднеарифметическое значение из всех опытных данных.
Затем следует вычислить отклонение (α) опытных данных от их среднеарифметического:
α = Ŷ - Yi,
где Ŷ - среднеарифметическое значение из всех опытных данных; Yi – опытные данные.
После этого необходимо определить величину попарных отклонений (∆) опытных и теоретических значений по формуле
N
∆ = ∑ (Yi - Hi),
i=1
где H – теоретические данные по уравнению регрессии.
Результаты расчетов отклонений показаны в табл. 2.6 для уравнения прямой линии и в табл. 2.7 – для параболы.
Подставим полученные значения и в формулу, определим показатель меры выравнивания для уравнения I прямой линии
r1 = 
(107,22 − 5,3 )/107,22 = 0,975;
15
для уравнения II прямой линии
r1 = 
(107,22 − 4,7 )/107,22 = 0,978;
для уравнения 1 параболы
r1 = 
(107,22 −1,68 )/107,22 = 0,992;
для уравнения 2 параболы
r1 = 
(107,22 −1,43 )/107,22 = 0,993.
Поскольку в нашем примере
r1 = (0,975 ... 0,993) > 0,95,
все четыре анализируемых функции аппроксимируют исследуемую взаимосвязь Д/Н достаточно надежно.
2.3.3 Величина несовпадения исходных и вероятных значений Смысл указанного критерия заключается в том, что чем меньше
отличаются опытные (Y) и теоретические (H) частоты (или иначе, чем меньше величина ∑∆2, тем точнее «работает» уравнение регрессии.
Из итогов табл. 2.4 видно, что точность выравнивания опытных данных по прямой линии для уравнения 1, полученного методом координат избранных точек, выражается величиной ∑∆I2 = 5,3, а для уравнения II, полученного методом наименьших квадратов, - величиной ∑∆II2 = 4,7.
|
|
|
|
|
|
|
|
Таблица 2.6 |
Вспомогательные расчеты для получения меры выравнивания (r1) |
||||||||
|
опытных данных по уравнению прямой линии |
|
||||||
|
|
|
|
|
|
|
|
|
Высоты, м |
|
|
Отклонения, м |
|
||||
Y |
H |
α |
|
α 2 |
|
∆ |
∆2 |
|
(опытные) |
(теорети- |
|
|
|
|
|
|
|
|
ческие) |
|Y-22,34| |
|
|
|
|
Y-H |
|
|
|
УРАВНЕНИЕ I |
|
|
||||
16,00 |
16,60 |
6,34 |
|
40,20 |
|
-0,60 |
0,36 |
|
18,00 |
18,00 |
4,34 |
|
18,84 |
|
0,00 |
0,00 |
|
20,15 |
19,40 |
2,19 |
|
4,80 |
|
+0,75 |
0,56 |
|
22,14 |
20,80 |
0,20 |
|
0,40 |
|
+1,34 |
1,80 |
|
23,48 |
22,20 |
1,14 |
|
1,30 |
|
+1,28 |
1,64 |
|
23,65 |
23,60 |
1,31 |
|
1,72 |
|
+0,05 |
0,00 |
|
24,62 |
25,00 |
2,28 |
|
5,20 |
|
-0,38 |
0,14 |
|
26,00 |
26,40 |
3,66 |
|
13,40 |
|
-0,40 |
0,16 |
|
27,00 |
27,80 |
4,66 |
|
21,72 |
|
-0,80 |
0,64 |
|
∑201,04 |
|
|
|
∑107,22 |
|
|
∑5,30 |
|
|
|
|
|
= |
|
|
|
|
|
|
|
r1 |
(107,22 |
− 5,3 |
)/107,22 |
= 0,975 > 0,95 |
|
Ŷ=201,04/9=22,34 Вывод: уравнение аппроксимирует опытные данные достаточно точно