Материал: 2655

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

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 Вывод: уравнение аппроксимирует опытные данные достаточно точно

Источник: https://studfile.net/preview/16412687/