Файл: Автоматизированная система обработки и интерпретации результатов гравиметрических измерений..pdf

В зависимости от необходимой точности вычислений, от вариации {t°} u {th} и величины заданных констант l'f, 0 * , поправки

тяжести на поверхности невращающегося шара, вмещающего в себя

зенитное расстояние Солнца и Луны, т. е. угловое расстояние по

* Численные значения констант берутся пз паспортных данных прибора (например, а3 — заданная среднеквадратпческая погрешность) или опреде­

ляются предварительными просчетами.

18

часовой угол Солнца и Луны, т. е. угловое расстояние по небесному

где п — номер часового пояса, в котором производится наблюдение; % — долгота пункта наблюдения; т|с — уравнение времени, т. е. разность часовых углов среднего экваториального солнца и истин­ ного солнца; S — звездное время в 0f t мирового времени; а л — прямое

ратуру, если выполняется условие (П.4).

Поправка ô 4 за сползание нуль-пункта вводится всегда. Как показали авторы работ [15, 18, 76, 80], кривую сползания нульпункта можно аппроксимировать многочленом не выше третьей степени и использовать для этого метод наименьших квадратов.

Рассматриваемый алгоритм построен таким образом, чтобы по нему можно было обрабатывать как {R}, так и {R}, причем неза­ висимо от методики (однократной либо с повторением). Все пункты наблюдений, в которых проводились измерения в течение рейса, подразделяются на опорные, повторные и рядовые. В опорных пунк­ тах либо известно gH (при обработке рядовой сети), либо принимается g0 = 0 (в начальной опорной точке в каждом рейсе при обработке опорной сети). Программным путем (просматриваются соответству­ ющие признаки в массиве исходной информации) выбирается степень

кривой, коэффициенты которой определяются методом наименьших квадратов:

на повторных пунктах.

Решая систему линейных уравнении, которая получается из (11.12), находим коэффициенты редукционной кривой. Получив аналитическое выражение редукционной кривой для любого вычисляем искомую поправку за сползание нуль-пункта прибора.

При вычислении редукционной кривой производится браковка измерений с грубыми погрешностями. Для этого для каждой пары измерений проверяется условие

многочлены сравниваются по (11.13) и выбирается тот, который лучше удовлетворяет условию (11.13). Забракованная точка исклю­ чается из расчетов. По изложенному алгоритму проводится обработка наблюдений, по рядовым и опорным рейсам.

Итак, методы и алгоритмы первоначальной обработки в соответ­ ствии с поставленными в гл. I задачами и сформулированными требо­ ваниями к алгоритмам позволяют выполнять следующие процедуры:

1.Автоматически выбирать совокупность необходимых поправок при первоначальной обработке рейсов.

2.Осуществлять браковку грубых погрешностей на стадии обра­ ботки рейсов и при уравнивании.

3.Проводить первоначальную обработку и уравнивание при любом виде построения потоков исходных данных. При поточной обработке можно проводить расчеты по единичным массивам (рейсам),

при непоточной обработке — по совокупности единичных массивов. 4. Обрабатывать наблюдения, полученные при любой методике

съемки.

2.МЕТОД И АЛГОРИТМ УРАВНИВАНИЯ ГРАВИМЕТРИЧЕСКИХ ИЗМЕРЕНИЙ

Цель действия оператора 5 2 на множество {ЬС1, П Р П К Ь ПРПК/} состоит в том, чтобы найти множество {gH (х, у, z)} абсолютных значений

Для этого нужно произвести три преобразования: 1) уравнять

Известно, что при наблюдениях с гравиметром сумма приращений измеренных значений силы тяжести по замкнутому полигону равна не нулю, а некоторой невязке.

Задача уравнивания сводится к минимизации квадратов невязок полигонов методом наименьших квадратов [92, 109]. Гравиметристы,

20