Глава V
ДИФФЕРЕНЦИАЛЬНЫЕ ФОРМУЛЫ
§ 33. Дифференциальные формулы первого рода
После обработки триангуляции и вычисления геодезических координат пунктов может оказаться, что начальные данные (длина и азимут исход
ной стороны, координаты начального пункта), принятые при обработке, под лежат небольшим изменениям. Это, естественно, вызывает необходимость испра
вления всех вычисленных широт, долгот и азимутов триангуляции. :Конечно,
можно заново решить треугольники и вычислить широты, долготы и азимуты
р
|
А |
Рис. 63 |
Рис. 64 |
на основе новых исходных данных, однако проще исправить координаты пунк
тов путем вычисления поправок к ним.
Формулы, выражающие поправки геодезических координат пунктов и ази-· мутов направлений за изменение начальных данных триан гуляции, называются д и ф ф е р е н ц и а л ь н ы м и ф о р м у л а м и п е р
в ого род а.
Бывают случаи, когда необходимо изменить параметры принятого рефе
ренц-эллипсоида. Это может случиться при использовании старых триангу
ляций, которые относились в России к эллипсоидам Вальбека, :Кларка, Бесселя,
тогда как в настоящее время в СССР принят эллипсоид :Красовского. :Кроме того, в связи с переходом к эллипсоиду :Красовского возникает необходимость
перевычисления координат пунктов на этот эллипсоид, так как в старых ката
логах приведены координаты пунктов, вычисленные с использованием пара
метров эллипсоида Бесселя, который был принят в геодезических работах
СССР до 1942 г. :Конечно, координаты, отнесенные к новому эллипсоиду, могут
быть получены путем перевычисления координат пунктов с использованием новых значений параметров референц-эллипсоида. Однако и в этом случае проще
получить новые координаты пунктов путем вычисления и введения поправок
за изменение параметров эллипсоида.
Формулы, выражающие поправки геодезических координат з а и з м е -
н е н и е п а р а м е т р о в э л л и п с о и д а, |
называются д и ф ф е р е н - |
ц и аль н ы ми ф о р мул а l\I и в то рог о |
род а. |
150
Выведем дифференциальные формулы первого рода.
Пусть в результате ранее выполненных вычислений получены геодезические
координаты конечных точек стороны триангуляции АВ (см. рис. 62), ее длина,
:прямой и обратный азимуты.
Обозначим: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
В 1 , |
L 1 |
- |
координаты |
пункта |
А; |
|
|
|
|
|
|
|
||||||
В,., |
L 2 |
- |
координаты |
пункта |
В; |
|
|
|
|
|
|
|
||||||
А |
1 |
2 |
- |
азимут |
с |
А |
на |
В; |
|
|
|
|
|
|
|
|
||
А |
2 · 1 |
- |
азимут |
с |
В |
на |
А ; |
|
|
|
|
|
|
|
|
|||
|
· s - |
расстояние АВ. |
|
|
|
|
|
|
|
|
|
величину dB 1 , |
|
|||||
Пусть широта пункта А изменилась на |
малую |
азимут |
||||||||||||||||
и длина линии АВ - |
на малые величины dA 1 |
2 и ds соответственно. Найдем |
||||||||||||||||
выражения для поправок в координаты пункта· В и в обратный азимут А 2 . 1 , |
||||||||||||||||||
т. е. dB 2 , |
dL 2 |
и dA 2 . |
1 как функции изменений dB 1 , |
dA 1 . 2 и ds. -Учитывая, что |
||||||||||||||
dB 1 , dA 1 . |
2 |
и ds - малые величины, |
имеем: |
|
|
1 |
|
|||||||||||
|
|
|
|
dB |
|
дВ2 |
dB |
|
|
|
дВ2 |
d |
s + |
дВ2 |
dA |
|
|
|
|
|
|
|
2 |
= дВ1 |
|
1 + fis |
|
дА1. |
2 |
1. 2 |
|
|
|||||
|
|
|
|
dL2 = |
:~: dB1 + 8t2 |
ds + аа::.:. |
2 dA 1 • 2 + dL 1 j, |
(33.1) |
||||||||||
|
|
|
|
dA |
= дА2.1 dB + дА2.1 ds+ дА2.1 dA |
|
||||||||||||
|
|
|
|
2, 1 |
|
дВ1 |
|
1 |
|
|
дs |
|
|
дА1. |
2 |
1. 2 |
|
|
или |
|
|
|
|
|
|
|
|
|
+dB: +dBf 1. |
|
|
1 |
|
||||
|
|
|
|
|
dB2 = d.в1; |
1 |
2 |
|
|
|||||||||
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dL2 = dL 1+dLf 1 +dL~ +dLf 1. 2 |
• |
(33.2) |
|||||||||||
|
|
|
|
dA 2. 1 =dAf.\ +dA~. 1 +dA:\ 2 |
|
|
||||||||||||
1. Вы в од вел и чин |
dBf1, dLf 1 |
и dAf\. Пусть точка А', |
лежа |
|||||||||||||||
щая на меридиане точки А (рис. 63), |
имеет широту В |
1 + dB 1 • Будем повора |
чивать геодезическую линию БА вокруг точки В до тех пор, пока она не пройдет через точку А'. Точка А переместится в положение А {. Перемещаем точку А {
(в положение А', сохраняя при этом длину линии АВ, равной s. Тогда точка В
переместится в положение В~ и, очевидно, А'В~ = s. Будем поворачивать ли
вию А'В' вокруг точки А' до тех пор, пока ее азимут не сделается равщ,1м А 1 2 ;
при этом В{ переместится в положение В'. Очевидно, разность широт точек·в
и. В' получится
dBB1 = дВ2 dB |
||
2 |
дВ1 |
1, |
т. е. это будет часть поправки dB 2 , обусловленная изменением широты точки А
ва величину dB 1 .
Из рис. 63 имеем:
вв~ = А1А' = М1 dB1 cos А1, 2,
АА~ =М1 dB1 sinA 1 • 2·
Применяя формулы для решения прямой геодезической задачи, получаем
разность широт точек В{ и В
М1 dB1 cos А1. 2 cos А2• 1
М2
151
а разность широт точек В' и Вi |
равна |
|
|
|
|
М1 dB 1 sin А1. 2 sin А2. 1 |
|
||
|
|
vf2 |
|
|
Так |
как АА; ~ В~В', следовательно, |
|
||
|
.llf1 dB1 cos А1. 2 cos A:i. 1 |
М1 dB 1 sin А1. 2 sin А |
2. 1 |
|
или |
М2 |
|
М2 |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
(33.3) |
Из треугольника АВР (рис. 64), который мы рассматриваем как сфери- |
||||
ческий, |
|
|
|
|
|
cos l = -cos А1. 2cos А2. 1- |
sin А1. 2sin А2. 1 cos а. |
(33.4) |
|
Полагая в выражении (33.4) |
cos а = 1, на основании (33.3) |
и (33.4) полу- |
||
чаем |
|
|
|
|
|
|
|
|
(33.5) |
Для |
вывода dLf 1 заметим, |
что |
|
|
|
L2 = L1 + l; dL2 = dL1 + dl; |
dl = dlBt + dl 5 + dzA1. 2. |
||
По аналогии с предыдущим |
|
|
|
В1 |
M1dB 1 |
cosA1. :isinA 2 |
.1 |
d l2 |
= - ------::------- |
||
|
|
N 2 cos В2 |
|
или
+ M1dB 1 sinA1. cosA .
----=--------2 2 1 '
N2 cos В2
dlf 1 = - |
м1 |
{cos А |
1 |
. 2sin А2. 1 -sin А |
1 |
. |
2cos А2. |
1 |
} dB1. |
|||
|
N2 cos В2 |
|
|
|
|
|
|
|
||||
Но из рис. 64 имеем: |
|
|
|
|
|
|
|
|
|
|
||
sin В2 sin l |
= --sin А2. |
1 cos А 1 . 2 +cos А2. |
1 sin А1• 2 cos а, |
|||||||||
-sin В2 sin l = cos А 1• |
2 ~in А2• |
1 -sin А 1. |
2cos А2. 1 cos а. |
|||||||||
Положим, что cos а = |
1, тогда формула (33.6) |
примет вид |
||||||||||
|
|
|
dlf 1 = |
|
~ 1 |
sin l tg В2dB • |
|
|
|
|
||
|
|
|
|
|
1V 2 |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Для вывода |
dA 2 _1 вспомни:м, что |
|
|
|
|
|
|
|||||
следовательно, |
|
|
А2. 1 = А 1. |
2 ± |
180° + t, |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dA 2 . |
1 = dt. |
|
|
|
|
|
|||
Из треугольника |
сРЬ (рис. |
65) |
|
|
|
|
|
|
|
|||
|
|
|
sinB2 =ctgltgt, |
|
|
|
|
|
||||
|
|
|
tg t = tg l |
sin В2. |
|
|
|
|
|
|||
Дифференцируя (33.9), |
находим |
|
|
|
|
|
|
|
||||
|
~ |
2 |
= __!}_! l sin В2+ tg l cos В |
2 |
dB2 • |
|
|
|||||
|
cos· t |
cos2 |
|
|
|
|
|
|
|
(33 .6)
(33. 7)
(33.8)
(33.9)
152
Полагая cos 2 |
t |
|
= 1 |
и cos 2 |
l |
= 1, получаем |
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
dt = dl sin В2+ tg l cos В2 dB2• |
|
|
|
|
|||||||||||||||
Принимая во внимание (33.8), (33.7) и (33.5), последнее уравнение примет |
||||||||||||||||||||||||||
вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
:: cos l} dB1 , |
||||
|
dA :.1 1 = { ; 1 |
sin l tg В2sin В2+ tg l cos В2 |
||||||||||||||||||||||||
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ М1 |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
B 1 |
|
sinl |
( |
М1 |
, 2в |
2 |
|
|
|
2в}dВ |
|
|
|||||||||
|
|
dA 2.1=--в-{-N Slil |
|
|
|
|
-М cos |
2 |
1 . |
|
|
|||||||||||||||
|
|
|
|
|
|
|
cos 2 |
t |
2 |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
Полагая М1 |
~ М2 |
и имея в виду, что по формуле (5.12) |
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
м2 |
~ 1-e2 cos2 B2 |
' |
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
N2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
_получаем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dА |
в |
1 |
|
sin Z {( |
1 -е |
2 |
cos |
2 В ) |
|
· |
2 В |
+cos |
2 В} dB. |
, |
|||||||||||
|
2 |
· |
1 = - ~ |
|
|
|
2 |
|
sш |
|
2 |
|
2 |
1 |
||||||||||||
|
|
|
|
cos л 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
или |
окончательно |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dA |
B1 |
sinZ {1 |
-е |
2 . |
|
|
2в |
|
|
2B}dB |
1" |
|
(33.10) |
||||||||
|
|
|
|
|
2. |
1 = cosJJ |
2 |
|
S1ll |
|
2 cos |
|
2 |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2. Вы в од |
вел и чин |
dB~, |
dls и dA;.,. Оставляя прежние обозначе |
||||||||||||||||||||||
ния, |
положим, что длина геодезической линии АВ = s |
(рис. |
66) изменилась |
на величину ВВ 1 = ds, причем широта начальной точки и азимут линии АВ
остались без изменения. Таким образом, dB~ - изменение широты точки В,
обусловленное изменением в длине геодезической линии АВ на ds, равно раз
ности широт точек В и В 1 .
..р
с |
ь |
|
А |
|
|
|
|
||
|
Рис. 65 |
Рис. 66 |
|
Рис. 67 |
Азимут линии ВВ 1 |
равен А 2 . 1 - 180°, поэтому по формула,м для решения |
|||
прямой геодезической |
задачи найдем |
|
|
|
|
|
dB~ = -cos А2• 1 (1)2 ds. |
(33.11) |
|
Так как |
ds = s Л~s, то формула (33.11) |
примет вид |
||
|
|
s |
d lg s |
(33.12) |
|
|
dв2 = -cosA 2 • 1 (1)2s - µ - · |
153
Рассуждая та:к же в отношении долготы второй точ:ки и обратного азимута,
получаем |
|
|
|
|
|
|
|
|
|
dLs = _ sin А2. 1 |
( |
). |
ds |
|
|
||||
2 |
|
cos В2 |
|
2 2 |
' |
|
|
||
или |
|
|
|
|
|
|
|
|
|
dLs = _ |
sin А |
2• 1 |
( |
) |
s~ } |
|
|
||
2 |
cos |
В<> |
2 2 |
|
µ. |
|
(33.13) |
||
dA~. 1 = -sin А |
2 |
.:tg В2 |
|
|
|
||||
(2)2 ds |
' |
|
|||||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(33.14) |
3. Вы в од вел и чин |
dBf 1 |
• |
2 , |
dLf 1 • 2 , |
dA{~· 2 • Оставляя |
прежние |
|||
обозначения, предполагаем, что изменился азимут А 1 2 |
геодезической линии АВ |
||||||||
на малую величину dA 1 . 2 , в результате чего точ:ка В переместилась в положение |
|||||||||
В 1 (рис. 67). Определим длину линии, |
соединяющей точ:ки В и В 1 . |
Если бы |
линия АВ располагалась на плоскости, то, очевидно, :кривая малой длины
ВВ 1 , будучи элементарной дугой окружности, равнялась бы длине геодезиче |
||
ской линии s, умноженной на угол dA 1 • 2 , т. е. существовало бы равенство ВВ 1 = |
||
= dA 1 . 2s. |
Но, пос:коль:ку линия АВ расположена на эллипсоиде, то последнее |
|
равенство будет неточным. Напишем |
|
|
|
ВВ1 = du =·mdA 1 • 2, |
(33.15) |
где т - |
функция длины и азимута геодезической линии, при :которой справед |
|
ливо написанное равенство. Величина т называется |
п р и в е д е н н о й д л и - |
|
н о й г е о д е з и ч е с к о й л и н и и. |
|
Для определения приведенной длины геодезической линии, учитывая близость земного эллипсоида :к шару и небольшую величину дифференциальных поправок, примем эллипсоид за сферу с радиусом, равным среднему радиусу
кривизны.
Тогда, рассматривая треугольник АВВ 1 как сферический и выражая его стороны в угловой мере, находим
. du |
. s |
sш~ |
sшR |
|
|
sin dA 1 . 2 - sin 90° |
' |
|
|
или, по малости величин |
dи |
и dA 1 . 2 , |
|
|
|
R |
|
|
|
||
|
|
du = R sin _!_ dA |
1 . |
2· |
(33.16) |
|
|
л |
|
|
|
Сравнивая (33А.5) и |
(33.16), находим |
|
|
|
|
|
|
. • |
|
|
(33.17) |
|
|
m= R sш7r• |
|
||
Изменение широты точки В, вызванное изменением азимута А 1 • 2 |
на dA 1 • 2, |
будет, очевидно, равно разности широт точек В и В1 • |
Имея в виду, что азимут |
|
линии ВВ1 равен А 2• 1 |
+ 270°, получаем |
|
|
dBf1• 2 = т sin А2• 1 (1)2 dA 1• 2 , |
(33.18) |
154
Рассуждая аналогично этомуI находим
|
|
dL:1. 2 = - т cos А2. 1 |
(2) dA |
1. 2· |
(33.19) |
|
|
|
|
cos В2 |
2 |
|
|
Для вывода dAf.\2 заметим, что |
поправка в обратный азимут вследствие |
|||||
изменения прямого азимута должна состоять из двух членов: |
|
|||||
1) из |
поправки |
dA 1 . 2 , отнесенной к |
|
|
|
|
приведенной длине |
геодезической |
линии, |
|
|
|
|
эта часть |
поправки имеет вид |
|
|
|
|
2) |
из поправки, обусловленной изме |
s |
/ |
|
|||
нением сближения меридианов при переме |
Рис. 68 |
|
|
щении |
конечной точки в результате изме |
|
нения начального азимута.
На рис. 68 через В обозначено положение конечной точки линии с азиму
том А1 • 2 ; |
если азимут А 1 . 2 |
изменяется на dA 1 . 2 , то точка В перемещается в В1 • |
|||||||||
Очевидно, |
отрезок |
ВВ1 по-прежнему |
будет |
равен |
mdA 1 • 2 , |
а |
азимут его |
||||
А 2•1±180°+90°. Следовательно, изменение сближения |
меридианов |
в конечной |
|||||||||
точке или сближение меридианов между точками В и В |
1 |
равно |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
dt = dLf1. |
2 sin В2= -тcos А2. 1 (2)2 tg В2dA 1• 2• |
|
|
|||||||
Таким образом, |
полная поправка в обратный азимут будет иметь вид |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
(33.20) |
Но |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
' |
|
|
|
|
|
|
|
|
|
m= R Slll"]ft |
|
|
|
|
|
|
||
поэтому |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dm |
s |
|
|
|
|
|
(33.21) |
|
|
|
|
-- = cos-. |
|
|
|
|
|
|||
|
|
|
ds |
R |
|
|
|
|
|
|
|
Следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
dA:\ 2 = cos ~ |
dA 1 • 2 - R sin · ~ |
cosA2. 1 (2)2 tg B2 dA 1 • 2• |
(33.22) |
|||||||
|
|
|
|
|
в |
|
|
8 |
2 |
|
|
Для упрощения формулы заменим |
cos |
R |
через 1 - |
2R 2 ' |
пренебрегая при |
этом различием между R и N, и положим во втором члене, что R sin; =s;
тогда получим окончательно
(33.23)
f55
Таким образом, на основании формулы (33.1) дифференциальные формулы
первого рода в окончательном виде примут вид:
dв |
2 |
= -1- М1 |
|
, zdB |
1 |
_ |
|
cos А |
2. 1 |
р |
" |
s |
d lg s + т sin А |
2. 1 |
dA |
1. 2 |
|
|
|
||||||||||||||
|
|
1 м |
cos |
|
|
|
|
|
|
м |
|
|
|
µ |
|
|
|
м |
2 |
|
|
|
|
|
|
||||||||
|
|
|
• |
2 |
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
dL 2 = dL1+ М1 |
sin l t |
В dB |
|
- |
|
sin А2. 1s |
р |
" |
d lg s |
- .т cos А2. 1. dA |
1· 2 |
|
|
||||||||||||||||||||
|
|
|
N2 |
|
|
g |
|
2 |
|
|
1 |
|
|
N2 cos |
В2 |
|
µ |
|
|
N2 cos |
В2 |
|
|
|
|||||||||
dА2.1=--в-{1-е |
sш |
2 |
B2cos |
2 |
B2}dB1-sшA . |
|
|
tgB -N р |
---+ |
>• |
(33.24) |
||||||||||||||||||||||
|
|
sin |
Z |
|
|
2 . |
|
|
|
|
|
|
|
|
|
|
|
. |
|
2 |
1 |
|
2 |
|
s |
" d Ig s |
|
||||||
|
cos |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
µ |
|
|
|
||||
|
|
|
|
+ |
{1- |
|
s2 |
|
|
|
|
|
1 |
cosA 2. 1 tg В2 |
} |
dA 1 • 2 |
|
|
|
|
|
|
|
||||||||||
|
|
|
|
2N: |
- s N 2 |
|
|
|
|
|
|
|
|
Эти формулы пригодны для расстояний до 200-250 км.
Точные формулы, справедливые для любых значений s, имеют вид*.
dB2= :: [ cos А1• |
|
2 cos А2. |
1+ ( d:;; |
) 2 sin А1. 2 sin А2. |
1] |
dB1 - |
|||||||||||||||||||||||
|
|
- |
|
|
cos А2 |
1 |
d + т . |
А |
|
. |
dA |
1. 2 |
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
М2. |
|
s |
М2 |
|
s1n |
2. 1 Slll |
|
|
|
|
|
|
|
|||||||||
dL2 = dL 1 - |
N |
М1 |
|
в [cos А1 |
2 sin А2 |
• |
1 - |
( |
ddm) |
sinA 1 |
· |
2cosA 2 |
1 ]dB1 - |
||||||||||||||||
|
2 |
cos |
|
2 |
|
|
|
|
|
|
|
• |
|
|
|
|
S |
2 |
|
|
|
|
|
· |
|
||||
|
|
|
|
|
_ |
|
|
sin А2. |
1 |
ds- |
т cos А2. 1 dA |
2 |
|
|
|
|
|
|
|
, (33.25), |
|||||||||
|
|
|
|
|
|
|
N2cosB 2 |
|
N 2 cosB2 |
l. |
|
|
|
|
|
|
|
||||||||||||
dA |
1 = - |
м |
tg |
А |
|
|
tg В1 |
|
tg В2 |
[ |
cos |
А |
2 cosA 2. 1 |
+ |
|
||||||||||||||
2. |
|
|
1 |
|
|
2. 1 {~+----у;г;- |
|
1 . |
|
||||||||||||||||||||
|
dm ) |
|
|
. |
|
|
А |
|
|
|
. |
]) dB |
|
|
tg R2 . |
|
А |
2. 1 |
d |
s |
+ |
|
|||||||
+(-;fs" |
|
2 sш |
|
|
1. 2 sш А2. 1 J |
1 ---л;:- sш |
|
|
|
|
|||||||||||||||||||
|
|
+ |
|
|
-- |
|
- |
т tg В2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
[ ( |
|
dm ) |
|
|
|
|
] |
dA1. 2 |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
ds |
|
1 |
|
N2 |
|
cosA2.1 |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ргде
dm ) |
т tg В1 |
tg А1. 2 |
( ds 2 |
= N1 cos А1. 2 |
+ tg А2. 1 • |
Для топографических и картографических целей,
когда расстояния не превышают 40-50 км, а поправкю
Акоординат достаточно иметь с точностью 0,001-0,002",
формулы (33.24) можно упростить. Пренебрежем сфе
Рис. 69
1. Принимая Mi
М2
роидичностью Земли и заменим приведенную длину гео
дезической линии т длиной геодезической линии s (:&
данном случае дугой большого круга), тогда:
= 1, согласно формуле (33.5), получаем
dBf 1 = cos l dBi- |
(33.26) |
* См. [44, стр. 278].
156
2. Аналогично предыдущему, учитывая формулу (.33.18), получаем
dBf1. 2 = + ; sin А2. 1 dA 1 . 2· |
(33.27) |
Построим сферический треугольник АВР (рис. 69) и на меридиане точки А выберем точку С, как это делалось при выводе формул для решения прямой геодезической задачи.
Обозначим l = а. Из треугольника АВС имеем
с
а= S.Ш А1. 2 '
из треугольника СВР получим sin с = eos В2 sin Z или с = cos В 2 sin Z. Имея
ввиду последние выражения, можем написать:
|
dBf1. 2 = а sin А2 |
• |
1 dA |
1 |
2 = |
. : |
|
|
siп А2. 1 dA 1. 2, |
||||
|
|
|
|
|
' |
Slll |
1.2 |
|
|
|
|||
|
|
dBf1. 2 = -sin Zcos В2dA 1 . 2 • |
|
(33.28), |
|||||||||
3. |
Полагая в формуле (33. 7) М1 |
= 1, получаем |
|
|
|||||||||
|
|
|
|
N2 |
|
|
|
|
|
|
|
|
|
|
|
dL: 1 = sin Ztg B2 dB1 • |
|
|
(33.29) |
||||||||
4. |
Согласно |
формуле (33.19), |
полагая, |
что |
т |
|
с |
||||||
-а---- и с= |
|||||||||||||
sin l cos В 2 , получаем |
|
|
|
|
|
|
|
|
N 2 - |
- |
sin А1.2 |
||
|
|
|
.,=- |
|
|
|
|
|
.,= |
||||
|
dLA1,2=-<JcosA2.1aA |
С |
|
cosA?.1dA |
|||||||||
|
2 |
cos В2 |
|
l,w |
siп А1. 2 |
ccs В2 |
|
I,w |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= - |
|
sin l cos В2 |
cns А 2 |
.1 |
|
|
|
||||
|
|
|
sin А 1.2 |
cos В2 |
|
dA1.2· |
|
|
|||||
|
|
|
|
|
|
|
|
|
|||||
Окончательно находим |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
A1,2 _ |
|
, |
[ C(JS А2.1 dA |
|
(33.30), |
||||||
|
|
dL 2 |
|
- |
- |
Slll |
. |
А |
1. 2 |
1,2, |
|
||
|
|
|
|
|
|
|
Slll |
|
|
|
|
||
5. |
Считая эксцентриситет е равным О и учитывая формулу (33.10) получаем, |
||||||||||||
|
|
|
dAB 1 |
|
=~dB |
1· |
|
|
(33.31), |
||||
|
|
|
|
2.1 |
cosB2 |
|
|
|
|
|
|||
6. |
Согласно |
формуле (33.22) |
|
|
|
|
|
|
|
|
|
||
|
|
dAf~ 2 = ( cos |
|
~ - |
R sin ~ cos А2•1 |
(2)2 tg В2) dA1,2, |
Придерживаясь принятых обозначений и по-прежнему пренебрегая сфе
роидичностью Земли, получаем
йли
dA Ан= ( |
cos а cos B 2-sin а cos А2.1 sin В2 |
}' dA . |
2.1 |
cos В2 |
1,2 |
Из треугольника АВР (см. рис. 69) |
|
|
cof: В1 cos l |
= cos а cos В2- si.n а cos А2,1 sin В2, |
157
поэтому |
|
|
|
|
dA А1•2 _ |
cos В1 cos l |
dA |
|
2.1 - |
cos В2 |
1.2, |
но |
|
|
|
cos В1 |
sin (90° -В1) |
sin А2.1 |
|
cos В2 |
= sin (90° -В2) - |
sin А 1• 2 • |
~ледовательно, окончательно
(33.32)
Выражения ав;, dL1 и dA~. 1 остаются такими же, как и в формулах (33.24).
Сделав все эти преобразования, получим дифференциальные формулы первого
рода в упрощенном виде |
|
|
|
|
|
|
|
|
|
|
|
dB2= cos l dB1 - |
|
. |
|
|
|
(1)2cos |
A2.1s |
л lg s |
) |
||
cos В1 sш l dA 1•2- |
-µ- |
|
|||||||||
dL2 = dL1+ si11 l tgB2dB |
1 |
-sinl c?s~ |
2· 1 dA . - |
|
|||||||
|
|
|
|
|
|
s1n |
1,2 |
1 |
2 |
(33.33) |
|
|
|
|
|
|
|
|
|
|
|||
|
- ( |
) |
sinA 2.1 |
Лlgs |
|
|
|
||||
|
|
|
|
|
|||||||
|
|
2 2 |
cos В2 8 |
µ |
|
|
|
|
|||
sin l |
sinA 21 dA |
( |
. |
А |
|
В |
Л lgs |
||||
dA2.1 = --в- dB1 - |
cos l |
. А |
. |
|
1.2 - |
|
2)2 sш |
2.1 tg |
2S -- |
||
cos 2 |
Slll |
1.2 |
|
|
|
|
|
|
|
µ |
Для нелогарифмического вычисления поправок формулы (33.33) перепи
шутся: |
|
|
|
|
|
|
|
dB2 = cos zdB1 - cos В1 sin l dАы- 0,03234" cos А2.1 ds 1 |
|||||||
dL2 = dL 1 +sin l |
tg В2 dB1- |
sin l |
~osAA 2 · 1 dA 1•2 - |
||||
|
|
|
|
|
s1n |
1 . 2 |
|
|
- |
О 03234" sin А2· |
1 ds |
(33.34) |
|||
|
|
' |
|
cos В2 |
|
|
|
_ |
sin l |
dB |
cos |
z sin А2. 1 |
dA |
||
dA 2. 1 - |
cos В |
i - |
sin |
А |
2 |
1. 2- |
|
|
2 |
|
|
|
1. |
|
-О,03234" sin А2• 1 tg В2ds
·где принято с подстановкой приближенных числовых значений:
1 |
s |
== (2)2 s d |
1 |
= О,03234" ds. |
|
(1)2 s d : |
: s |
|
|||
Из сравнения упрощенных |
формул |
с |
более точными |
формулами (33.25) |
|
,следует, что в первых ошибочны коэффициенты при dA 1 |
и dB 1 на величины |
||||
порядка e2cr. Следовательно, |
при расстояниях в 40-50 км и значениях dA 1 |
и dB 1 в несколько секунд погрешности, обусловленные приближенностью фор
мул, могут достигать величины порядка О,001~'. Формулами (33.33) следует
пользоваться лишь при вычислении координат пунктов как опорных для съемок.
R дифференциальным формулам I рода следует отнести формулы, служащие для решения обратной, по сравнению с рассмотренной, задачи: определить
·i.58
изменения длины геодезической линии s, прямого и обратного азимутов ее
'.А 1 • 2 и А 2 . 1 , вызванные изменением координат конечных точек данной линииt
т. е. В 1, L 1, В 2, L 2 на dB 1, dL 1, dB 2 и dL 2.
Искомые формулы получаются путем алгебраических преобразований полученных выше дифференциальных формул.
Опуская вывод, приведем вывод формул в окончательном виде*:
p"ds= -M2 cosA 2• 1 dB2 - r2 sinA 2• 1 dL2 |
1 |
|
т dA 1. 2 = М2 sin А2• 1 dB2- r 2 cos А2• 1 dL 2 |
(33.35) |
|
тdA,. 1 =М,( ~7 ),sinA,. 1 dB2+r1 cosA 1• ,dL, |
||
' |
Написанные формулы соответствуют случаю, когда изменились координаты
второй точки, т. е. В 2 , L 2 • Если изменить индексы <<1» на <<2>>, то формулы будут
справедливы для случая, когда изменяются координаты первой точки, т. е.
В1, L1.
При одновременном изменении координат точек 1 и 2 формулы примут вид
р" ds = -М1 cos А 1. 2 dB1- М2 cos А2• 1 dB2 - r2 sin А2• 1 (dL 2 -dL 1 ) )
mdA 1• 2 =M1 |
(dd: ) sinA 1• 2 dB1 +M2 sinA 2• 1 dB2 - |
|
|
2 |
|
|
- Г2 cos А2. 1 (dL2-dL1) |
(33 36} |
т dA 2• 1 =М2 |
( ~7) 1 sin А2. 1dB2+M1sin А1 |
• 2 dB1 + |
-л- |
+r1 cosA 1 .2(dL2-dL1) |
|
f; • |
|
|
r;це r - радиус параллели под данной широтой.
;,1;,
!(:
§ 34. Дифференциальные формулы второго рода
J,:j
•• ,., Пусть некоторая триангуляция вычислена на поверхности эллипсоида,
1
~~з:меры которого определяются большой полуосью а1 и сжатием а.1 (или экс центриситетом е1). Rак известно из гл. IV, геодезические координаты пунктов lРИангуляции вычисляют путем последовательной передачи разностей коор
~ат смежных пунктов. При вычислении разностей координат используют
8~11,?Вные геодезические величины (1), (2), которые являются функциями боль
р,хои полуоси, эксцентриситета эллипсоида и широты. Если возникает задача
~еревычисления триангуляции u на поверхность нового эллипсоида, размеры
~~торого определяются большои полуосью а2 и сжатием а.2 (или е2), то, оче JИдио, изменение размеров эллипсоида вызывает изменение разностей коор
-риат. Отсюда следует, что вывод формул для поправок за изменение парамет
ров эллипсоида должен заключаться в нахождении формул для поправок в раз
~-с()tти координат пунктов ходовой линии, по которой вычислялись координаты.
~ормулы этих поправок легко найдутся путем дифференцирования главных
Членов известных выражений для разностей широт, долгот и азимутов.
• См. [2, стр. 235, 236].
159'