Закатов Вища геодезія 1
.pdfдЕ _ |
1 |
дМ |
1— е •cos Е * |
|
|
|
|
|
|
|
|
|
|
|
|
дМ |
= т - т 0, |
|
|
|
|
|
|
|
|
||||
в итоге |
|
|
|
|
|
|
|
|
дп |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
96 |
_ |
sin і •cos (со+ 6і)•sin Ф (Т —То) |
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
Подобннм же образом |
дп |
|
|
|
|
cos 6 -sin і? (1— |
е • cos Е) |
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
- f - = |
|
----— sin І sin ((О + |
'0') = |
sinx |
COS (со + '&) |
|
(OJ -f- 0) = |
sin І |
COS U |
2 |
|||||||||||||||
|
де |
------- r |
де * |
||||||||||||||||||||||
де |
|
cos 6 |
де |
|
4 |
|
|
1 |
' |
|
cos 6 |
|
|
4 |
1 |
' |
v |
r |
|
cos 6 |
|||||
если обозначить u — co + |
|
ft и учитьівать, что dO = |
dco. |
|
|
|
|
|
|||||||||||||||||
|
Из |
(110.30) |
получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
1 |
90 |
|
|
. |
E 9 т / 1 + e |
, і / 1+ e д . |
E |
|
|
||||||||||||
|
|
|
COS2 А |
|
|
*8— -вГ |
І7 1 = 7 |
+ |
|
K T = 7 |
ИГХ2 — |
|
|
||||||||||||
|
|
|
96і |
0 |
о |
О . |
|
Е |
1 + Є |
-- |
|
COS" |
|
|
1 + е |
|
1 |
|
9Е |
|
|||||
|
|
|
9е |
2 cos2-т- tg |
|
---- л |
|
|
|
|
1 — е |
|
„ Е |
|
де |
|
|||||||||
|
|
|
|
|
2 & 2 / |
і _ е2 |
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
~ cos2 — |
|
|
|
|
Так |
как |
|
|
9 |
|
|
|
|
|
|
|
|
|
1+■ е |
|
|
1 |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
і / 1+ |
Є. = ] / . |
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
де |
|
V |
1- е |
|
|
|
1— е |
|
1- е 2 ’ |
|
|
|
|
|
|||||
|
|
|
|
|
|
дМ |
|
|
дЕ |
sin Е —е cos Е ■4 ^ = ° , |
’ |
|
|
|
|
||||||||||
|
|
|
|
|
|
де |
|
|
|
де |
|
|
|
|
|
|
|
|
де |
|
|
|
|
|
|
|
|
|
|
|
дЕ |
|
|
е cos Е) ==sin Е, |
дЕ |
__ |
sin Е |
|
|
|
|
|
|||||||||
ТО |
|
|
|
|
де (1- |
|
|
|
|
4 |
а ппс тг |
|
|
|
|
||||||||||
90’ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
sin Е |
|||
|
^ c o s + tg + |
j Z |
+ |
f + |
|
. H |
- / ^ |
|
COS" |
|
|
||||||||||||||
|
де |
|
|
|
cos2 E |
1— |
e cos E |
||||||||||||||||||
|
|
|
|
|
|||||||||||||||||||||
и окончательно |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
О |
|
б1 |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
/ Т |
+ |
Г |
" |
|
|
|
|
||
96 |
|
|
sin і cos u |
In |
0 |
|
|
• |
0 |
|
І |
|
|
|
cos T |
cos T |
sin E |
|
|||||||
де |
|
|
cos 6 |
,2 c o s^ -sm -r T I I 7 |
|
> |
1 |
|
|
£ |
|
E |
l |
—e cos E |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
cos - y cos - y |
|
|
|
|||
|
|
|
|
sin І COS U |
I |
|
. |
a |
|
1 |
|
|
s m _0 |
|
cos —0 |
sin E |
|
|
|
|
|||||
|
|
|
|
|
cos 6 |
|
|
sin гг |
1 — e2 |
|
* |
. |
E |
|
E |
1 —e cos E |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
w |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s m _ |
|
cos — |
|
|
|
|
|
|
||
|
|
|
|
|
|
sin u cos u sin |
. { _ ! _ |
+ |
____ 1____ |
|
|
|
|
||||||||||||
Далее |
|
|
|
|
|
|
|
|
cos 6 |
|
|
\ l |
— e2 |
~ |
1 —e cos E |
j |
|
|
|
|
|||||
|
|
|
96 |
|
|
cos іsin u |
|
96 |
|
96 |
sin і cos u |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
9i |
|
|
|
cos 6 |
|
|
d(0 |
|
/9Л |
|
|
|
|
|
|
|
||||
56 |
|
96 |
90 |
9.E |
9M |
|
_ |
sin іcos u |
sinO |
|
1 |
|
^ _ |
sin і cos u sin ■6’ |
|||||||||||
дМо |
|
90 |
9E |
9M |
9M0 |
|
|
cos 6 |
|
sin E 1 —e cos E |
|
sin E (1— e cos E) * |
|||||||||||||
|
|
|
9a _ |
9a |
90 |
дЕ |
|
dM |
_ |
________ cos і__________ sin 0 |
T — T 0 _ |
||||||||||||||
|
|
|
dn |
90 |
дЕ |
дМ |
|
дп |
|
[1 + |
cos2 і tg2 u] cos2 u sin E 1— e cos E |
|
490
|
|
|
________ cos і sin |
___________ T — T o |
|
|||||||||
|
|
|
[1+ |
cos2 і tg2 u] cos2 u sin E |
1 - е |
cos E * |
||||||||
|
|
da |
da |
5'fr |
|
|
|
|
|
|
|
|
||
HO |
|
де |
50 |
де |
|
[1 + |
cos2 і tg2 u\ cos2 e |
де |
||||||
|
|
<90 |
|
- |
a |
/ |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
1 — e cos Er) |
|
|||||||
|
|
|
де |
|
sin 0 |
(T |
|
|
1 |
|||||
|
da |
|
|
cos і sin |
» ______ |
1 |
|
|
||||||
|
де |
(1 + cos2 і |
|
( — L _ 4 |
-e cos E |
|||||||||
|
tg2 u) |
cos2 |
і |
\ 1 —e'1 |
||||||||||
|
|
|
|
|
|
|
|
da |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5£2- = 1, |
|
|
|
|||
|
|
|
|
|
da |
|
|
sin І tg І |
|
|
|
|||
|
|
|
|
|
di |
|
1 + cos2 і tg2 u |
|
|
|
||||
|
|
|
da |
|
|
da |
|
|
|
|
cos і |
|
|
|
|
|
|
5co |
|
50 |
|
(1 + |
cos2 і tg2 u) cos2 u |
’ |
|||||
da |
da |
50 |
dE |
|
dM |
|
|
|
|
|
cos і sin ' |
|||
дМo |
50 |
dE |
dM |
|
dM0 |
|
(1 + cos2 і tg2 u) {1 - е cos E) cos2 e sin E ’ |
|||||||
|
|
dRdn |
= (1 |
|
є cos E) |
dnda |
+, a dn5 |
(i —e eos E), |
||||||
|
|
|
|
|
|
|
da |
|
|
|
2a |
|
|
|
|
|
|
|
|
|
|
dn |
|
|
|
n |
|
|
|
d |
= |
, A |
T?\ |
dn |
(1— e cos E) = |
т . e
dR dn
<9 /л |
r?\ |
dE |
dM |
e sin E (T —To) |
(1 —£ cos Zi) |
dM |
On |
1 - е cos E |
|
2a (1— e cos E) |
|
ae sin E (T —T 0) |
||
3n |
' |
1 - е cos E |
|
d R |
5 , |
T?\ |
<9 / |
j-,, |
a cos E + ae sin E |
|
— |
= - d r ( a — aecosE) =:z- t e ( - aecosE) |
|||||
|
dE
(M + e sin E) = sin E + e cos E
де
dE де
|
|
dE |
|
sin E |
|
|
|
|
|
|
де |
1 — e cos E |
* |
|
|
||
|
dR |
= — a cos E |
ae sin2 E |
|
|
|||
|
де |
|
|
1— e cos E |
|
|
||
dR |
dR |
dE |
dM |
|
|
1 |
|
|
dM r |
dE |
dM |
dM0 |
ae sin E 1 — e cos E |
|
|||
Уравнения (110.31) можно записать как |
|
|
|
|||||
( 56 |
56 |
56 |
56 |
56 |
56 1 ( dn 1 |
|
||
dn |
де |
5£2 |
di |
5(0 |
дМ0 |
de |
|
|
da |
da |
da |
da |
da |
da |
dQ |
Ddu. (110.32) |
|
dn |
де |
5Q |
di |
5(0 |
5Mo |
di |
||
|
||||||||
dR |
dR |
dR |
dR |
dR |
dR |
5(0 |
|
|
. dn |
де |
5£2 |
di |
5(0 |
dMo |
d M 0 |
|
491
IIодставив (110.32) в (110.28), получим
-рд- 'Q - Q - D-du — Q dRM. |
(110.33) |
Если наблюдения производить со станций с известньїми координатами» то
Q - Q - D - d u . |
(110.34) |
В зтом случае измереннне по крайней мере с двух твердих станций значе |
|
ння , а ' , г' дают возможность определить поправки du |
к елементам прибли- |
женной орбитн и 0. Исправив атими поправками приближенное значение орбитн,
вичислим (см. 110.25) координати R c на моменти наблюдений ИСЗ сопределяемьіх станций. Таким образом, измерения с определяемнх станций будут произведенн на известнне положення ИСЗ, что дает возможность вичислить коорди нати определяемнх станций. Если спутник наблюдался с нескольких твердих
станций или наблюдалось более двух положений ИСЗ с одной станции, то du определяется по способу наименьших квадратов.
В аллиптическую орбиту можно ввести вековне возмущения, внзваннне сжатием Земли и сопротивлением атмосфери. Сжатие Земли приводит к вращению линии узлов в акваториальной плоскости (регрессия Q ) H K вращению
линии апсид по плоскости орбитн (регрессия со).
Ранее уже били полученн формули, определяющие мери атих движений:
Я == — 10 |
) /2 • cos і, |
|
\ Г Ср / |
to == — 5 |
/2' (5 sin2і — 1). |
Влияние атмосфери внзнвает изменение форми и размеров орбитн (а и е) в ее плоскости при очень малом изменении перигея. Вращение атмосфери внзн вает вековне изменения в наклонении, которне можно считать линейннми на коротком интервале времени. Если взять только вековую часть возмущений» то орбита может бить представлена в общем виде степенннм рядом
и — uQ-\-du-\-du (Т — TQ)-\- -Tf-du (Т — T0)2-f . . ., |
(110.35) |
где
du = d (п, е, Я, і, со, М 0)
(110.36)
du = d(n, е, Я, і, w, М 0)
492
Число членов зтого уравнения зависит от величини Т — Т 0. Для определения орбитш второго приближения используем метод, аналогичннй предьідущему, но уравнение (108.53) заменяют на интеграл
|
|
|
|
|
|
т |
|
|
|
|
|
|
|
Л/ = А/„ + |
j ndT. |
|
|
(110.37) |
|||
|
|
|
|
|
|
о |
|
|
|
|
В соответствии с новим определением орбитш (как орбитн второго при- |
||||||||||
ближения) нужно искать поправки du, du, |
du и т. д. |
|
|
|||||||
Пусть |
Т — Т 0 достаточно мало, чтобн ограничиться первшм членом раз- |
|||||||||
ложсния, |
т. е. |
|
|
|
|
|
|
|
|
|
|
|
|
и — UQ-j—du —|—du (Т — Уд). |
|
(110.38) |
|||||
Тогда искомьіми параметрами будут вектори |
|
|
||||||||
|
|
|
du — (dn, |
de, |
dQ, |
di, |
doo, |
dM0) |
|
(110.39) |
|
|
|
du — (dn, |
de, |
dQ, |
di, |
dco, |
dMQ) |
|
|
|
|
|
|
|
||||||
її уравнение (110.35) примет вид |
|
|
|
|
|
|
||||
|
|
|
db |
\ |
|
|
( du |
|
|
|
|
|
|
cos da |
= {Q, |
|
|
(110.40) |
|||
|
|
|
|
|
|
|
||||
|
|
|
dR |
|
|
\du-y |
|
|||
|
|
|
У |
|
|
|
|
|||
|
|
|
R |
|
|
|
|
|
|
|
где матрица D равна |
|
|
|
|
|
|
|
|
||
|
' |
дд |
дЬ |
дЬ |
дЬ |
|
дд |
дд |
|
|
|
|
ди |
де |
dQ |
di |
|
<9(0 |
дМ о |
|
|
|
D |
да |
да |
да |
да |
|
да |
да |
( T - T 0) |
(110.41) |
|
ди |
де |
dQ |
ді |
|
<9(0 |
дМ0 |
|||
|
|
|
|
|
||||||
|
|
dR |
dR |
dR |
dR |
|
dR |
dR |
|
|
|
. |
ди |
де |
dQ |
ді |
|
д(о |
дМ0 ) |
|
|
В итоге получим |
|
|
|
|
|
|
|
|
||
|
|
|
|
Q '{& |
D) |
|
Q dRM. |
(110.42) |
Если учитнвать не только вековне, но и периодические изменения злементов орбитьі, возникающие из-за несферичности гравитационного поля и изме нения плотности атмосфери с внсотой, то орбиту следует записивать как
и Ьи — щ-\- du + du (Т— У0) + • • ч |
(110.43) |
•где Ьи — короткопериодические изменения орбитн.
493
Приведем вираження для короткопериодических изменений четьірех злементов, которне бьіли полученьї Козаи:
б; = ~ |
2sin і jcos 2 |
+ |
to) + |
е cos (Ф + 2со) + |
|
^ |
|
|
|
cos (Зг*)+ |
2со)| |
|
|
|
|
бй = — J |
cos * |
— Щ + е sin ^ |
sin 2 (й + со)— |
|
|||
---- Y е sin (&-f 2со)---- е sin (ЗФ+ 2(о)j |
|
|
|||||
б R = - ^ - J ~ ~ У Р |
( і — si n2 і ) |
1 — |
( і — 1/~1 — е 2 ) c o s,0' + |
|
|||
+ |
. г І + J |
(-J-)а Р s i n 2 1 [ т cos3(# + “>] |
(110.44) |
||||
|
|||||||
би = / ( ~ - ) 2 |( 2 |
— ^ sin2 |
ф - М - \ - е sin -0) -f- ^1 — |
sin2 |
x |
|||
X |j-|-e ( l ---- 1----- V 1 |
) sin ft + |
|
—] / 1 —e2 ) sin 20J — |
||||
— ------ 1- sin2 |
e sin (0 -j- 2co)— |
------ sin2 ij sin 2 (0 + |
to) — |
|
|||
|
---- g- (1 — sin2і) e sin (30 -f 2oo)j |
|
|
||||
где |
|
|
|
|
|
|
|
|
ц = |
|
( l - - i s i n 2i ) l / T ^ " . |
(110.45) |
Из сказанного следует, что применение орбитального метода требует достаточно сложньїх вичислений. Кроме того, орбитальньїй метод менее точен, чем метод синхронних наблюдений, из-за неточности знання параметров, определяющих движение ИСЗ, или законов изменения зтих параметров во времени (гравитационное поле Земли, плотность атмосфери и т. д.). Короче говоря, основним препятствием для широкого применения орбитального метода для геодезических целей является трудность високоточного прогнозирования орбитн на большие отрезки времени Т — Т 0. В то же время применение зтого метода очень заманчиво, так как он не требует синхронних (или квазисинхронннх) наблюдений.
§ 111. Использование наблюдений ИСЗ для определения параметров, характеризующих гравитационное поле Земли
Если считать Землю симметричной относительно оси вращения, торазложение V земного потенциала по сферическим функциям можно записать в следующем виде:
F = Ж [ і + С20 ( 4 - ) 3 Р т+ Сш (£■)* і ?40 + |
(1И.1) |
Если при зтом ограничиться определением только 6*20, с 30, С40, то, как показано И. Д. Жонголовичем, формули для возмущений узла и перигея за один драконический период будут иметь вид
6 Q — р2/ + |
Рзб'зо |
р 4П + Р22>^2 |
бсо = q2J + |
^3 6 * 3 0 |
(111.2) |
q±D -f- q22J 2 |
494
где
|
|
|
|
|
|
|
J |
— |
|
____ r |
|
) |
|
|
|
|
|
||
|
|
|
|
|
|
|
J — |
|
2 Ь2° |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
n |
— |
35 |
c |
I |
|
|
|
|
|||
|
|
|
|
|
|
|
|
u |
g |
^40 J |
|
|
|
|
|||||
И |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
p2 — |
|
2яйд cos і |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
a2(1- е 2)2 » |
|
|
|
|
|||||||
|
|
|
|
Рз |
|
2лаі cos і |
p |
|
45 |
sin2 і |
3 1 |
e sin ft) |
|
||||||
|
|
|
|
|
a2 (1 —e2)3 |
|_“8 |
2 J |
|
sin і |
|
|
||||||||
|
|
|
2"°or-°si |
г (4 -sm2 1- 4 - ) |
f 1 + - 5 - 4 |
- |
|
|
|||||||||||
|
|
|
a4 (1 — Є2)4 |
L W |
|
|
|
|
‘ / |
\ |
|
Z |
J |
|
|
|
|||
|
|
|
— |
|
sin2і ---- e2cos 2(0J |
|
|
|
|
|
|
||||||||
P2: |
|
2ла?cos |
[ ( - r |
sin! j - |
T |
) |
|
+ |
( " f - sini J ~ |
T |
) e 008 “ + |
||||||||
|
a 4 ( l _ |
e2)4 |
|
||||||||||||||||
|
|
-f- |
sin2 і + |
|
|
e2 — ( “I" sin2 * — Y |
) e2 cos 2coJ |
|
|
||||||||||
|
|
|
?2 |
|
23X02 |
[~2 — |
^ |
sin2 і |
1 |
|
|
|
|
||||||
|
|
|
|
|
\2 |
L |
|
|
|
|
|
J |
|
|
|
|
|||
|
|
|
|
|
a2 ( 1 - е 2)2 |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
= |
2лп3 |
Г / |
3 |
. o . , |
15 |
|
|
|
|
|
3 |
|
105 |
. » . |
, |
|||
|
^ ao |
|
( --- 5- sm21+ -7T—sm |
|
|
|
-- ------g -S m 3>+ |
||||||||||||
4z |
|
a3(1— e2)3 *- ' |
|
|
|
|
|
|
in4i) + ( |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
105 |
in4 i) |
|
|
|
sm to |
|
|
|
|
|
||||
|
|
|
|
|
|
|
, 21 — |
|
і |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
} |
|
|
J |
sin |
|
|
|
|
|
||
|
|
|
2л |
|
Г |
12 |
93 . |
. . . |
21 . |
4 A |
|
|
|||||||
|
|
Qi |
a 4 (1 •— e2)4 |
|
— |
« |
|
sm , + |
T |
sm |
v |
|
|
||||||
+ (-^j- sin2І---- 4 sin4 i'j |
COS2ft> + |
|
|
------ j - |
sin3І-f- 4 4 sin4 2) ’e24- |
||||||||||||||
|
|
+ ( _ TT+ |
|
15 |
smйі |
|
■Щ- sin4 іj e2cos 2(0J |
|
|
||||||||||
|
Я22 |
2ла; |
^ |
— 2 + |
sin2і ----1- sin4 |
e-1 COS (0+ |
|
||||||||||||
|
a4 (1_ K»2)4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
95 |
2 . |
;445 |
sin4 г -j- ( — 2 4- |
|
|
sin2i + |
sin4i j |
cos 2(0 -f |
||||||||||
+ (l12 |
sirr l |
48 |
|
|
|||||||||||||||
+ ( |
25 |
sin2і ------- sin4 |
|
e cos (o + |
^ |
~ |
|
sin2 |
X |
||||||||||
|
|
|
|||||||||||||||||
X e cos l |
|
~ 4 ~ s i n 4 ~~'W sinU) e2+ |
( і 2 ~ Ж |
sin2i + |
|||||||||||||||
|
|
3t0+ ( l 2 |
|||||||||||||||||
|
|
|
|
|
|
45 |
sin4 |
|
e2cos 2oo~I |
|
|
|
|
|
|||||
|
|
|
|
|
|
16 |
|
|
|
|
|
|
|
|
|
|
|
|
|
(111.3)
(111.4)
495
Величини 6Q и бю в вираженнях для козффициентов р и q заданн для мо менте прохождения спутником восходящего узла. Значение большой полуоси а определяется из уравнений
TQ — Tn — J T |
[ З ------|- s i n 2 *— е cos ш (1 — 5 sin 2 і) + |
|
-h~Y~ е2 ( і — |
sin2 + “|"е2 cos2оо ( і — sin2 i^j J |
, (111.5) |
Тп = TQ -f /_^е2 ("^ -)2 (2 — ^ -sin 2i) ( l — 2ecosco + - |- e 2 cos 2to) |
|
|
которне решаются методом последовательннх приближений. Здесь |
Га и Тп — |
соответственно наблюденнне значення драконического и аномалистического периода. Величини 6Q и бсо внчисляются по результатам наблюдений а, б и г. Значення Q и со в узле и перигее определяют интерполяцией (по Лагранжу или по способу наименьших квадратов) по неравномерно расположенннм во
времени значенням Q и со, которне определенн по результатам наблюдений. |
|
Далее по формулам (111.4) можно определить козффициентн р 2, |
р 3, р 4, р22 |
и g2, q3, q4, g22 по средним значенням а. е, со, и, і. Из решения |
уравнений |
типа (111.2) определяют неизвестнне / , С30 и D и по ним три первнх параметра |
гравитационного потенциала Земли С20, С30, С40. Сжатие Земли может бить
определено через |
/ , to, а, /, М как |
|
|
“ = 1 + Т п + Т |
Л - Т Jn — п%+ |
( т р + Ж 'Р п + Ж - Jn*+ -ЩІГ ) ’ |
|
|
|
|
( 111.6) |
где |
|
|
|
|
п |
(о2а3 |
(111.7) |
|
~ ж |
||
|
|
|
Практически в формуле (111.6) во всех членах, кроме первого, можно вместо./ подставить его приближенное значение 0,001625. Положив также
(0 |
2зт |
729 212 • 10-10 рад/с2, |
|
86164,09 |
|
||
|
|
|
|
|
fM — 398 600 км3/с2, |
|
|
получим |
|
|
|
а = J + |
1728,88 • 10~6= 0,003354 ^ |
. |
|
Для зллипсоида Красовского а — onQ1 . |
|
||
|
|
2Уо о |
|
Определение внсших гармоник в разложении потенциала по сферическим |
функциям — задача достаточно сложная, и изложение ее внходит за рамки данной глави
496
§ 112. Связь различньїх геодезических систем с помощью ИСЗ
Геодезическая прямоугольная система координат R связана с геодезическими координатами В, L, Н г уравнениями
X = (N + Нг) •cos В •cos L )
Y = (N Нг)‘ cos В • sin L . |
(112 1 > |
|
Z = |
[N (1 —е2) -f- Я г] - sin В ) |
|
Система координат В |
задается координатами В°, |
L 0 и Н гс начального» |
пункта и параметрами расчетного аллипсоида. Параллельность между геоцентрической системой (общий земной зллипсоид) и геодезической (референцзллипсоид) обеспечивается не только применением астрономических коорди
нат ф, X и азимута а, но также и внесением поправки за уклонение |
отвеса |
в астрономические координатні при условии, что астрономический |
азимут |
тоже исправлен за уравнение Лапласа. Пусть ф, X и На будут координати исходной точки геодезической системи, полученние астрономическими мето
дами и нивелированием, и |
— абсолютнне уклонения отвеса. Тогда, |
|
В = ф — І |
|
Ь = Х— Ц • cosec ф |
|
(112.2) |
А — а — (X— £)• sin ф
Нт= На- і
Погрешности в ф Д и Я практически не влияют на непараллельность геодезической системи относительно геоцентрической. Ошибка же в азимуте dA,. измеренном из начального пункта системи, внразится в повороте системи координат вокруг вертикали на величину dA. Так как ось, вокруг которой поворачивается система координат, не исходит из начальной точки, то одновременно будет иметь место и перенос. Матрица поворота может бить получена. как произведение пяти матриц поворота:
|
|
|
|
|
|
|
sin dA |
М, |
|
|
|
|
|
|
cos dA |
|
СОвф |
|
|
Sin ф/ |
|
|
0 |
0 |
— COS ф\ |
|
/f |
cos X |
sin X |
O' |
(112.3> |
1 |
0 |
' |
|
— sin X |
cos X |
0 |
|
0 |
'\ |
0 |
0 |
1 |
|
||
sin ф/ |
|
|
Ошибки в І и т] приведут егце к двум вращениям геодезической системи координат относительно геоцентрической, причем матрицн поворота будут равнн M s и М ц, если считать положительннм направлением запад для d Н. и север для d\і :
|
1 |
0 |
— dH cos X0S |
М р = \ |
0 |
1 |
— dl sin X0 |
\d^ • cos к0 |
dl • sin Хо |
1 |
|
0 |
— dv • cos ф0 |
(112.4). |
|
— drI • sin ф0*sin Xo' |
|||
Mrj = [ dr\ • cos ф0 |
|
l |
dr\ • sin ф0cos Xo |
\dr\ • sin ф0 |
— dij • sin ф0• cos X, |
1 |
^2 п. C. Закатов |
497' |
Общий поворот от совместного влияния d A , d |
dr\ будет |
М аіч — Ма • |
(112.5) |
Внполнив все перемноження и учитнвая лишь членьї 1-го порядка малости, ЗІолучим
/«11 |
« 1 2 |
«1з \ |
М аіч =~ 1 «21 |
«22 |
«23 1 > |
'«31 |
«32 |
«33' |
где
«11 == 1, «21 = - — d A • sin ф + dr\ ■COS ф |
|
«зі = d A • COS ф • sin J^ |
■cos X 4" d*\ • sin ф- sin X |
«із = dA • sin ф-— dt\ • cos ф |
|
«33 = |
1 |
#32 = |
— dA • COS ф • sin l + d l • sin X — dr\ • COS ф • cos X і. |
|
& CO II |
! |
• COS ф • sin X - d l • cos X — dr\ • sin ф • sin X |
«23 = СІА • COS ф . s i n ;k - d l :sin X 4~ • sin ф • cos X
«33 — 1
(112.6)
(112.7)
Обозначив через С — (СІ5 С2, С3) координати центра масо Земли в гео-
дезической системе координат R' = (Xr, Y ' , Z') и через R = (X , Y, Z) — геоцентрические координати, получим
R' = МАІп ■(R + С) = R + |
С+ ( М - Е) - {ЇЇ+С) . |
(112.8) |
|
Произведя умножение М — £ на й |
+ С и сгруппировав подобнне члени |
||
относительно d A , |
йц, получим |
7 |
|
X' = X -j- ^ |
-f [sin cp • (Y + C 2) — cos cp • sin X• (Z + C 3)\ dA — |
|
|
— cos X-(Z-\- C 3) • dІ — [cos ф • (Y + C2) + sin cp • sin X-{ZA- C3)] • dr\ |
|
||
Y* = Y -)- C2-f [—sin ф• (X -f C-L)-f cos ф• cos X• (Z + C3)] dA — |
|
||
— [sin X’ (Z-{- C3)] • d%+ [cos ф • (X + Сг) -j- sin ф• cos X- {Z-\- C3)] • dц |
. # (112.9) |
||
Z' = Z -f- C3-(- [cos ф • sin X (X + Сг)— cos ф• cos X• (Y -f C2)] dA + |
' |
||
[cos X■(X -j- Ci) -4- sin X"{Y -j- C2)] ’ d \ 4“ [sin ф• sin X • (X 4~ Ci) — |
|
||
или |
— sin ф • cos X • (Y |
C2)] Лц |
|
|
|
|
|
|
R' = R-\-C + |
G |
( 112. 10) |
Если имеются не связаннне между собой триангуляции (для каждой свой расчетннй зллипсоид и своя начальная точка), с которнх синхронно наблюдались положення ИСЗ, то при различии геодезических координат ИСЗ, полу-
498
ченннх в разннх системах (пока не учитьіваются ошибки измерений), геоцентрические координати в разннх системах должнн бнть одинаковьі, т. е.
|
|
fdAA |
ІR■2= + |
С2 + |
/dAo\ |
І |
||
Ri |
+ Сх-\- G\ |
(ЩхІ |
G2 1^£2 |
|||||
или |
|
\dr]1J |
|
|
\dr}2/ |
|
||
|
|
|
/d A x\ |
|
/dAo\ |
|
|
|
|
|
|
|
|
|
|
||
R 2- R x = £c + Gx(dgi \ - G J d l 2 |
I. |
(112.11)> |
||||||
|
|
|
|
\dv)x/ |
|
\dr\2J |
|
|
Имея определение |
двух |
положений ИСЗ, можно |
отнскать вектор Ас и |
|||||
d A x d \ x dx\, если на |
dA, d |
dr\ наложить какое-нибудь |
условие, например' |
|||||
или |
|
dA2 = d^2= dr\2—0 |
|
|
(112.12) |
|||
|
|
d^)X— d'E.2, |
dy\x — d^]2. |
(112.13) |
||||
dAx — dA2, |
Если имеется больше двух одновременньїх наблюдений ИСЗ, то параметри Асі5 Дс2, Дс3, d A x, d | 1? dy] х получают по способу наименьших квадратов.
Знание зтих параметров и дает возможность связать между собой различнне триангуляции. Совместная обработка наблюдений ИСЗ с ряда различннх триангуляций дает возможность получить довольно хорошеє приближение к общему зллипсоиду своеобразннм осреднением референц-зллипсоидов раз личннх триангуляций.
32*