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

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

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

осей будем иметь различнше системьі координат, которьіе, оставаясь сфероидическими, будут иметь свои особенности. Система сфероидических координат, являясь системой криволинейньїх координат на поверхности зллипсоида, родственна системам, указашшм в п. З и 4.

Дадим описание наиболее простой системьі прямоугольньїх сфероидических

координат.

І

Примем некоторую точку А

(рис. 7), геодезические координати которой

известньї, за начало координат. Меридиан, проходящий через точку А, примем за нервую координатную ось — ось абсцисс. Абсциссш будем считать положительньїми для точек, лежащих севернее точки А, и отрицательннми для точек, лежащих к югу от точки А. Для определения положення точки М проведем

У

£, Е

 

Р,

Рис. С

Рис. 7

через М нормальнеє сечение таким образом, чтобьі оно пересекло меридиан начальной точки А под углом 90°. Пусть кривая зтого нормального сечения (точнеє, геодезическая линия — кривая кратчайшего расстояния на поверхности зллипсоида) изобразится на рис. 7 линией М М х. Тогда положение точки М в рассматриваемой системе координат определят длиньї следующих двух кривих на поверхности зллипсоида, которне и будут сфероидическими прямоугольннми координатами точки М:

АМі = р,

ММі = q.

Зти криволинейние координати р и q полностью определяют положение точки М на поверхности зллипсоида, если известньї геодезические координати Б и L (или другие, им зквивалентнне) начала сфероидических координат А. Система координат (р , q) имеет много общего с прямоугольной системой коор­ динат на плоскости.

Возможнн и другие системи сфероидических криволинейньїх координат в зависимости от вьібора координатних осей и порядка счета координат р и q.

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

20

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

В настоящее время в СССР принята п р о е к ц и я Г а у с с а — К р ю -

г е р а и л и

с и с т е м а п р я м о у г о л ь н ь ї х п л о с к и х п р я -

м о л и н е й н н х к о о р д и н а т в к о н ф о р м н о й

п р о е к ц и и

Г а у с с а ,

в которой производят вьічисления всех пунктов опорной геодези-

ческой сети.

 

 

 

 

§ 4. Связь между некоторьіми системами координат

1.

С в я з ь м е ж д у г е о д е з и ч е с к о шй и р о т о й В и к о -

о р д и н а т а м и х и у, о т н е с е н н ь ї м и к п л о с к о с т и м е р и -

д и а н а о п р е д е л я е м о й

то ч к и . Возьмем меридианньїй зллипс, проходящий через точку М (рис. 8). Напишем уравнение зтого зллипса

+

У - = І .

(4.1)

а2

62

 

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

єсть первая производная

сле-

 

 

 

 

довательно,

 

 

 

 

 

dy

tg (90° -f В) — —ctg В.

(4.2)

сіх

 

 

 

 

 

Виразим нервую производную

 

в функции прямоугольньїх

координат х,

у . Дифференцируя (4.1), находим

 

 

 

 

37

.

ydy

0

 

 

а2

1“ 62dx

 

откуда

dy

 

b2

х

 

 

 

(4.3)

 

dx

 

a 2

у

 

 

 

Сопоставляя вьірая^ения

(4.2)

и

(4.3),

находим

 

у_

(4.4)

х

 

Уравнение (4.4) дает внражение для геодезической широти как функции прямоугольннх координат х , у.

Чтобш найти обратную зависимость, т. е. виразить х и у как функции гео­ дезической широти В , вспомним (2.7).

На оснований (4.4) напишем

t g £

а2

 

У

1

У

(4.5)

а'1 (1

е)2

х

1 — Є2

X

 

 

и

 

у — х (1 е2) tg В.

 

(4.6)

 

 

 

Перепишем (4.1), заменив у, согласно уравнению (4.6), получим

 

CF2

Х2(1__е2)2 tg2 В

_ .

 

 

а2

'

«2 (1_ Є2)

L-

 

Решая зто уравнение относительно х, находим:

 

откуда

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

 

 

 

 

 

 

 

* =

™osB, .

(4.7)

 

 

 

У 1—e2sin2B

 

 

Для нахождения у подставляем в уравнение (4.6) найденное значение х

согласно (4.7). Получим окончательно

 

 

 

 

 

у =

а ( 1 - е 2) sin В

.

,, оч

 

 

>— .

 

(4.8)

 

 

 

у 1 —Є2 sin2 В

 

 

Из рис. 8 следует, что абсцисса точки М

 

 

 

 

 

х — ОМх — МС

 

 

 

в то же время является радиусом

г параллели, проходящей

через точку М

и имеющей широту В. Следовательно,

 

 

 

2.

С в я з ь м е ж д у г е о д е з и ч е с к о й щ и р о т о й В и г е о -

ц е н т р и ч е с к о й ш и р о т о й

Ф. Из рис. 9 легко находим вьіражение

для геоцентрической широтьі в функции прямоугольньїх координат X жу

 

 

 

 

 

 

(4.10)

На оснований формули (4.5) имеем

 

 

 

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

tg Ф = tg В (1 — е2).

(4.11)

 

 

Найдем вьіражение для

р а з н о с т и геодезической и

геоцентрической

широт В — Ф. Из формули (4.11) имеем:

 

 

 

 

 

tg В t g 0 — e2tgB,

 

 

 

sin (В — Ф) = е2sin В cos Ф

(4.12)

22

Полученная формула еще непригодна для практического употребления, так как sin Ф) вьіражается в функции В и Ф. Однако вследствие незначительности разности — Ф), не превшпающей, как увидим далее, 11,8', можно в правой части уравнения (4.12) cos Ф заменить через cos В. Рассмотрим, какая погрешность будет допущена при такой замене. Для зтого (4.12) перепишем так:

sin —Ф) = е2 sin В cos [В— (В— Ф)].

Раскладьівая cos — Ф)] по строке Тейлора, получаем:

sin (В — Ф) — е2sin В [cos5 + (В —Ф) sin 5],

sin (В — Ф) = е2 sin В cos В-{- е2 (В — Ф) sin2 В.

Второй член в правой части полученного вьіражения представляет собой малую величину порядка е4 (величина Ф), согласно формуле (4.12^. является малой величиной порядка е2).

Позтому, если в правой части уравнения (4.12), заменим cos Ф через cos В, то пренебрежем членами порядкам4. Сатойточностью

sin (В— Ф) = у е2sin 2В.

Раскладьівая sin — Ф) в ряд и ограничиваясь первьім членом, получаем приближенную формулу

 

( £ _ ф ) ^ -L p V sin

2Б,

(4.13)

 

 

допустив

снова при атом погрешность по­

Р,

 

рядка е4.

Нетрудно

видеть, что максималь­

Рис. 9

 

неє значение — Ф)" будет

при В = 45°.

 

 

В атом случае Ф)' =

11,8'.

 

 

 

Более точная формула для — Ф) имеет вид [27, стр. 24].

 

Ф)" = р" [ - ^

sin 2В -

2 (21 е2)2 sin 4В +

sin 6 В - . . . ] .

(4.14)

3.

С в я з ь м е ж д у г е о ц е н т р и ч е с к о й ш и р о т о й и к о ­

о р д и н а т а м и х и у, о т н е с е н н н м и к ц е н т р у и о с я м а л -

л и п с а .

В н р а ж е н и е

 

р а д и у с а - в е к т о р а . Обозначая

радиус-

вектор ОМ через р, на оснований рис. 9 напишем:

 

 

 

р —- У х 2-} у2;

х рсовФ;

у = рзіпФ.

(4.15)

Подставляя ати внражения в уравнение меридианного аллипса (4.1),

получаем

 

р2 cos2 Ф

. р2 sin2 Ф __ л

 

 

 

 

 

 

 

 

а2

а2 (1 —е2)

А*

 

Решаем ато уравнение

относительно р:

 

 

 

 

~0г"(1-

Є2)

[С082Ф(1 —e2) +

sin2 Ф] = 1,

 

 

 

 

аг (Се!) ( 1 - « 8С08«Ф)=1,

 

23

откуда

 

 

а У 1 — е-

 

 

 

 

.

 

 

р = ..____

 

На оснований (4.15)

имеем

У 1 — е2 cos2 ф

 

 

 

 

X=

а / і —е2 cos Ф

■•

ц =

а V 1 — е2 sin Ф

----

 

---—.,

 

У 1—

Є2COS2 Ф ’

 

У Ч — е2 соз2 Ф

(4.16)

(4.17)

Внражение для радиуса-вектора в функции геодезической широти определяется из (4.7), (4.8), (4.15), т. е.

р2 — х2-\-у2

a- cos2 В

а2 (1 —е2)2 sin2 В

 

і — е- sin2 В

 

1 — е2 sin2 В

 

 

 

 

 

Решая зто уравнение относительно р и удерживая членьї с е4, после преоб-

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

 

 

 

 

р — а ^1 —

sin2 В-\-

sin2 В — -|-е4 sin4 В — .. .^ .

(4.18)

4. С в я з ь м е ж д у п р и в е д е н н о й ш и р о т о й и и г е о д е *

з и ч е с к о й ш и р о т о й

 

В. На рис.

10 изображеньї меридианньїй зллипс

У

 

P E XP XE

и полуокружность EQEX,

необхо-

 

димая для

построения угла, являющегося

 

 

приведенной

широтой

и.

 

 

 

 

Предварительно установим связь между

 

 

ординатами точек зллииса и окружности, име-

 

 

ющими одну

и ту же

абсциссу. Например,

 

 

для точки М установим связь между отрез-

 

 

ками М 2М и М 2-М 1.

О М хМ 2 следует, что

 

 

 

Из

треугольника

 

 

 

 

(ОМ.2)2 +

(М2М 1)2 = аг,

(4.19)

 

 

 

Поскольку точка М принадлежит мери-

 

 

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

 

 

удовлетворять уравнению аллипса (4.1), т. е.

Рис. 10

 

 

 

 

(0М2)2

(М2Л/)2_ ,

 

 

 

 

 

а2

Т

Ь2

 

или

 

 

 

 

 

 

 

 

 

(0А/2)2+

(М2М)2

- а2.

 

 

(4.20)

Сопоставление виражений (4.19) и (4.20) дает

М1М2= М 2М у ,

или

ММ 2 = у = МхМь Y .

(4.21)

Для получения связи между приведенной широтой и и геодезической В имеем из рис. 10

х — a cos и

(4.22)

24

и на оснований (4.21)

 

 

 

 

у = М 1М

 

НО

 

МгМ2= a sin и,

 

позтому

 

 

 

 

 

 

 

у = a sin и ~ — b sin и.

(4.23)

Заметим,

что вираження (4.22) и (4.23) являются уравнениями зллипса

в параметрической форме.

 

Из (4.22) и (4.23) легко получаем виражение для приведенной широти и

через прямоугольнне

координати х и у.

 

 

 

— = tg U = V 1 — е2 tg и,

(4.24)

 

 

tg и

(4.25)

Но на оснований

V 1-

 

(4.5)

 

 

 

f = ( i - e 2)t gB,

(4.26)

следовательно, из (4.24) и (4.26) имеем

 

 

 

tg В (1 — е2) tg и У 1 — е2,

 

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

 

 

 

tg M = ]/l —e2tg В.

(4*27)

Получим

еще дополнительнне зависимости, которне будут

необходимьі

в дальнейшем.

 

 

Из (4.23),

принимая во внимание (2.7) и (4.8), получаем

 

sm и =

Из (4.27) пишем

1^1 <?2 sin В

(4.28)

V i — е2 sin 2 В

tg2 В

 

tg2 и

( 1 - е2)

 

 

 

1

В

= 1

tg2 и

COS2

 

1

(1 — е2) ’

откуда

cos В

V 1 — I?2 COS и ,

(4.29)

У 1 - е2 COS2 W

 

На оснований (4.9) и (4.22) можем также написать для радиуса параллели r — a cos и.

Виведем приближенную формулу разности (В — н), удобную для подсчетов tg В — tg и tg В l —e2tg В,

sin (В — и) = tg В [1 — (1 — е2),/2]. cos В cos и

25

Раскладьівая (1 — е2)1/* в ряд и заменяя cos и на cos В (допуская тем самьім ошибку на малую величину порядна е4), получаем окончательное внражение для (В и)

 

 

(В— и)" =

р"е2sin 2В.

 

 

 

 

 

(4.30)

Более точная формула для и) имеет вид

 

 

 

 

 

 

(В — u)" = p "^ns m2B

sin YB-)—

sin 6 5 — ^-sin8J5-f .. .J, (4.31)

где

_

а b

tg В — tg и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П а + Ь

tg В -j- tg u

 

 

 

 

 

 

 

5. С в я з ь м е ж д у

с и с т е м о й п р я м о у г о л ь н ь ї х

п р о -

с т р а н с т в е н н ь ї х к о о р д и н а т X, Y, Z и д р у г и м и с и с т е -

^

 

 

 

м а м и .

На рис.

11

Р В гР гВ — мери-

 

 

 

 

дианньїй зллипс, в плоскости которого

 

 

 

 

находится

точка

G начала

счета долгот

 

 

 

 

и,

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

в

зтой

плоскости

 

 

 

 

располагается

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

 

ось ОХ\

 

 

 

 

Р Е гР ХЕ — меридианньїй зллипс, на кото-

 

 

 

 

ром расположена

данная точка М и ко-

 

 

 

R,

ординатньїе

оси

Ох и

Оу.

Угол между

 

 

 

плоскостями зтих

меридианньїх

зллипсов

 

 

 

 

равен геодезической долготе L.

На рис. 11

 

 

 

 

имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X = х cos L

 

 

 

 

 

 

 

 

 

 

 

Y = ж sin L

 

 

 

(4.32)

 

 

 

 

 

 

 

Z = V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 11

 

 

 

 

 

 

ж = і/‘Х 2+ У 2

 

 

 

Далее, на оснований (4.32), (4.22) и (4.23) получим:

 

 

 

 

 

X — a cos и cos L

 

1

 

 

 

 

 

 

Y = a cos и sin L

 

1

 

 

 

 

(4.33)

 

Z = Ьsin и — a

1— е1sin и ]

 

 

 

 

 

На оснований

(4.32), (4.7)

и (4.8)

напишем

 

 

 

 

 

 

 

Х

=

а cos В

 

 

т ї

 

 

 

 

 

 

 

- ---- ■

 

 

cos L

 

 

 

 

 

 

 

 

 

У 1— Є2 sill- В

 

 

 

 

 

 

 

 

 

 

 

a cos В

 

.

т

 

 

 

 

 

(4.34)

 

 

 

У 1 — е2 sin2 В

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z а (1 —еі) sin в

 

 

 

 

 

 

 

 

 

 

 

У і —є2 sin2 В

 

 

 

 

 

 

 

 

Если заменить

в (4.34)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а2 —Ь2

 

 

 

 

 

 

 

 

 

 

 

 

а2

 

 

 

 

 

 

 

 

26

получим

 

 

a2 cos В cos L

 

 

У a2 cos2 В -f- b2 sin2 В

 

Y =

 

a2 cos В sin L

(4.35)

V a2cos2 В 4- Ь2 sin2 В >*

Z =

 

 

Ь2 sin В

 

V а 2cos2 В -f- b2 sin2 В

 

обозначая

 

 

 

 

 

 

]/ а2cos2 В +

b2sin2 В — р .

(4.36)

Формульї (4.35) перепишутся:

a2 cos В cos L

 

 

X

 

 

 

 

 

 

 

 

Y ■

a2 cos В sin L

(4.37)

 

 

р

 

Z

 

Ь2 sin В

 

 

 

 

 

 

§ 5. Главньїе радиусьі кривизни в данной точке зллипсоида

 

Через нормаль к поверхности зллипсоида можно провести бесчисленное

множество плоскостей. Зти плоскости, перпендикулярньїе к касательной пло-

скости к поверхности зллипсоида в данной точке, назьіваются

н о р м а л ь ­

н и м и .

Кривне,

образуемьіе

от

пересече-

 

 

ния нормальних плоскостей, проведенньїх в

 

 

данной точке, с поверхностью зллипсоида,

 

 

назнваются

н о р м а л ь н и м и

с е ч е -

 

 

н и я м и.

В

каждой точке существует два

 

 

взаимно перпендикулярних

нормальних се-

 

 

чения,

кривизна

которнх имеет максималь­

 

 

неє и минимальное значення;

зти нормаль-

 

 

нне

сечения

назнваются

г л а в н н м и

 

 

н о р м а л ь н и м и с е ч е н и я ми.

 

 

В

некоторой

точке М поверхности зл­

 

 

липсоида вращения главннми нормальними

 

 

сечениями, как известно из дифференциаль-

 

 

ной геометрии, являютея:

 

 

 

Pi

1) меридиональное сечение, проходящее

Рис. 12

через данную точку М и оба полюса зл­

 

представляется

липсоида Р и Р х (на рис. 12 меридиональное сечение в точке М

зллипсом РМЕгРіЕ);

2) сечение первого вертикала, проходящее через точку М и перпендику­ лярнеє к меридиональному сечению точки М. Сечение первого вертикала изображено на рис. 12 кривой WME, представляющей собой также зллипс.

Обозначим через М и N радиусьі кривизни меридиана и первого вертикала соответственно. Найдем вираження для радиусов кривизни главннх нормаль­ них сечений в функции геодезической широти В. Радиус кривизни плоской кривой, вьіражаемой уравнением вида у = / (#), определяетея формулой

к е т 1

d2y dx2

27

Применив зту формулу к меридианному зллипсу, напишем

М =

н ю т

(5.1)

 

d2y

 

dx2

(знак минус взят потому, что ^ < 0 ) . Из (4.2) имеем

Рассматривая В как функцию х, дифференцируем формулу еще раз по х

и полупаєм

1

dB

 

d2y

(5.2)

dx2

sin2 в

dx

 

 

dB

 

 

 

 

 

Для вьічисления — воспользуемся формулой (4.7)

 

a cos В

 

a cos В (1 —е2sin2 В)-1В.

 

V 1 — е2 sin2 В

 

 

 

 

Дифференцируя последнюю

формулу, находим:

dx = а { — sin Б (1 — е2 sin2 В)-1В е 2 sin В cos2 В (1 е2sin2 В)~*В} dB,

-Ц- =

a sin В (1 —<?2 sin2 В)~3в { — (1 — е2 sin2 JS)+ <?2 cos2 5},

 

= — a sin Б (1 — е2 sin2 В)~*В (1 _ е2).

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

 

 

 

 

 

 

 

с?2у

 

(1 — е2 sin 2 i?)3/z

 

 

 

dx2

 

a s i n 3 jE?(l — e2)

 

Подставляя

полученньїе вираження для

CLdC

и —— в (5.1), находим

 

 

 

 

 

CLduм

 

М:

(1 -j-ctg2 В)312 a sin3 В (1 — б2)

 

 

( 1 - е 2 sin2 В)3В

 

 

 

 

 

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

 

 

а (і —<?2)

 

 

 

 

М:

(5.3)

 

 

 

 

 

( 1 - е 2 sin2 В)

Из (5.3) ясно, что М возрастает при изменении Б от 0 до 90°.

Радиус кривизни меридианного зллипса в полюсах (при В = 90°) обозначим через с, тогда

с а ( ! —е2)

д _

 

( 1 - е 2)aB

' У 1 _ С2'

 

Принимая во внимание (2.7) и (2.5), находим

 

V 1 - е2

Ь

(5.4)

 

28

Обозначив

 

 

 

напишем

 

 

IV = у гі е2sin2 В,

(5.5)

 

 

 

а (1 —е2)

 

 

 

 

М =

(5.6)

Введем еще функцию

 

W* *

 

 

 

 

 

 

 

у = і / і + е*2соS 2B .

(5.7)

Так

как,

согласно (2.5) и

(2.6),

 

 

 

 

е.2

Є2

 

 

то

 

 

1— Є2 •

 

 

 

 

 

l+ e '2 COS2 Д

 

 

1 - е 2 sin2 5 = 1

sin25:

 

 

1 + е' 2

 

 

 

 

1 + е "

и

 

 

 

а (і+ е ,2)І/г

 

 

 

 

М*=

(5.8)

 

 

 

 

у з

 

или, согласно

(5.4),

 

 

 

 

 

 

 

 

(5.9)

W и

V — соответственно

назьіваются первой

и второй основними функ-

циями геодезической широтьі; они имеют большое значение в теории сфероиди-

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

полуоси

Заменяя в (5.3) первьій ексцентриситет его вьіражением через

и используя обозначение (4.36), формула (5.3) для М перепишется

 

М = аЧ2

(5.9')

Р

 

Для определения радиуса N первого вертикала заметим, что если сечение первого вертикала WME (рис. 12) — нормальное сечение, то параллель MQS — наклонное сечение, поскольку нормаль не лежит в плоскости зтого сечения. Указанньїе два сечения в точке М имеют общую касательную. Для доказательства зтого положення проведем в точке М касательную к параллели МТ\ зта касательная, лежащая в плоскости MQSC, перпендикулярной к меридианной плоскости ME уР гЕ Р , перпендикулярна к прямой МС, образованной пересечением зтих плоскостей. Таким образом, касательная М Т перпендикулярна к пло­ скости меридиана Р М Е гР у, позтому плоскость первого вертикала будет содержать прямую МТ; если Мп — нормаль к поверхности зллипсоида в точке М, то угол ТМп равен 90°, следовательно, М Т будет касательной и к кривой EMW.

Имея зто в виду, воспользуемся теоремой: е с л и ч е р е з т о ч к у п о в е р х н о с т и п р о в е д е н а д в а с е ч е н и я — н о р м а л ь ­ н о е и н а к л о н н о е , п р и ч е м в р а с с м а т р и в а е м о й т о ч к е з т и д в а с е ч е н и я и м е ю т о б щ у ю к а с а т е л ь н у ю , то р а д и у с к р и в и з н и н а к л о н н о г о с е ч е н и я р а в е н р а д и - у с у к р и в и з н и н о р м а л ь н о г о с е ч е н и я , у м н о ж е н н о м у

на к о с и н у с

у г л а м е ж д у п л о с к о с т я м и з т и х д в у х

с е ч е н и й.

 

29

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