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

Литература / Закатов П.С. - Курс высшей геодезии (1976)

.pdf
Скачиваний:
732
Добавлен:
10.06.2017
Размер:
34.58 Mб
Скачать

Глава 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 А21

М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 А12 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в}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

 

 

2

 

2

1

 

 

 

 

cos л 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dA

B1

sinZ {1

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 А21 (1)2 ds.

(33.11)

Так как

ds = s Л~s, то формула (33.11)

примет вид

 

 

s

d lg s

(33.12)

 

 

2 = -cosA 2 1 (1)2s - µ - ·

153

Рассуждая та:к же в отношении долготы второй точ:ки и обратного азимута,

получаем

 

 

 

 

 

 

 

 

 

dLs = _ sin А2. 1

(

).

ds

 

 

2

 

cos В2

 

2 2

'

 

 

или

 

 

 

 

 

 

 

 

 

dLs = _

sin А

21

(

)

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°

'

 

или, по малости величин

и dA 1 . 2 ,

 

 

 

R

 

 

 

 

 

du = R sin _!_ dA

1 .

(33.16)

 

 

л

 

 

Сравнивая (33А.5) и

(33.16), находим

 

 

 

 

 

. •

 

 

(33.17)

 

 

m= R sш7r•

 

Изменение широты точки В, вызванное изменением азимута А 1 2

на dA 1 • 2,

будет, очевидно, равно разности широт точек В и В1

Имея в виду, что азимут

линии ВВ1 равен А 21

+ 270°, получаем

 

 

dBf12 = т sin А21 (1)2 dA 12 ,

(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 ,

а

азимут его

А 21±180°+90°. Следовательно, изменение сближения

меридианов

в конечной

точке или сближение меридианов между точками В и В

1

равно

 

 

 

 

 

 

 

 

 

 

 

 

 

dt = dLf1.

2 sin В2= cos А2. 1 (2)2 tg В2dA 12•

 

 

Таким образом,

полная поправка в обратный азимут будет иметь вид

 

 

 

 

 

 

 

 

 

 

 

(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) дифференциальные формулы

первого рода в окончательном виде примут вид:

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-е

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

 

 

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 .

(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

 

 

(33.31),

 

 

 

 

2.1

cosB2

 

 

 

 

 

6.

Согласно

формуле (33.22)

 

 

 

 

 

 

 

 

 

 

 

dAf~ 2 = ( cos

 

~ -

R sin ~ cos А21

(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 А 12

~ледовательно, окончательно

(33.32)

Выражения ав;, dL1 и dA~. 1 остаются такими же, как и в формулах (33.24).

Сделав все эти преобразования, получим дифференциальные формулы первого

рода в упрощенном виде

 

 

 

 

 

 

 

 

 

 

 

dB2= cos l dB1 -

 

.

 

 

 

(1)2cos

A2.1s

л lg s

)

cos В1 l dA 12-

-µ-

 

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

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 12 -

 

 

 

 

 

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 А21 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 А21 dB2- r 2 cos А21 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 А21 dB2 - r2 sin А21 (dL 2 -dL 1 ) )

mdA 12 =M1

(dd: ) sinA 12 dB1 +M2 sinA 21 dB2 -

 

2

 

 

- Г2 cos А2. 1 (dL2-dL1)

(33 36}

т dA 21 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'