Файл: Разумов, О. С. Пространственная геодезическая векторная сеть.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