· как вычислить веса величин, используемых в обработке;
· как выполнять оценку точности параметров преобразования и преобразованных координат, так как Гаусс разработал метод оценки точности для другой целевой функции.
К этому же классу задач относится и методика анализа деформаций колец тоннеля. Современные алгоритмы обработки результатов измерений предусматривают вычисление положения вероятнейшей окружности под условием минимума суммы квадратов расхождений реального положения колец тоннеля от вероятнейшего. Новая целевая функция не позволяет использовать при обработке результатов измерений все точностные характеристики измеренных и приближенно известных величин и выполнить объективную оценку как результатов измерений, так и их функций.
С развитием дальномерной техники в геодезии стали широко применять полигонометрию и линейно-угловые сети. При уравнивании таких сетей возникли трудности: как вычислять веса для угловых и линейных измерений? Вес – величина размерная или безразмерная? Этот вопрос станет понятен, если целевую функцию (5) записать в виде:
,(6)где
– веса угловых измерений; – поправки в угловые измерения; – веса линейных измерений; – поправки в линейные измерения.Если вес – величина безразмерная, то в целевой функции (6) будут складываться, например, квадратные секунды с квадратными миллиметрами. В результате остро стоит вопрос о соотношении весов в угловых и линейных измерениях. Для того чтобы устранить возникшую парадоксальную ситуацию, можно записать формулу (6) в следующем виде:
. (7)Средняя квадратическая ошибка единицы веса μ2, стоящая перед знаком суммы, не влияет на отыскание минимума, следовательно, целевую функцию (7) можно представить в окончательном виде:
. (8)По сути, эта та же формула Гаусса, лишь записана она в другом виде. Однако такая форма записи снимает все трудности поиска соотношения весов между разнородными измерениями, так как под знаком суммы стоят безразмерные коэффициенты, если средние квадратические ошибки и поправки вычислены в единой размерности. Более того, целевая функция (8) позволяет вычислять поправки в любые величины, которые измерены или известны приближенно при совместной их обработке. В целевой функции (8) роль веса выполняет величина, обратная квадрату средней квадратической ошибки, что и рекомендовал Гаусс. Вводить в эту целевую функцию понятие веса бессмысленно, так как при обработке результатов измерений это ничего нового не прибавит и не убавит.
Аналогичная ситуация сложилась и при анализе результатов наблюдений за деформациями колец туннеля. По результатам угловых и линейных измерений вычисляют координаты нескольких точек по периметру тоннеля в условной системе координат, а затем вычисляют положение аппроксимирующей окружности под условием:
, (9)где Δi– отклонение радиуса аппроксимирующей окружности от реального расстояния от оси тоннеля до обделки.
Целевая функция (9) позволяет вычислить интересующие параметры тоннеля, но не допускает учета точностных характеристик измеренных величин, и невозможно воспользоваться алгоритмом Гаусса для оценки точности вычисленных параметров тоннеля. К сожалению, целевая функция (9) нашла неоправданно широкое применение при решении инженерно-геодезических задач, несмотря на ее откровенные недостатки. Автором разработана методика обработки результатов измерений при наблюдениях за деформациями колец тоннеля с использованием целевой функции (8) при сохранении возможности выполнения оценки точности всех вычисляемых параметров тоннеля.
Обычно метод решения задачи в геодезии состоит из самостоятельных этапов.
Этап 1. Формулировка основной цели работы.
Выполнить анализ деформаций колец тоннеля c заданной средней квадратической ошибкой 3 мм. Под термином "деформация колец тоннеля" может подразумеваться: отклонение размеров тоннеля от проектного; отклонение размеров тоннеля от вероятнейшей окружности.
При анализе отклонений размеров тоннеля от проектного значения все проектные размеры при обработке входят как константы, и к ним не требуется вычислять поправки. В зависимости от поставленной задачи могут встречаться оба варианта анализа деформаций. В большинстве случаев в инженерно-геодезической практике задается несколько точностных характеристик, например, допуск на радиус тоннеля и допуск на отклонение от вероятнейшего радиуса тоннеля. В таком случае целесообразнее вычислять действительный размер собранного тоннеля. В дальнейшем примере будем рассматривать именно этот вариант, как наиболее характерный.
Этап 2. Выбор метода измерений, который решает поставленную задачу.
Одновременно с выбором метода измерений необходимо записать математические зависимости между измеряемыми и вычисляемыми величинами (в данном случае деформационные характеристики колец тоннеля). Строгая математическая зависимость между измеряемыми и вычисляемыми величинами полностью исключает дальнейший выбор каких-либо иных "независимых параметров". В том случае, если измеряемые и вычисляемые величины связаны нелинейными уравнениями, то приведение данной функции к линейному виду возможно лишь в том случае, если удастся найти приближенные значения именно вычисляемых величин, и в таком случае нет места другим "независимым параметрам".
Этап 3. Предварительная оценка точности с использованием метода наименьших квадратов, по результатам оценки точности выбор метода измерений, а также обоснование точности полевых измерений.
Этап 4. Полевые измерения.
Этап 5. Обработка результатов полевых измерений, вычисление уравненных значений искомых величин с оценкой их точности.
Процесс измерения заключается в следующем. В некоторой точке А устанавливают инструмент и измеряют углы наклона βi, и расстояние Si до стенок тоннеля в нескольких точках, расположенных в вертикальной плоскости, перпендикулярной оси тоннеля. Зная проектные размеры тоннеля и выполнив дополнительные измерения можно определить приближенные координаты оси тоннеля относительно оси теодолита (рис.3) со средней квадратической ошибкой 3 - 4 см.Используя полярные координаты Si и βi и их точностные характеристики, необходимо вычислить положение оси тоннеля, радиус тоннеля и деформационные характеристики тоннеля с объективной оценкой точности. Как видно из рис.4, уравнения, которые связывают измерения и интересующие нас величины, имеют вид:
,(10)гдеR – радиус тоннеля; Δi– отклонение фактического положения стенок тоннеля от окружности; Si - расстояние от прибора до наблюдаемой точки; X – расстояние от прибора до центра тоннеля по оси Х; Y– расстояние от центра тоннеля до горизонтальной оси прибора по оси Y; φi – угол между направлением на центр тоннеля и наблюдаемой точкой. Учитывая, что
; (11) ,(12)где α – угол между горизонтом инструмента и направлением на центр тоннеля;
, если ,(13)где βi – измеренный угол между горизонтом инструмента и визирным лучом на точку Iи
, если . (14)При этом необходимо выбрать знак координат X и Y. В дальнейшем будем считать величину Y положительной, если центр прибора расположен ниже оси тоннеля, X – величиной положительной, если центр прибора расположен слева от оси тоннеля, как показано на рис.3.
Уравнению (10) будут удовлетворять лишь уравненные значения, причем измеренные или приближенно известные величины (далее выделены их волнистой чертой сверху) и уравненные связаны следующими равенствами:
.С учетом этих представлений приведем уравнение (10) к линейному виду относительно поправок в измеренные величины, но в начале определим зависимость между поправками в βi и φi. С учетом (12) из уравнений (13) и (14) получим:
, если ; (15) , если . (16)В свою очередь, поправку Vγ получим из уравнения (11), представив его в виде:
,(17)