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

Закатов Вища геодезія 1

.pdf
Скачиваний:
9
Добавлен:
28.06.2022
Размер:
23.3 Mб
Скачать

Щ

Исходнне условия

для получения

зтих функций: а) конформ-

ность проекции, б) внбор зон и осей

координат

описанньїм

вьіше

способом, в) изображение

осевого меридиана

и зкватора

на плоскости

пря­

мими линиями, г) масштаб изображения по осевому меридиану равен единице. Пусть на рис. 76: ОЕ — зкватор зллипсоида, а РО — осевой меридиан.

Дана точка А, имеющия геодезические координати — широту В и долготу L; точка А х расположена на бесконечно малом расстоянии от точки А в произволь-

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

от нее;

координатами точки А х будут В +

dB и L -J- dl.

Через X назовем

длину

дуги осевого меридиана от зкватора

до параллели

с широтой В.

Пусть прямне ОхжОу — ось абсцисс и ось ординат соответственно (рис. 77), расположеннне перпендикулярно одна к другой, суть изображения зкватора

Р

Злпипсоид

ОсеВой

Рис. 78

и осевого меридиана на плоскости. Пусть точка а — изображение на плоскости точки А; точка ах, находящаяся на бесконечно малом расстоянии от точки а, соответствует точке А х. Обозначим прямоугольнне координати точки а на пло­

скости через х и у, а точки ах— через х +

dx и у +

dy. Обозначим расстояние

A A j (рис.

76) на зллипсоиде через ds, а расстояние аахна плоскости — через dS.

Построим

вспомогательннй треугольник

А А ХС и

будем рассматривать его

как злементарннй, тогда найдем, что А ХС равно злементу дуги меридиана MdB,

АС — злементу дуги параллели

N cos Bdl.

 

Отсюда

 

 

 

ds=*V{M dBf +

(N cos В d l f ,

 

ИЛИ

 

 

 

& =

 

 

(37.2)

Введем обозначение

dB

,

 

M

 

N

cos В ~

dq'

 

тогда

 

 

 

ds - N cos В V (іi q f + (dlf.

(37.3)

Для плоскости

 

 

 

d£= V ( d x f + (<i y f ..

(37.4)

Масштаб изображения внразится так:

 

170

или

m

V(dx)24-(dy)2

 

(37.5)

 

N cos В V(dq)2-f- (dl)2 '

 

 

 

 

Перепишем уравнение

(37.5) в виде

 

 

 

(dx-\-i dy) (dx і dy)

 

(37.6)

 

N2 cos2 В (dq-\- і dl) (dq і dl)

где

 

 

 

 

i = V = T -

Уравнение (37.6) вьіражает масштаб изображения для произвольного закона изображения. Для обеспечения конформности изображения масштаб т по всем направленням должен бить одинаков и зависеть только от координат данной точки, но не от направлення злемента ds. Следовательно, вид функции, определяющий зависимость между координатами на плоскости и геодезическими координатами на зллипсоиде, должен бить таков, чтобьі внражение

(37.6) для масштаба не зависело от отношения дифференциалов ^или-^у-,

которне определяют направление отрезков dS или ds.

В теории функции комплексного переменного доказивается, что Р + iQ только тогда является аналитической функцией комплексного переменного

„ І

І

P +

1Чі

село от

когда

d ( P + i Q )

dp

наоборот, чтобьі

d(P + iQ)

^

' не зависит от

dq

, , '

. . -не зави-

 

d (Р +

*?)!

^

d

гд)

необходимо, чтоби P-f- iQ бьшо аналитической функцией комплексного

переменного р + iq. Имея

зто в виду, переписиваем уравнение

(37.6) так:

2 __

1

d(x + iy)d(x — iy)

І

/о 7 7 Ч

 

N*cos*B

d (q + il)d {q— il)

*|

\ОІ-П

Согласно изложенному вшше свойству функций комплексного переменного, для того чтоби масштаб изображения не зависел от направлення отрезка ds

или от ^|, необходимо, чтоби х + іу било некоторой аналитической функцией

от q + il или х іу — аналитической функцией от (q il), т. е. чтоби

х-\- іу = /(<? + U).

(37.8)

Справедливость и достаточность последнего условия можно еще доказать следующим образом. Если имеет место условие (37.8), то

d (ж + іу)

d (q-\~il) = Ґ (я+іі)-

Заменив в (37.8) + і на —і, получим:

х—іу— f(qа)

и

d (х іу)

f ( q - i l ) .

d (qil)

(37.9)

(37.10)

(37.11)

Следовательно, внражение (37.11) имеет место, если поставлено условие (37.8). На оснований формул (37.9) и (37.11) внражение (37.7) примет вид

щ2 =

=

+

(37.12)

171

Производньїе /' (q -{- il) и /' (q il) при наличии условий (37.8) и (37.10) зависят только от координат х, у или q, І, но не зависят от их дифференциалов, позтому не зависят от них и вьіражения (37.12) и (37.7), что и требовалось до­ казать.

Условие (37.8) обеспечит конформность изображения при любом произвольном виде аналитической функции /.

Для получения конформного изображения зллипсоида на плоскости по Гауссу вид функции / определяют путем введення следующих дополнительньїх

условий:

 

 

 

1)

осевой меридиан ОР (рис. 78) зллипсоида изображается на плоскости

прямой, являющейся осью

абсцисс; следовательно, в уравнении (37.8) при

І =

0

ординатьі у должньї

бьіть равньї нулю;

 

 

2)

для точек осевого меридиана абсциссьі х должньї бшть равньї соответ-

ствующим дугам X, т. е. дугам меридиана от зкватора до данной точки с широ-

той

В.

(37.8) при у = 0

 

 

Отсюда в уравнении

 

 

 

 

x = f { q ) = X .

(37.13)

Разложим / (q + il) в ряд Тейлора, предполагая, что величина іі сравнительно невелика. Будем иметь

f { q + i l) =

f {q )+ il df (?)

 

(ііу

d2f (q)

 

(il)

3 d 4 ( q )

 

ИЛИ

 

dq

 

 

dq2

 

 

dq%

 

 

 

 

 

 

 

 

 

 

 

x -j- iy X -j- il

df (q)

___ 1_ /2

d2f (q)___ 1_ 73

dZf (q)

+

dq

2

dq2

6

 

dqs

і

1

l i dbfjq)

1

t75 d4 (q ) ___ 1

,e d4 (g)

 

^

24

dq*

•“

120

dqb

720

 

dqe

 

Производньїе в зтом ряде должньї бшть внчисленьї при І = 0, чего они обращаются в производньїе

dX .

d*X .

 

.

 

 

 

~dq~ ’

dq2 ’

dq*

 

 

 

позтому

d X

12

d^X

i P

d *X

l*

d*X

x —j—iy X —|—il

dq

2

dq2

(J

dqs

"T 24

dq* T

 

il*

dbX

 

d *X .

 

 

 

120

dq*

720

dqG

+

' • • •

 

(37.14)

вследствие

(37.15)

Приравнивая порознь действительньїе и мнимьіе части в полученном вира­ жений, находим в общем виде основньїе уравнения, о п р е д е л я ю щ и е з а к о н и з о б р а ж е н и я т о ч е к з л л и п с о и д а на п л о с к о с т и в проекции Гаусса:

12

d*X

.

I*

d*X

d е х

(37.16)

2

dq*

 

24

dq*

720

dq*

 

 

dX

 

tPX

Z&

d*>X _

(37.17)

dq

 

6

dqs

*

120

dqb

 

 

 

 

Очевидно, чем больше Z, тем больше должно бшть удержано членов вуравнениях (37.16) и (37.17), т. е. чем больше будут искажения на краях зоньї, тем учет зтих искажений сложнее.

Как указьівалось вьіше, в настоящее время наибольшая протяженность зон по долготе принята 6°. Таким образом, удаление пунктов от осевого меридиана по долготе не превосходит 3°, позтому последующие рабочие формули для вичислений при применении проекции Гаусса — Крюгера и рассчитанн

на интервал в долготе І

3°.

§ 38. Формульї для определения конформних плоских координат х н у

по геодезическим координатам В и X

Вьівод рабочих формул для внчисления координат Гаусса — Крюгера по геодезиче ким координатам, очевидно, сводится к нахождению производннх

dnH и подстановке найденннх их

значений в

уравнения

(37.16)

и (37.17).

Введем обозначения:

 

 

 

 

t g B — t\

ц2 — е'г cos3В;

1 т)2 =

V 2.

(38.1)

Так как dX — дифференциал дуги меридиана, то

dX — М dB,

или, согласно (5.9),

 

 

dX = - 4

dB.

(38.2)

При виводе основних формул проекции ми обозначили

 

dq — М

dB

 

N

cos В 9

 

или, учитнвая формули (5.9) и (5.10), можно написать

dq = — - - „ .

(38.3)

^

V2 cos В

v

'

Получаем внражение для первой

производной

 

 

f

=

f c o s B .

(38.4)

Применяя правило дифференцирования сложной функции, пишем в общем виде внражение для второй производной

 

d*X

d*X

dB

 

(38.5)

 

dq2

dq dB

dq *

 

 

 

 

Дифференцируем (38.4) по В , имея в виду, что V єсть функция J5,

 

d?X

с dV

D

c

. D

(38.6)

dq dB

v * - d B cosB-

~

s m B -

 

Для нахождения dV

дифференцируем

формулу V2 = 1 + e'2cos2 В:

dB

2V d V = — 2e'2 cos В sin В dB, V dV — — e*2 cos2 Big В dB,

173

или

 

 

dV

_

ч*

 

 

(38.7)

 

 

dB

~

V

 

 

 

 

 

 

 

Подставляя последнее вьіражение в (38.6), получаем

 

d*X

=

-уз Т]2* cos В —-^-sin В =

уз sin В.

 

dq dB

 

Так как из формули (38.3) dB =

 

V2 cos В, то на

оснований формули

(38.5)

получаем вторую производную

 

 

 

 

 

 

 

d *X

 

с

sin 2В

 

(38.8)

 

 

dq1*

~2

V

 

 

 

 

 

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

=

N, значення первой и второй производних напишем

так:

 

d X

 

 

 

 

 

 

 

 

N cos В

 

 

 

 

dq

 

 

(38.9)

 

 

 

 

[•

 

 

d*X

 

N

 

 

dq2

 

2

sm 25

 

 

Пользуясь формулами (37.16), (37.17) и (38.9), находим главньїе члени виражений для х н у как функций геодезических координат:

х — Х-\ — 1"%sin В cos В + ..

2р"2

(38.10)

у = -^-N cos В . ..

 

Вираження (38.10) для х н у , полученние без учета второго и последующих членов рядов (37.16) и (37.17), для точних вичислений непригоднн. Для получения формул, позволяющих (вести внчисления требуемой в триангуляции

точности, необходимо вичислить производнне] до шестого порядка включи-

тельно и после подстановки их значений в формули (37.16) и (37.17) получить вираження для х н у .

Не приводя викладок, связанннх с получением указанннх производних внсшего порядка, и последующих преобразований при подстановке зтих про-

изводннх в (37.16) и (37.17), напишем окончательний

результат

X = X + X X N sin В cos В | l

+

^

0 .s. (5_,2^9>і2 +

4т,‘) +

+ ' 1

^

(6 1 -5 8 (2 + (4 )} ’

 

(38.11)

у = X N cos В {і +

 

(1 _

А +

п*) +

 

+ ‘" 1 ° ° ^ - (5 -

18(2 +

і*+ 14л2 -

58г]2(2) | .

(38.12)

174

§ 39. Формульї для вьічисления геодезических координат

по координатам Гаусса — Крюгера

Вьівод настоящих формул дадим по методу, предложенному Г. В. Багратуни *.

Решая уравнения в общем виде относительно q -f- іі и q іі, получаем

q-\~ИF (х+ іу)

(39.1)

q — и —F (х— іу)

Применяя ряд Тейлора к (39.1), находим

Ч ІІ- F (х+іу) ^F(x) + iyF1 (х) —j - F n (х) —

- i f F m (x) + ^ F l V {x) + . . .

 

 

q— i l = F ( x iy)=F (z)— iyF1 (х)— ~

F11 (х) +

 

 

 

 

+ i f - ^ m

(z) + - i ^ IV( x ) + ...

 

(39.2)

 

Беря полусумму и полуразность виражений (39.2), получаем

 

 

 

 

q = F (х) -

Щ- F11 (х) + -g- ^Іт (х) — .J

.

(39.3)

 

 

 

 

l ~ y F l ( x ) - £ - F m {x)+....

 

(39.4)

 

Вираження

(39.3) и (39.4) внтекают из условия конформности проекции

с зллипсоида на

плоскость.

 

 

 

 

 

 

ции

Для вьічисления рядов (39.3) и (39.4) по-прежнему вводим условия проек­

Гаусса — Крюгера

плоскости на

зллипсоид:

х

Ппоскость

 

 

1) при

у = 0 должно бьіть

1

— 0 и

2) F (х) =

 

 

 

 

= 0і-

79 заданная

абсцисса

х точки а опреде-

 

У

 

 

На рис.

 

 

 

ляется прямой Оех, которая, согласно данному ус-

 

^3o6pafe0 U

 

ловию, должна равняться длине дуги меридиана от

 

 

 

паролі

 

зкватора до

некоторой точки Е х (рис. 78), широту

х<

 

 

которой обозначим через

В г;

следовательно, при

 

 

У =

0 и І =

0 вьіражение (39.1)

должно

обратиться

 

 

 

 

 

в равенство

 

F(x) = qі ,

 

 

(39.5)

 

Изображение зкдатора

 

 

 

 

 

 

где

qx вьгчисляется как

функция широтьі В х\ ши­

 

 

 

рота В х может бьіть сразу легко получена по х из

 

Рис. 79

 

таблиц длин дуг меридианов.

 

 

 

 

 

 

 

Таким образом, формула (39.5) внражает рассматриваемое второе условие

проекции Гаусса — Крюгера.

 

 

 

 

 

 

* Труди МИИГАиК. Вші. 24, стр. 57—59.

175

С учетом поставленньїх условий уравнения (39.3) и (39.4) примут вид:

д

У%( d24 \

i V L ( di(i \

(39.6)

 

2 l < t e 2 j i + 2 4 l W i '

 

 

/ =

«

( J l )

_ l i

(d^L)

! *

(39.7)

 

у

\ d x ) і

6

V. dxs /

 

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

M d B

dq

A COS JO

M d B

УJ N cos В *

На оснований последнего в общем виде можем написать

B = y(q) = y[q1+ ( q - q 1)\,

(39.8)

в і = ц Ш -

Применяя ряд Тейлора и заменяя (q qx), на оснований (39.6) полупаєм

D _D

Г У2 ( а ч \

у і / diq \ -1

СІВ

(39.9)

1_—

24 V d ^ / i J

dq

 

(при зтом учтено, что всегда

В х > Б ) .

 

 

 

Формульї (39.4) и (39.9) в общем виде решают задачу перехода от прямо-

угольньїх координат Гаусса — Крюгера

к геодезическим.

 

Для получения рабочих формул находим необходимьіе производньїе.

Имеем

 

 

 

7

M d B

 

 

 

 

 

 

 

 

 

(39.10;

 

 

 

 

4 ~

N cos В '

 

 

 

 

 

 

 

Имея в виду,

что

dx =

M dB ,

полупаєм

 

 

 

/

dq \ __

 

1

 

/

dq \

1

(39.11)

 

\

dx )

N cos и*

\

dx ) і

Лті cos Вх

или

 

 

 

 

 

dq

 

 

 

 

 

 

 

 

 

 

т»

 

(39.12)

 

 

 

 

 

dx

 

 

из (39.10)

 

 

 

 

 

 

 

 

 

 

 

 

N\ cos Bi

 

 

 

 

 

( — )

 

 

(39.13)

 

 

 

 

 

M~x

 

 

 

 

\

dq J і

 

 

 

Переходим к

вьічислению

вторьіх

производньїх

 

dlq

_ d

/ dq \

dB

d2q

dB

dx2

dB

\ dx )

dx

dB dx

dx

d

^ dcl } _

d

1_

1

dr

dB

dx )

dB

r

r'l

dB

Внчисляем ^ dB

(39.14)

(39.15)

dr

d

(^N cosB) = N sin5 + cos В ~ - \ a { 1 — e2sin2#)"'/*] =

~dB

dB

 

 

 

= —N sin B-\- cos В ae2 sin В cos В

 

 

(I —*>=sin? B)

176

откуда

после простьіх преобразований

долучаєм

 

 

 

-^L =

—м sin В.

(39.16)

 

 

 

 

аВ

 

 

 

Следовательно,

 

 

 

 

 

 

 

 

d2q

 

1 П/Г .

D 1

sin В

 

 

* г = 7і Л

 

іі£ ї

 

или

( & д \

sin Ві

 

sin Ві

 

 

 

 

 

\ dx2 / і

 

 

 

N2 cos2 Вх

N \ cos L 1

Принимая во внимание (39.11), (39.13) и (39.16), искомьіе вьіражения (39.7),.

(39.9)

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

и І приближенно напишем так:

 

 

 

І

 

у

 

 

 

 

 

Ni cos Вх

 

 

 

 

 

 

 

 

В = ВХ

У2Н

(39.17)

 

 

2NXM X

 

 

 

 

 

 

Вьічисляя в вираженнях (39.7) и (39.9) производньїе вмсшего порядка, долучаєм формули, удовлетворяющие по точности все случаи практики.

Опустив внчисления зтих производннх и последующие преобразования, напишем точние формули в окончательном виде:

Г

N і cos В і Р

6N 0 + 2 < ї+ л ї) +

 

 

 

 

— (5 -f 2&t\ -j- 2kt\ -f- 6т]ї -f- 8т]і^)

(39.18)

 

120Л'^

 

 

В = ВХ

У~г

(5+ 3^ + —9т)^ї) +

 

 

2МХМ

 

 

 

360N*

(61 + 90^ + 45^4)

(39.19)

 

 

 

 

§ 40. Формули для вьічисления

 

 

сближения

меридианов на плоскости

 

Сближение меридианов может бить виражено в функции геодезических координат и в функции прямоугольннх координат.

1. С б л и ж е н и е м р и д и а н о в на п л о с к о с т и в ф у н к ­ ц и и г е о д е з и ч е с к и х к о о р д и н а т . Пусть точка а (рис. 80) — изображение на плоскости точки А зллипсоида; ns — изображение на плоско­ сти меридиана, проходящего через точку A; aw — изображение на плоскости параллели, проходящей через ту же точку А ; ab — прямая на плоскости, параллєльная оси абсцисо; ае — прямая, параллельная оси ординат; у — сбли­

жение меридианов на плоскости,

равное углу при точке а между ирямой аЬ

и кривой an.

 

 

 

 

Очевидно, угол равен

также

углу

между изображением

параллели aw

и ирямой ае. Возьмем на

изображении

параллели aw точку

а х, бесконечно

12 п. С. Закатов

 

 

 

1г/7

близкую к точне а. Разности координат точек а и аг будут равнн dx и dy. Из еле­ ментарного треугольника ааха 2 (рис. 81) имеем

tgy =

dx

(40.1)

~dy

dx

причем щ находится из уравнения кривои aw, для которои, как для дуги па-

р ал л ел и , широта В должна бьіть принята постоянной величиной (можно взять

злементарннй треугольник, используя злементьі кривой an и прямой аЬ\ однако зто невьігодно, так как в зтом случае широту В нельзя считать по­

стоянной величиной, что усложнило бм последую щ ее дифференцирование).

dy а

У

Уравнение (40.1) можно переписать так:

 

dx

 

tgy

ЧГ

(40.:

dy

 

dl

 

Производньїе ^

и

найдем путем дифференцирования вьіражений (38.11)

и (38.12). Из (38.11) с ошибками на величини второго порядна имеем

 

 

 

dx

lnN sin В cos В

 

 

 

чи

~я2

 

и из (38.12)

 

 

Р

 

 

dy

N cos В

 

 

 

 

 

 

dl

 

 

На оснований

(40.2)

с ошибкой на величину третьего порядки

 

 

 

 

у

(40.3)

 

 

tg 7 = -r r sin В.

 

 

 

г

 

В результате дифференцирования формул (38.11) и (38.12) по І с учетом всех членов и после подстановки в (40.2) получим

 

І" sin В

 

//З

 

tg?

\-------о- sin В cos2 В (1 -f-і2+ Зг)2

2ц4) -J-

 

Р"

*

Зр'

 

 

Г'5

sin В cos4 В (2 + At24- 2г4).

(40.4)

 

15р"

 

 

178

Так как

 

 

 

 

 

у = t g у

f y

+ -^tgby,

 

то из (40.4) имеем окончательно

 

 

 

у" = Г sin В +

sin В cosі2 В (1 +

Зт)2+

2г)4) +

sin В cos4 В (2 t2). (40.5)

2. С б л и ж е н и е м е р и д и а н о в н а п л о с к о с т и в ф у н к ~

ц и и п л о с к и х

к о о р д и н а т .

Внражение для сближения меридианов

в функции плоских координат получится путем заменьї в (40.5) І через прямоугольнне координатьі и заменьї аргумента В через В х. Последняя замена осуществляется путем предварительного разложения по строке Тейлора виражений:

sin В

sin {Вх — (Вх— В)} = sin Вх— (ВхВ) cos Вх

sin Вх (ВхВ)2

 

 

 

 

sin В =

sin Вх----- sin Вх(і + ї]і ) -f- ■5уі sin Вх

 

 

 

 

2N\

14

'

 

(40.6)

 

 

cos2 В = cos2 Вх^1 +

t\ j

 

 

 

 

 

После указанннх подстановок получим окончательно

 

У

- к п

 

 

 

 

 

(40.7):

 

 

 

 

 

 

или в логарифмическом виде

 

 

 

 

 

 

 

W2

(і — г]і cos2 Вх — 2r]J C O S2 Вх)-\-

 

 

 

 

 

 

 

 

і

(13— 6 cos2 Вх).

 

(40.7')

 

 

~

90iV{ cos* Вх

 

 

 

 

 

 

 

Нетрудно видеть,

что главнне члени для у соответствуют разности прямого

и обратного азимутов в главной геодезической задаче, если у заменить

через

і sin А 1Л.

 

 

 

 

 

 

 

§41. Формули для вьічисления масштаба изображения

Для получения формул масштаба изображения в функции геодезических координат согласно формулам (37.2) и (37.4) напишем

 

dx2 4- dt/2

 

 

/ dx \ 2

 

 

т

+

\~dt/ )

(41.1)

* ■ - ( £ ) ■

M2dB2-fjV2 cos2 В dl*

JV2 cos2■2в [ і

 

\ 2\

 

 

 

\

N cos В dl J

)

На оснований (40.1) будем иметь

 

 

 

 

dx V =

1 -j- tg2 у = sec2 у.

 

 

 

* + ( f

 

 

 

 

12*

179

Соседние файлы в предмете Вища геодезія