УДК 528.11

В.А. БОНДАРЕНКО, В.И. МИЦКЕВИЧ, Н.С. СЫРОВА (Белоруссия)

Уравнивание нивелирных сетей с поиском грубых ошибок измерений

Геодезия и картография. – 2003. - №5. – с. 26-28

 

Программа NIWА2, разработанная в Полоцком государственном университете для ПЭВМ, выпол­няет уравнивание и оценку точности нивелирных сетей параметрическим способом с поиском гру­бых ошибок в измерениях. Последнее осуществ­ляется по методике, которая основана на предложении, опубликованном в [3]. Нивелирная сеть уравнивается как свободная (опирающаяся на один исходный пункт)  с анализом отношения

                                                           (1)

где vi – поправка в измерение из уравнивания, м;

                                       (2)

n — показатель степени (при n = 1,0 уравнивание выполняется по методу наименьших модулей; при n = 2,0 — по методу наименьших квадратов и т. д.). В подкоренном выражении используют­ся диагональные элементы матрицы

                                           (3)

в котором A — матрица коэффициентов парамет­рических уравнений  поправок;

Выражение (3) для корреляционной матрицы поправок справедливо в том случае, если диаго­нальные элементы матрицы весов измерении вы­числяют по  формуле

                                                (4)

где si — стандарт  измеренного  превышения, м.

За одно уравнивание исключается из обработки одно ошибочное измерение, для которого vi/di, максимально, и одновременно выполняется неравенство (1).

В исходной информации к программе NIWА2 веса нивелирных линий могут задаваться по од­ной  из  четырех  формул:

где  L0 — длина линии, км, для которой Рi´ = 1; k0 — число нивелирных  станций в линии.

Если нивелирная сеть многоклассная, то Рi´ для превышений старшего класса умножают, а Рi" делят на коэффициент а [4], определяемый по допускам инструкции. Например, при совме­стной обработке сетей нивелирования III и IV  классов

где   в скобках  даны  допуски   инструкции  на  ве­личину  свободных   членов условных  уравнений.

Программный переход от заданных Рi" к Рi´ осуществить легко. Для перехода от Рi´ к Рn предлагаем формулу

                                              (5)

где s0 — стандарт измеренного превышения, мм, для которого L0 или k0 назначалось так, чтобы Рi´ = 1. Как определить s0 для одной любой ни­велирной линии, сказано в [4, с. 350]. В край­нем случае, можно s0 назначить по h, указанной в табл. 1, с последующим применением выраже­ний s0 = hÖ`L0 или s0 = hÖ`k0. Но эта методика дает весьма приближенный результат для весов измерений.

Таблица 1

h, мм

Класс

I

II

III

IV

для L0

2,0

4.0

8,0

20,0

для k0

0,5

1

2

5

 

Применим программу NIWА2 для поиска гру­бых ошибок в измерениях для тестового примера из [1]. В табл. 2 указываются n — показатель степени; s0 — погрешность линии 17 мм и но­мера линий в порядке выбора наибольшего и недопустимого vi/di.

Таблица 2

n

1.0

1.5

2,0

3,0

s0

50

50

50

5.0

Номера линий

92

92

92

92

91

91

91

91

60

6

99

 

 

76

76

74

5

45

 

 

77

77

32

51

96

 

 

 

47

81

97

103

 

 

 

 

 

 

и т. д.

 

По данным табл. 2 видно, что поиск ошибоч­ных превышений зависит от n и s0. В табл. 3 приведены номера ошибочных линий не только при разных n, но и при разных s0. Избранный тестовый пример удобен тем, что при n = 2,0 варианты исключаемых превышений зависят от начального номера превышения: 92 или 63, для которых vi/di в начале счета практически оди­наково  и  недопустимо.

n

1,0

1,5

2,0

2,0                                       

3,0

s0

7,0

13

20

20

200

Номера линий

92

92

92

63

92

91

91

91

64

91

78

77

76

76

47

76

76

77

77

74

47

79

79

79

60

81

75

75

75

54

75

81

81

81

34

61

47

47

47

70

79

73

73

73

6

73

61

61

61

72

51

51

51

51

66

100

100

100

91

65

53

 

 

 

 

38

 

 

 

 

44

 

 

 

 

 

Сравнивая данные табл. 3 с результатами, опубликованными в [2],   можно сделать вывод, что найдены практически одни и те же ошибоч­ные линии. При этом нами в обработку взято превышение h461-462 = 2,451 м для линии 64, указанное в [1]. Отметим, что в публикации [2] принято h461-462 = 2,541 м для линии 64 и спра­ведливо утверждается, что оно содержит грубую ошибку.

СПИСОК ЛИТЕРАТУРЫ

1. Дьяков Б.Н., Рудикова М.П. О контроле, поиске и учете грубых ошибок измерений // Геодезия и картография. — 1997. — № 6. — С.  21—24.

2. Дьяков Б.Н., Рудикова М.П. Поиск гру­бых ошибок при обработке нивелирной сети // Геодезия и картография. — 1998. — № 8. — С.   13—14.

3. Коугия В.Л. Сравнение методов обнаруже­ния и идентификации ошибок измерений // Геодезия и картография. — 1998. — № 5. — С.  23—27.

4. Яковлев Н.В. Высшая геодезия. — М.: Нед­ра, 1989. — 445 с.