УДК 528.11
В.А. БОНДАРЕНКО, В.И. МИЦКЕВИЧ, Н.С. СЫРОВА (Белоруссия)
Программа 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 с.