Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Разумов, О. С. Пространственная геодезическая векторная сеть

.pdf
Скачиваний:
5
Добавлен:
19.10.2023
Размер:
5.69 Mб
Скачать

В итоге, искомая ошибка mw оказывается равной

т%- = а2 cos2 ф А-

Ьгsin2 ср=

„2 ѵс

cos2

ßa~ ß

^cosec2

- Р а -

р— -|-

 

+

sec_1£- sin2

 

 

(2.37)

 

,

о

Эг

 

 

 

 

Для более общего

случая,

когда щф щ ,

задачу

также решим

графоаналитическим путем.

На рис. 35 изображена вспомогательная сфера единичного ра­ диуса, построенная вокруг точки 1 хорды. На этом же рисунке показаны направление хорды 12, плоскость синхронизации Р (в плоскости чертежа), единичные векторы щ и я2, направленные на спутник, и линии положения Д0, и La2 конечной точки век­ тора W, определяемые следами пересечения сферы плоскостями, перпендикулярными векторами щ и я2.

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

воздействием

этих

ошибок

(перпендикулярных

к

плоскости

рисунка)

меняется

угловое

положение плоскости

синхронизации

и точка

W может

изменить

свое положение на

W

(изменение

направлений

векторов ö,_jв

плоскости синхронизации

не влияет

на направление вектора №). Отклонение вектора ДИ7 по направ­ лению хорды при этом будет равно

ДЯ7

(2.38)

 

sin (ßi -f ß.)

Перейдя по общим правилам теории ошибок от конечных при­ ращений к средним квадратическим ошибкам, найдем

ти? = к

sin2 ß2

2 S'n°- ßi

(2.39)

sin2 (ßLAr ß2)

- sin2 (ßi + ßo)

 

 

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

стояния спутника. Так, в [38] показано, что в среднем

 

ца =

,

(2.40)

 

cos г

 

где цо— ошибка определения координат в зените. При равноточных измерениях на станциях

т 'ѵ = sin2 (ßi + ß.) ^

+ sin2 ^

( 2-41>

Результаты, получаемые по формулам (2.41) и (2.37), одно­ значны, но для анализа функции тлѵ при предельных значениях углов ß формула (2.37) более удобна. Уравнение (2.41) легко при­ водится к виду

т\ѵ =

+ 4),

(2.42)

90

]

где /'j и /'2 — расстояния до спутника от наблюдательных станций, L — длина хорды.

Остановимся на особенностях найденных соотношений.

1. При параллактическом угле ßc, равном 90°, и равноточных измерениях на станциях погрешность направления вектора W отражается на сфере круговой ошибкой радиуса р. Если и вторая пара синхронных наблюдений выполнена с той же точностью, а плоскости синхронизации перпендикулярны друг к другу, то ошибка направления искомого вектора ё будет равна значению р.

2. С геометрической точки зрения объект наблюдений (ИСЗ) выгоднее всего располагать на хорде или ее продолжении. Хоть

это и невозможно

осуществить на

практике, но в этом случае

ошибка_ /птг была , бы равна малой

полуоси эллипса ошибок век­

тора W, предельное_минималыюе значение которой может дости­

гать величины р /|/

2. И хотя при этом большая полуось эллипса

стремится к бесконечности, это не влияет на точность искомого направления хорды, так как смещение векторов W в плоскости, нормальной к вектору ё, не влияет на угловое положение этого вектора.

3. Влияние несинхронности наблюдений на ошибку углового положения вектора W уменьшается по мере уменьшения угла между траекторией спутника и плоскостью синхронизации.

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

Во-первых, наблюдать низкие спутники опасно из-за искажаю­ щего влияния рефракции (полагают, что предельное допустимое значение г0=75°); во-вторых, при заданной величине z0 область расположения подспутниковых точек, доступных наблюдению, резко снижается с уменьшением высоты спутника над . поверх­ ностью Земли; в-третьих, при небольшой высоте Н спутника на­ блюдать его будет трудно из-за большой угловой скорости дви­ жения.

Итак, для практических расчетов погрешностей направления хорды, при необходимом минимуме измерений, могут быть реко­ мендованы формулы (2.37), (2.41) и (2.32), а для неравноточных измерений на станциях — формулы (2.39) и (2.32). Эти соотноше­ ния в самом общем виде решают задачу оценки точности.

В недавней работе Г. А. Устинова [75] предложены более про­ стые формулы для оценки точности направления хорды для част­ ного случая равноточных измерений, когда углы ßi = ß2=ß, и плос­ кости синхронизации расположены симметрично относительно хорды. Вывод формул Г. А. Устинова можно повторить на основе

более общих

закономерностей (2.32) и (2.37), если положить

ф = 0°, it = 90°

; /2 = 90°+ ■— .

При условиях, поставленных Г. А. Устиновым,

т\Ѵі = т\г„ =

2|і2 sin2 ß

Ц2

=

тлѵ

sin2 2ß

2 cos2 ß

 

 

 

91

и при значениях

ßc== l8 0 ° -2 ß ,

тогда ni\v

/

ßc = 90° - ß,

2

В2

2

В свою очередь, учитывая значения углов і, на основании (2.32) получим

„ о

У

О

О V

 

 

 

 

 

 

 

 

 

 

2"1\Vcos“ V

LI“

COS------

 

 

 

 

fl“

 

 

 

 

 

гщ, =

 

h

2

 

 

 

 

 

 

=

mii’

(2.43)

sin2 7

. ,

Pc

7

 

,

. J e

. о V

 

 

sm -

----- sin-

 

4 sin2

—— sm- —

 

 

 

 

 

 

9

'

 

 

 

2

 

2

 

 

 

 

л У . n

I

j.1- sin2

Y

 

 

 

 

 

 

 

 

2mjp sm2 —7"

 

 

 

 

fl“

 

 

 

mx cos3

=-.

 

 

ßC .

 

 

 

 

=

mi

 

sin- у

 

sin-5

 

. .

ßc

Y

 

 

 

 

-----sin-1 7

 

4 sin2

----- cos- —

 

 

 

 

 

 

2

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.44)

 

m - =■■m- +

тл2 cos- -ф=

 

2mw

 

 

 

B"

 

 

(2.45)

 

 

 

 

 

sm2 у

 

sm2

ßc

 

 

 

 

 

 

 

 

 

 

 

—— sin- Y

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

По мнению Г. А. Устинова, формулы (2.43) и (2.44) должны отражать горизонтальную та и вертикальную гпп ошибки направ­ ления хорды. Однако пользоваться этими формулами можно толь­ ко при соблюдении трех условий:

а) равноточности измерений на станциях; б) равенстве расстоянии от наблюдателя до ИСЗ (ПВЦ);

в) симметричном (зеркальном) расположении плоскостей син­

хронизации относительно геоцентрического сечения

хорды.

§ 10. УРАВНИВАНИЕ МНОГОКРАТНЫХ НАБЛЮДЕНИЙ

 

ПРИ ОПРЕДЕЛЕНИИ НАПРАВЛЕНИЯ ХОРДЫ

 

Сп о с о б Ве й с а

 

Среди нескольких известных методов определения вероятней­ шего значения направления хорды из серии многократных наблю­ дений наибольшее распространение в работах Смитсонианской астрофизической лаборатории (SAO) получил способ Г. Вейса [71]. В этом методе решение трехмерной задачи о пересечении плос­ костей синхронизации сведено к двухмерной задаче пересечения линий — следов этих плоскостей — в плоскости со, перпендикуляр­ ной к хорде (рис. 36).

92

При этом предусматривается следующий порядок действий,

1.По приближенным координатам точек 1 и 2 хорды опреде­

ляется единичный вектор этого направления

Т3=

.

(2.46}

IR t - R i I

2.Составляется уравнение плоскости геоцентрического сече­ ния хорды

(R Wi Яа) = АX + BY + CZ = 0

(2.47}

и определяется угол наклона этой плоскости к экватору

cos і0

С

(2.48)

Л2 + В"- + О '

/

 

Рис.

36. Геометрия

плоскости

син­

Рис. 37. Следы

пересечения пло­

 

 

хронизации

 

 

 

 

скостей синхронизации с пло­

 

 

 

 

 

 

 

 

 

скостью со

 

3.

Вводится

специальная

система координат TLT2T3 с нача­

лом в точке 2 и направлениями осей, определяемыми взаимнопер-

пендпкулярнымн единичными векторами

 

 

 

 

Т3;

 

Тг =

і 1Х *3

и

Т\ = Тя х Т г .

 

(2.49}

 

 

 

 

 

 

I R i

X ^2 I

 

 

 

В этой системе

координатная

плоскость (Ти Т2) = со перпенди­

кулярна к хорде.

 

 

 

момента

0j синхронных наблюдений опреде­

4.

Для

каждого

ляется вектор нормальной к плоскости синхронизации

 

 

 

 

 

 

 

U-7/ = К X a2)j

 

(2.50}

и его направляющие косинусы

 

 

 

 

 

cos aw. =

wx .

;

cos ß,,., =

W„

cos yw

Wz .

(2.51).

 

I wf I

W,-

zj

 

;

 

 

 

 

 

Wi I

 

93-

5. Определяются углы между векторами Wj и осью 7\

cos а

wI

 

sina,-=

Wi

(2.52)

——

-Ti,

п j

 

IV/

1

 

1

1IV; 1

 

и отклонения векторов Wj от перпендикуляра к ветору Тг при­ ближенного направления хорды

~ еі =

(2.53)

6. Полагая, что все плоскости синхронизации при их нача ном положении проходят через точку 1 и пересекают координатную плоскость со по линиям, отстоящим от точки 2 на удалении ejL12 (рис. 37), составляют уравнения этих линий

cos а

+

sin a,jT2 — — e;L12,

(2.54)

где а — угол между положительным

направлением

вектора Wj

я осью Ті.

 

 

на L12, переходят от линей­

Разделив обе части равенства (2.54)

ных величин к угловым

 

 

 

 

cosg^Tx+

sina^a = — ej,

(2.55)

и составляют для каждой пары синхронных наблюдений уравнений поправок

cosajTx + sina^Tj +

ej =

е^,

/'-=1,2, . . ., п

(2.56)

•с весами

 

4L],

 

 

Pj =

 

 

(2.57)

9

9 , 9

9

 

ИГ

П 4- Pi

 

 

где гі и r%— удаление ИСЗ от наблюдательных станций' 1 и 2. Значение весов здесь согласовано с уравнением (2.39), выражаю­ щем ошибку углового положения вектора Wj по направлению хорды.

7. Решив систему уравнения (2.56) по способу наименьш квадратов, получают поправки ті и т2 к приближенному направ­ лению хорды и вычисляют поправки 6 к начальным значениям ее направляющих косинусов

/$ а \

/Т'іх Тіу

ТіД

/ ті

\

(2.58)

= Т 2Х Т 2у

T 2Z

т3

.

\ А /

\ т з< Т 3у

Т 3 2 ) \ о )

 

Поправки ті и тг можно считать приближенно поправками к азимуту и вертикальному углу хорды.

8. Оценку точности окончательных результатов выполняю помощью обратной матрицы коэффициентов нормальных уравнений

М2

p2Q =

(2.59)

ah

 

 

Э4

где ошибка единицы веса

[ре2]

(2.60)

п —2

 

а іпа и nih — суть горизонтальная и вертикальная составляющая ошибки направления хорды.

Для отыскания средних квадратических ошибок и mvcos ij)

направляющих углов хорды тензор [2.59] необходимо подвергнуть преобразованию

фА

а Н

“ f W ,

(2.61)

\ Шак

m l)

 

 

 

 

где

 

sin Aq

COS Ag \

 

Ат.

=

(2.62)

 

 

COSy40 -sin^0j ’

Ао—-угол отстояния плоскости геоцентрического сечения хорды от

плоскости

часового

круга,

проведенного

через

хорду,

л

cos г0

,

 

 

 

 

sm А0=

-----1

 

 

 

 

 

и тогда

cos tJj

 

 

 

 

 

 

"4 = т2аsin2 А 0 -f ml cos2 А 0 -f 2 sin А 0 cos А 0 mh,1.

 

 

(2.63)

 

т2\ cos2 ф =-• m2acos2 A0-j- ml sin2 A 0 — 2 sin A 0 cos A üm ah

 

Метод Вейса обычно применяют тогда, когда в дальнейшей математической обработке пространственной сети используются значения направляющих косинусов хорды. Для непосредственногоотыскания сферических координат ф и Л хорды из серии синхрон­ ных наблюдений более удобны другие методы.

Сп о с о б И . Д. Ж о и г л о в и ч а

В основу этого способа [5], [19] легло уравнение компланар­ ности векторов в треугольнике, образованном наблюдательными станциями и ПВЦ

COS фcos Л

COS фsin Л

sin Ф

F — а2 Е) — COS l \ COS 4

COS 6Xsin tx

sin ö;

COS 6 2 COS to

COS Ö2 sin 4

sin 6;

или

F= cos ф cos ба cos б2 [ tg ф sin (t2— 4) Н-

+tg 6Хsin (Л — 4 )-f- tg б2 sin (^i — Л)} = 0.

(2.64>

(2.65)

Получив тем или иным

способом приближенные значения ф

и Л направляющих углов

хорды, от условного уравнения (2.65)

можно перейти методом линеаризации к уравнению поправок функ-

95-

ции ряда измеренных величин и составить такие уравнения для каждой пары синхронных наблюдений

( 2. 66)

где

с 5чВ0^°'^і1ыв ІІлен уравнения, получаемый после подстановки в (2.64) приближенных значений направляющих углов хорды и из­

меренных топоцентрических координат ПВЦ, />, = —

и в свою

mrF

 

очередь

 

 

 

 

(2.67)

где Л4бі t

корреляционная

матрица

непосредственных нзме>

рений.

 

 

 

При независимых результатах измерений имеем

mF =

dF_

 

 

2

 

 

(2.68)

 

dt,

 

Частные производные в (2.66)

и (2.67)

равны

(ihjr) =

sin ^ C0S sin s‘n

~ О + sin 'Фcos 4 sin öx sin (4 — Л) +

-г cos я|з cos бх cos fi3 sin (4 — 4),

=cos яр cos Öx sin 62 COS (4 — Л) -f- COS яр cos ö2 sin 6XCOS (4 — Л),

~ S’n 4 COS ‘Ф sin 63 sin (4 — Л) — sin cos 62 sin 1)3 sin

(4 — 4 )

— cos <4cos яр cos ö3sin (4 — Л),

(2.69)

sin 0|) COS fix sin Ö2 sin (4 — 4) — COS яр cos öx cos ö3sin (A.— 4) + + cos яр sin 6Xsin ö2sin (4 — Л),

cos яр sin 8, cos 62COS (4 —1Л) -j- sin яр COS 8, COS Ö2COS (4 — 4),

= COS t cos 63 sin öx cos (4 — Л) — sin яр cos 6Xcos ö2 cos (4 — 4).

96

Практически,

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

функций

В. В. Бойков и II. Я. Плешаков в [5] рекомендуют составить мат­

рицу _частных

производных

направляющих

косинусов

векторов

я,, 5-2, Ё по углам ф, Л, 4, 6і, 4, 62

dcos yE

 

 

 

дcos аЕ

dcos

\

 

 

дф

Ö<|>

дф

 

 

 

d cos аЕ

dcos ߣ

dcos yE

 

 

 

ЗЛ

ö.Y

дЛ

 

 

 

 

дcos аа

dcos ßa

dcos у

 

 

 

 

дбі

dbl

Ö6i

 

 

 

 

дcos а„

dcos ßön

dcos y0i

 

 

 

 

а1

dli

dt1

 

 

 

 

Ö/i

 

 

 

 

дcos а Ял

dcosßac

dcos yai

 

 

 

<56о

döt

öS,

 

 

 

 

дcos а 0>

dcos ß0^

dcos v0>

 

 

 

д(2

dt0

dk

 

 

 

sin ф cos А

— sin ф sin A

cos ф sin А

COS ф cos Л

sin бх cos 4

— sin 6Xsin 4

COS бх sin ti

cos öx COS 4

sin б2 COS 4

— sin ö2sin 4

COS Öjlsin 4

cos 6.2 cos 4

COS ф \

О'

COS 61

О(2.70)

cos б, О ’ /

Для вычисления значений производных функций F по перемен­ ным ф, Л, б;, 4 достаточно заменить в определителе (2.64) соот­ ветствующую строку строкой матрицы (2.70) и вычислить этот определитель.

Вычисление весов по (2.67) и (2.68) требует громоздких рас­ четов; но если корреляционная матрица непосредственных измере­ ний, полученная после астрометрической обработки снимков, имеет диагональный вид, а сами ошибки топоцентрических коор­ динат ПВЦ равны между собой, так что іт = т а cos6= p,, то для получения весов уравнений поправок можно предложить более простое решение.

Дело в том, что веса Pj должны характеризовать собой точ­ ность определения свободных членов уравнений поправок. В рас­ сматриваемом случае определитель (2.64) формально равен объе-

му_параллелепипеда, построенного на

единичных векторах äi,

5%

и Е. Эту величину можно найти по формуле

 

 

7 = (ßI Xfl2)/£'o =

scosffii,

(2.71)

где s = sin(ßi + ß2) — площадь

параллелограмма,

построенного

на

единичных векторах аі и а9,

w- ■угол

между хордой и вектором

Wj.

 

 

 

 

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

97

Дифференцируя (2.71) и переходя к конечным приращен ням, получим

 

AF =

— s sin wAw + cos wAs.

(2.72)

Имея в

виду, что угол w при доброкачественных

наблюде­

ниях всегда

близок к

прямому, второй член уравнения (2.72)

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

квадратическим ошибкам, можно ограничиться учетом

только

первого члена этого уравнения

 

піг — s2mlv = sin2 (ß2 -f- ßo) miy/.

(2.73)

Подставляя в это уравнение формулу (2.39),

окончательно най­

дем

р2 sin2 ßt

(2.74)

m2F ^ р2 sin2 ßo +

и

 

 

Pj = ------------- —с --------- .

(2.75)

(р? sin2 ß-2 -I- F

sin2 ß1)/.

 

В [5] аналогичное выражение получено более сложным путем. После совместного решения уравнении (2.66) по способу наимень­ ших квадратов определяются поправки к приближенным значе­ ниям направляющих углов хорды. Точность окончательных ре­ зультатов характеризуется тензором ошибок

Ліф.ѵ = f Q-

(2.76)

В геодезической литературе известно несколько модификации способа И. Д. Жонгловнча. В некоторых из них [32], [99] вместо условия компланарности, выражаемого равенствами (2.64) или (2.65), иногда приравнивают нулю выражение, заключенное в фигурных скобках уравнения (2.65), пли, как сделано в [99], при­ равнивают нулю определитель

cos A0

sin A0

lg Фо

 

cos А

sin А

tg 6j = 0.

(2.77)

cos 4

sin t.

tg&

 

Это уравнение получено путем сокращения на постоянные множители направляющих косинусов единичных векторов в (2.64).

Лианеризируя (2.77) можно составить уравнение поправок функций ряда измеренных величин вида

а/Л + b/Jq -I-

Ej,

(2.78)

 

—sin A0

cos A0

0

 

ai ~

cos А

sin А

 

tgÖ!

(2.79)

cos 4

sin 4

 

tg62

 

 

 

 

 

Ьі = COS2 '1>о

sin (/о —А)-

(2.80)

 

 

 

 

 

98

Уравнения (2.78) решают относительно поправок пд и по

способу наименьших квадратов.

Здесь, однако, следует заметить, что подобные «упрощения» основного способа не всегда будут целесообразны, а в некоторых случаях могут привести к нежелательным последствиям. Дело в

том,

что

коэффициенты уравнений (2.78) содержат тангенсы углов

ф , 6 і

и

6 2 и могут быть получены с большими погрешностями.

И если угол ф принимает значения, близкие к л/2 в редких слу­ чаях, то топоцеитрические склонения 5і и 62 будут близки к 90° каждый раз, когда ПВЦ наблюдается на фоне зоны близполюсных звезд. Между тем при использовании для целей уравнивания основного соотношения (2.64) эти опасения не возникают.

При обработке результатов наблюдений по способу И. Д.- Жонгловнча или Г. Вейса отдельно не рассматривается влияние оши­ бок регистрации моментов наблюдения на окончательные резуль­ таты. Здесь предполагается, что при современной точности ра­ боты служб времени основное влияние на итоговую погрешность направления хорды оказывают ошибки определения топоцентрическнх координат ПВЦ. Однако если ошибкам регистрации мо­ ментов наблюдений на станциях свойствен систематический ха­ рактер, то они могут привести к нарушению синхронности наблю­ дений и к искажению окончательного результата. Методику учета таких ошибок рассмотрим при описании еще одного оригинального приема уравнивания наблюдений, предложенного И. Вяйсяля

[ 102].

Спо с о б В я й с я л я

Этот, пожалуй, наиболее простой способ уравнивания наблю­ дении рассмотрим с небольшими изменениями, внесенными в него автором [64].

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

1. Определив тем или иным способом приближенные значения направляющих углов фо и Л0 хорды, вычисляют по (2.50) и (2.51) для каждого момента синхронных наблюдений значения направ­ ляющих косинусов векторов Wj, нормальных к плоскостям синхро­

низации.

_

2. Вычисляются углы между векторами Wj и приближенным

направлением хорды

cos Wj =

cos aw. cos ф0 cos A0 -|- cos ß^. cos Ф0 sin Л0 -f- cos yw sin ф0.

 

(2.81)

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

4* 99

Соседние файлы в папке книги из ГПНТБ