Файл: Разумов, О. С. Пространственная геодезическая векторная сеть.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 29.10.2024

Просмотров: 54

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

ряда опорных пунктов, заданных в системе координат, отнесенной

кцентру масс и осп вращения Земли.

2.Наблюдая с течением времени эволюцию орбиты ИСЗ, в

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

численные характеристики внешнего

гравитационного поля

Земли

и параметры верхних слоев атмосферы.

Для решения этой задачи

 

 

 

 

дополнительно

необходима

 

 

 

 

детально

разработанная

тео­

 

 

 

 

рия движения ИСЗ.

 

 

 

 

 

 

 

 

3.

 

 

 

 

 

 

 

 

 

 

 

 

дачи по определению поло­

 

 

 

 

жения пунктов земной по­

 

 

 

 

верхности

в системе

коорди­

 

 

 

 

нат,

отнесенной

к

 

центру

 

 

 

 

масс

и осп вращения

Земли,

 

 

 

 

связано

с

предположениями,

 

 

 

 

что орбита

 

спутника,

вычис­

 

 

 

 

ленная

по

наблюдениям

с

 

 

 

 

твердых станций

наблюдения,

 

 

 

 

достаточно

точно аппроксими­

Рис.

12. Взаимное

положение

ИСЗ

рует

его движение

в некото­

 

и наблюдательной станции

 

рый

начальный

момент

вре­

ИСЗ

разработана

 

 

мени; что теория движения

хорошо и все необходимые

физические

пара­

метры (плотность воздуха, гравитационные

постоянные и т.

д.)

из­

вестны. Тогда па любой момент времени можно определять

поло­

жение ИСЗ в геоцентрической системе координат и спутник

пред­

станет как некий подвижной опорный пункт.

 

наблюдательных

Наблюдая такой спутник

с

тех

пли

иных

станций возможно

находить их положение в той же системе коор­

динат '.

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

Если же положение станций необходимо получить в некоторой геодезической системе координат, элементы ориентирования кото­ рой в теле Земли известны, то основное векторное уравнение с уче­ том (1,5) примет вид

ß0A * X ' = 7 —? — бг0.

(1.61)

1 Естественно предположить, что если результаты наблюдении на станции содержат только угловые координаты (а' и б') спутника, то для определения положения точки потребуются, как минимум, две серии наблюдении в моменты Л и І2, а при измерении одних топонеитрнческих расстояний г' — три серии.

30


Совокупное решение всех трех перечисленных задач входит в комплекс работ так называемого динамического метода использо­ вания наблюдений спутников Земли, и это решение может быть достигнуто только способом последовательных приближений.

В настоящее время для привязки отдельных островов, экспе­ диционных пунктов и в исследовательских целях используют так называемый «полудпиамическип метод коротких дуг» [40], [73] или

орбитальный метод, в котором прогно­

 

 

зирование положений ИСЗ осуществ­

 

 

ляется

на

небольшом

отрезке

траек­

 

 

тории,

в пределах одного или

двух обо­

 

 

ротов

спутника,

между участками, ох­

 

 

ваченными

наблюдениями

с опорных

 

 

пунктов на земной поверхности.

экстра­

 

 

Траектория ИСЗ

на участке

 

 

поляции вычисляется здесь по

элемен­

 

 

там оскулирующего возмущенного эл­

 

 

липса,

параметры

которого

определя­

 

 

ют численным интегрированием урав­

 

 

нений движения, или по уравнениям эм­

 

 

пирических орбит.

 

 

 

стан­

 

 

Положение

наблюдательных

 

 

ций

отыскивают

в этом

методе

путем

Рис. 13. Геометрия с'ии-

прямой

реализации

векторного

уравне­

хроппых

наблюдений

ния

(1.60)

или

отдельных

его

моди­

 

ИСЗ

фикаций.

Другую возможность геодезического (геометрического) исполь­ зования спутников доставляют синхронные наблюдения ИСЗ. Если, например, наблюдать спутник одновременно с двух точек земной поверхности (рис. 13), то взаимное положение этих пунктов, опре­

деляемое геодезическим

топоцентрическим вектором

A/?i2= j6i2,

можно получить из уравнения

 

ARn

- R , ~ R ! =~г\ — ~гІ

(1-62)

В этом случае спутник уже играет роль только вспомогатель­ ной визирной цели и его точное положение на момент наблюдений знать не обязательно. Уравнение (1-62) и его модификации поло­ жены в основу различных методов построения пространственных опорных геодезических сетей.

В 1969 г. Л. Б. Закиров, используя теорию «бликоидиого» пре­

образования пространства [21],

предложил

обобщающее

уравне­

ние,

связывающее все векторы трехгранника

(см.

рис. 13)

 

 

(/•[ щ /о) г =■ г'оЯі + г\ Rn — 2r\r'2p0cos

,

(1.63)

где

ßc — угол между векторами

и г'2, а р о— единичный

вектор,

направленный по биссектрисе угла ßc-

 

 

 

31


Учитывая, однако,

что

последний член равенства

(1.63) яв­

ляется функцией измеренных величин и равен

 

2г\ г 2р0cos- у = г\ г[ + г\ Го,

 

уравнение Закирова приобретает форму тождества

 

(г{ + Го) г =

Го (R1+ гі’) -f г\ {Я, + Го).

(1.64)

Для случая, когда с одной станции наблюдалось два положе­

ния спутника, тождество сохраняется

 

(/'і + Го) R = Го (і\ — fl) + Г\ (/‘о — Го).

(1.65)

При математической обработке наблюдении ИСЗ для решения геодезических задач основой для составления уравнений поправок измеренных величин служат дифференциальные формулы измене­ ния топоцентрических координат ИСЗ под влиянием тех или иных возмущающих факторов.

Структура таких формул, основанных на функциональной зави­ симости топоцентрических координат от взаимного расположения наблюдательной станции и ИСЗ, показана ниже. Согласно рис. 12, имеем

(г'Т = (Хс - Х,„Г -f (Yc -

Ymf 4- (Zc - zmf

 

 

tgö'

Zc

Zm

 

 

V(Xc-X,nY + iX c -Y m)

I

( 1. 66)

 

 

V -

V

 

 

 

tgt' =

X*Cc- Xл тm

t = а — Ѳ, t' — а! Ѳ.

Дифференцируя эти уравнения, получим

cos 6' cos t'

cos 6' sin t'

sin 6

sin 6' cos t'

sin 6' sin t'

cos 6'

r

 

 

r

 

sin t'

 

cos t'

0

r' cos 6'

r' cos 6'

 

 

/ dXcI

\

f dXm

= U\

dYc

 

I— U\ dYn

\ d i c )

\d Z n

В то же время, на основании

(1.7),

имеем

I dXc dX„ dYc - dY„ \ d Z c - dZ„

(1.67)

COS0

sin0

° \

/ dEc ^,

f dEc

 

—sin Ѳ

cos 0

°

dHc

= S0 dHc

( 1.68)

0

0

\d Z c j1

\d Z c

 

1/

 

32


Поэтому

(1.69)

Дифференциальные формулы (1.114), связывающие изменение геоцентрических координат Нс, Нс, Zc спутника с изменениями параметров начальной промежуточной орбиты и других возмущаю­ щих воздействий, даны на стр. 43. Если qt — матрицу частных производных геоцентрических координат ИСЗ по элементам оскулирующей орбиты, G — матрицу частных производных элементов оскулирующей орбиты по элементам начальной промежуточной ор­ биты, F — матрицы частных производных элементов оскулирующей орбиты по коэффициентам разложения гравитационного потен­ циала Земли и некоторым другим параметрам подставим в (1.69), то получим

dr'\

dö' I = US0q,G ( dt' )/

(1.70)

Это уравнение показывает в самом общем виде структуру урав­ нений поправок наблюдений в динамическом методе, когда в каче­ стве искомых неизвестных фигурируют поправки к элементам на­ чальной промежуточной орбиты, поправки к физическим парамет­ рам Земли и поправки к положениям наблюдательных станций в геоцентрической системе координат.

Если при решении этой задачи исходные координаты наблюда­ тельных станций были заданы в системе координат OrXrYrZr, отне­ сенной к центру референц-эллипсоида, и взаимное положение этих станций полагать заданным достаточно надежно, то в соответст­ вии с (1.5), в уравнения поправок (1.70) в качестве дополнитель­ ных неизвестных можно включить и элементы внутреннего ориен­ тирования референц-эллипсоида в теле Земли, полагая величины б*о, бг/о, 6z0, a.ij едиными для всей системы пунктов, принадлежа­ щих данной системе относимости.

Очевидно, что для определения всех неизвестных параметров уравнения (1.70) необходимы массовые наблюдения ИСЗ с наблю­ дательных станций, равномерно расположенных по земному шару и специальные программы наблюдений, в которых тот или иной параметр определяется наиболее надежно.

2 Разумов О. С.

33

Описанию и анализу динамического метода посвящена обшир­ ная литература; особенно глубоко и детально он рассмотрен в ра­ ботах У. Каула [27] н [71].

§ 6. НЕКОТОРЫЕ ВОПРОСЫ ОЦЕНКИ ТОЧНОСТИ ПОЛОЖЕНИЯ ПУНКТОВ

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

О ш и б к и п о л о ж е н и я п у н к т о в

Позиция наблюдаемого объекта в пространстве определяется пересечением трех поверхностей положения.

Так как полученные в процессе измерения скалярные величи­ ны и содержат ошибки, и соответствующие нм поверхности поло­ жения испытывают колебания относительно своей вероятнейшей позиции, то точность положения объекта в пространстве зависит как от этих ошибок, так и от взаимного расположения пересекаю­ щихся поверхностей положения.

По этой же причине, если число выполненных измерений пре­ вышает минимально необходимое, то поверхности положения не пересекаются в одной точке, а образуют вокруг искомого пункта

некоторую

фигуру

погрешностей.

 

 

 

 

 

 

Математическая обработка таких измерений преследует цели

отыскания

наилучших

 

приближений

к

неизвестным

значениям

функций измеренных

величин и решения

вопросов

оценки

точно­

сти. Обе эти задачи

решаются

методом

наименьших

квадратов,

который

может реализоваться

путем

совместного

решения

урав­

нений поправок s к измеренным величинам вида

 

 

 

(

диі

о,-.

duj

 

 

 

 

■+

// =

в/,

(1.71)

\

дхс

дуі

О

 

 

 

 

 

 

 

 

 

 

 

под условием

 

 

[/же] = min.

 

 

 

(1.72)

 

 

 

 

 

 

 

 

Записав

предыдущие

уравнения

в матричной форме, получим

 

 

 

 

 

A X +

L =

e

\

 

 

 

(1.73)

 

 

 

 

 

етРе = min

J

 

 

 

 

 

 

 

 

 

 

 

 

Здесь А — матрица

частных производных

скалярных

функций и

по искомым величинам

(координатам

пунктов) в счислимой точке,

34