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

книги из ГПНТБ / Проворов К.Л. Радиогеодезия учеб. пособие

.pdf
Скачиваний:
49
Добавлен:
27.10.2023
Размер:
26.9 Mб
Скачать

Прибавив к обеим частям первого равенства 2d\ - j - 2dxd%,

а к обеим

частям второго равенства — 2d\ + 2d2d3, после небольших

преобра­

зований

получим

 

 

 

 

 

 

 

 

с\ — r\ -= 2d2

г

cos a -f- гх);

 

 

 

4~rl=2d2(czcosß

 

+ r2),

(347)

где гj =

dx — d2 и r 2 = ds

— d2

— измеренные разности

расстоя­

ний. Разделив первое равенство,(347) на второе, найдем

 

 

ci

cos a +

r j

_

^

 

 

где

с 2

cos

ß +

r 2

~

'

 

 

 

 

,,2

r 2

 

 

 

 

 

 

 

 

 

 

 

К ~ ~

Г 2 - 7 - 2

 

 

 

 

 

2

'2

 

 

 

Обозначив дирекционный угол стороны ВМ

через Т, можем напи­

сать

 

а = Г-Ті;

 

 

 

(348)

 

 

 

 

 

 

 

&=т2-т,

 

 

 

 

где T j и Т г — дирекционные углы сторон AB

и ВС. После подста­

новки (348) в (347), в результате необходимых преобразований полу­ чим

 

m s i n r + ncosr-j-Z = 0, .

 

(349)

где

т = к(уя

— уг) + (у2

— у1)

 

 

 

 

 

 

 

 

n = k(xs

— х2) -f-(ж2

^і)

 

 

 

Обозначив

 

 

 

 

 

 

 

 

 

 

-

= tgÔ,

 

 

 

 

из (349) найдем

 

 

 

 

 

 

 

 

 

sin ( Т - f ô) =

— ~ cos ô;

 

 

 

T

= —ô — aresin

 

cos ô^j .

 

 

 

Далее по формулам (348) найдем а

и ß, после

чего из

формул

(347) получим

 

 

 

 

 

 

 

 

2

2 (ci cos а + ri)

2(c2 cosß + /-2)

'

^ , J U >

По найденным расстоянию d2

и дирекционному углу Г линии Б М

по формулам (344) найдем координаты

точки

М.

 

 

Обработку результатов радиогеодезйческих

измерений на участ­

ках, охватывающих

несколько

шестиградусных зон, и при длинах

сторон, превышающих 500 км, как правило, производят на эллип­ соиде в системе географических координат. Однако в ряде случаев

320

обработку радиогеодезических измерений целесообразно выполнять на сфере, так как в этом случае вычисления можно производить по сравнительно несложным формулам сферической тригонометрии. Для этого необходимо перейти по какому-либо закону от элементов на поверхности эллипсоида к элементам на сфере. Методы проектиро­ вания эллипсоида на шар рассматриваются в курсах сфероидической геодезии.

Р и с . 176

Вычисление географических координат u и А по заданным ази­ муту А и расстоянию s (прямая геодезическая задача) проводится на сфере по следующим формулам, которые легко получаются из треугольника AMP (рис. 176, а):

sin и = cos a sin их

-f- sin or cos ux

cos Ay;

 

 

 

 

л

cos

ux .

.

. •

sin o .

.

 

, 0

r n

sin Л =

cos

-sin

A ;

sin©—-

cos и

sin Ax\

 

(351)

*•

и

1

 

 

1

v

 

 

 

 

Я, = Я,1 + со.

 

 

 

 

 

 

На рис. 176, a A — исходная

и M — определяемая точки,

а

н ­

ш л ю с . Обратную геодезическую

задачу относительно

тех же

точек

можно решать по следующим формулам, полученным из того же тре­ угольника:

& — % — Кх,

 

 

cos a «=> sin и sin ut +

cos и cos ux cos со;

 

sin

і - ( и і - и )

 

 

t g { ( Ax + A2) =

f

ctg - ;

(352)

COS—(«і + м)

21 Заказ 129

321

 

 

 

 

 

1

/

 

 

 

 

 

tg

1

2

 

cos — (щ — и)

 

 

 

 

( Л - Л

) = _

|

 

ctg f . ;

 

 

 

 

 

 

 

S = O Ä .

 

 

 

 

 

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

координат

точки

M по двум азимутам Аг

т& А*

с исходных пунктов А и В найдем сначала

углы а и ß сферического

 

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

(см. рис. 176, 6")

 

 

 

 

 

 

 

 

 

а = А12

 

Ai,

 

 

 

 

 

 

 

ß =

^ 2

^ 2 1 .

 

 

 

 

 

После этого из того же треугольника по формуле четырех

элемен­

 

тов найдем стороны s2 и s2

 

 

 

 

 

 

 

 

 

etg at

sin o1 2 — cos o1 2

cos щ 4- sin wx ctg u2,

 

 

 

ctg o2

sin cr12 = cos o"12 cos u.2 4- sin u2

ctg wa;

 

(353)

 

 

 

 

a-Tf

 

?

 

 

 

 

 

Теперь остается получить искомые координаты точки M по форму­

 

лам (351). Координаты можно найти по любой

из сторон sx

или s2;

 

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

 

наты по каждой из этих сторон.

 

 

 

 

 

 

 

Для получения координат точки M по двум измеренным расстоя­

 

ниям st и s2 (рис. 176, б) необходимо сначала найти угол а и азимут

 

Ах стороны AM. Соответствующие формулы имеют вид

 

 

 

 

 

cos а =

COS 0"2— COS Ci cos 0"l2

(354)

,.

 

 

-.

 

sm

ffi2

;

 

 

 

 

 

sin

'

v

'

 

 

 

Ai =

Ai2

 

а.

 

 

 

 

 

После этого координаты точки M найдем по формулам (352). Решение на сфере задачи фазового зонда, когда измерены разности расстояний Гі = Si — s2 и r 2 = s3 — s2 с определяемой точки до трех исходных точек А, В и С (рис. 177), выполняется по способу, анало­ гичному для этой же задачи на плоскости. Из сферических тре­

угольников АМВ и ВМС найдем

cos Оі — cos о 1 2 cos а2 4- sin а1 2 sin а2 cos а;

 

 

cos 0 3 = cos а 2 3

cos ö 2

+ sin a23

sin a2

cos ß.

 

 

Подставив

в эти выражения О і = о 2

+ Рі И

Os — о 2 ~т- Рг> г Д е

о- =

р =

после соответствующих

преобразований

найдем

 

 

cos рх — cos 0 1 2

=

tg а2

(sin pi 4- sin o i

2

cos а);

(355)

 

 

cos p2 cos o 2 3

=

tg a2

(sin p2 4- sin a2

3

cos ß).

 

322

Положив

к-

COS Pi COS 0"i2

 

COS p2—COS 023 '

разделим первое из уравнений (355) на второе. После этого получим

 

 

sin pi +sin cos а

£

 

 

2

 

3 COS ß

 

 

 

Sin p

+Sin 023 cos ß

 

 

 

2

 

 

 

Положив

 

 

а = A — Аг;

(356)

 

 

 

где А — азимут

линии ВМ,

и,

 

сделав необходимые

преобразова­

 

ния, найдем

 

 

 

 

 

m sin А 4- п cos А 4-1 = 0,

(357)

 

где обозначено

 

 

 

 

 

т — к sin 0 2 зs u l ^ 2 s

i n °"i2 s

i

n -^îî

 

n = & sin 0 2 S cos Л 2

— sin cr12 cos

AL;

 

I — к sin p2 —sin p b Уравнение (357) с помощью под­

становки t g ô = -^- можно привести

к

виду

 

 

 

 

 

 

s i n ( ô - M H

cos ô = 0.

 

 

Р и с . 177

 

 

Откуда найдем

 

 

 

 

 

 

А =

arcs m ( i c

o s ô ) -

(358)

Из

(356) найдем а

и ß, а затем из (355) получим

 

 

=

COSp! — COS0l2

COS р2 — COS 023

 

 

° 2

sin рх - j - S i l l 012 cos

ce

Sin p2

+ Sin 023 cos ß

 

Теперь, зная o 2 и А из решения прямой геодезической задачи [формула (351)1, найдем координаты точки М.

Способы решения геодезических задач на поверхности эллип­ соида рассматриваются в сфероидической геодезии. Приведем один из способов. Прямоугольные геоцентрические координаты X, У, Z определяемой точки M по двум измеренным величинам — расстоя­ ниям или азимутам — от двух исходных точек А я В (см. рис. 176, б) можно найти из следующих трех уравнений:

 

z2 + z/24-ez2 —1 = 0;

 

 

ккх

+ 1{у 4- mxz 4- nxz% 4- рх

= 0;

(359)

к2х

4- l.j/ 4- m2z 4- re2z2,-f р 2

= 0.

 

21*

323

Здесь

где а и е' — большая полуось и второй эксцентриситет земного эл­ липсоида. Коэффициенты в написанных уравнениях для расстояния (длины хорды) si от точки А или В до точки M равны

ki

х і ' і h =

Уі> m i = zi'i

 

 

*i = Л Pi = - 1 + j ( 4 + e ' % z * ) •

( 3 5 9 a )

 

* =

(1,

2)

 

 

Для азимута A t нормального

сечения с точки А или В на

точку

M выражения для коэффициентов

следующие:

 

 

Ni

 

 

 

 

h =

Ус cos

А{

гХіУі sin Л,-;

 

li = —~- x i c o s

A-£ #A- s i n

A ;

 

7^ = (я?-)-у')sin Л,; raf = 0;

'

(3596)

pl

= (f*{x\

+

y\)zl&\n.Al.

 

 

Уравнения (359) получены путем преобразования выражений для

хорды и азимута

нормального

сечения,

имеющих вид

 

(x-XiY + iy-yif

+

i

z -

Z i Y - ^ Y ^ Q ;

 

cos Bi

\ ( l - e a ) t g Д +

 

 

 

'

sin 5; cos (L — Li)

 

ctgA- =

!=

s i

n

^

t )

 

( 3 5 9 B )

 

 

 

 

 

 

 

 

Выразив а: и у из второго и третьего уравнений (359) через z и под­ ставив найденные выражения в первое уравнение (359), получим урав­ нение, из которого найдем неизвестное z. Если обе измеренные вели­ чины являются азимутами нормальных сечений, то полученное урав­ нение будет квадратным, решение которого не представляет затруд­ нений. Если же одна или обе измеренные величины длины хорд, то полученное уравнение будет четвертой степени. Корни его можно найти численным методом (методом хорд или касательных). Выбор необходимых корней можно производить с помощью грубого графи­ ческого построения.

При вычислении на шаре радиуса R следует во всех формулах положить

N = a = R; e's = 0.

324

Тогда уравнения примут вид

 

x2 + j / 2 + z , 2 - l = 0 ;

 

кхх + lxy + mxz + px = 0;

(360)

к2х + l2y - f m2z + р2 = 0

 

и отыскание из них координат во всех случаях сведется к решению квадратного уравнения.

Приведенные выше формулы позволяют получить координаты определяемой точки при наличии только необходимого числа измере­ ний. Если для определения одной точки имеются избыточные (незави­ симые) измерения или же обработке подлежит несколько точек, соста­ вляющих геодезическую систему, то возникает необходимость урав­ нивания измеренных величин. Наиболее просто уравненные коорди­ наты определяемых точек можно получить посредственным способом уравнивания. Рассмотренные ниже методика и формулы являются общими, независимо от вида измеренных величин и от поверхности, на которой производится уравнивание. Способ сводится к отысканию из уравнивания поправок Ьх и ог/ (при уравнивании на сфере — попра­ вок ей и ÔÀ и при уравнивании на эллипсоиде — поправок и 8L)

к приближенным координатам х0

и у0 (соответственно и0

и Х0 на

сфере

или

В0

и L 0 на эллипсоиде).

Обозначим

поправки,

выраженные

в линейной мере, для определяемой точки с номером і через

и г|/,

т. е. положим:

 

 

 

 

 

 

1)

на

плоскости

 

 

 

 

 

 

 

 

1* = ож,;

г\і = Ьуі,

 

 

(361)

2)

на

сфере

 

 

 

 

 

 

 

 

h = jrbuü

ѣ

= ^рі8%1;

 

(361)»

3)

на

эллипсоиде

 

 

 

 

 

 

 

 

6і = !

§

; % = ^ б х -

 

(Збі)"

где В — радиус шара; (1) =

р" : М;

(2) =

р : N (М и N — радиусы

кривизны меридиана и первого

вертикала); щ и Bt — широта

опре­

деляемой точки (на сфере или на эллипсоиде). Исходя из известных дифференциальных формул, напишем уравнения погрешностей в сле­ дующем виде:

для расстояния, измеренного между точками 1 и 2,

 

— cos Ах\х — sin Axr\x — cos А2%2 sin А2т\2 + ls = vs;

(362a)

для направления или азимута, измеренного с точки 1 на точку 2,

— t 4- т ~ s i n Л г І і T -

c o s

+

T - s i n АъЪъ —

*12

*12

 

 

*12

 

— — cos А2ц2

+ lN

= vN.

(3626)

325

Здесь А Х и А 2 — прямой и обратный азимуты при уравнивании на сфере или на эллипсоиде, или же дирекционные углы при уравни­ вании на плоскости; £ — поправка ориентирования, общая для всех направлений на данном пункте. Если уравниваются азимуты (дирек­

ционные углы), то £ = 0. Свободные члены уравнений

погрешностей

будут

 

 

Zs = s*2 s12; In — А\

ALT

 

где s° 2 и А \ — приближенные длины и

азимуты

(дирекционные

углы), найденные из решения обратных геодезических задач по при­ ближенным координатам первой и второй точек, a s1 2 и A T — изме­ ренные длины и азимуты. Если измерялись направления, то для полу­ чения по ним «измеренных» азимутов достаточно к каждому направле­ нию прибавить приближенный азимут направления, принятого на станции за начальное (нулевое). Во втором уравнении (362) приведен­ ная длина геодезической линии для эллипсоида или шара заменена

длиной самой

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

линии, что несущественно повлияет

на

результаты

уравнивания.

 

 

Ясно, что если одна из точек (первая или вторая) будет исходной,

то

соответствующие слагаемые

в уравнениях (362) исчезнут. Если

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

t = vN.

Уравнения погрешностей для измеренных разностей расстояний (в способе фазового зонда) или разности направлений (в способе изме­ ренных углов) следует получать как разности соответствующих урав­ нений для расстояний и направлений. Веса для уравнений (362) следует находить по известной формуле

к

Р~ та '

где к — коэффициент пропорциональности, a m — средняя квадратическая ошибка измеренной величины: расстояния ms, направления или азимута mN, разности расстояний тг или угла та. Если уравни­ ваются однородные величины, например углы, измеренные с одина­ ковой точностью, то веса можно считать равными единице.

Составив уравнения погрешностей для всех измеренных величин, из решения соответствующих им нормальных уравнений можно найти величины I и и для всех определяемых пунктов. Тогда вероятнейшие координаты этих пунктов найдутся по формулам:

на плоскости

х = х0

+ 1; у = у0

+ ч;

(363)

на сфере

 

 

 

u = u0+Çb

І - ^

+ С ^ ч

(363)'

326

на эллипсоиде

В = В0 + (І)£, L = L 0 +(2) sec В r).

(363)"

По результатам уравнивания изложенным способом на эллип­ соиде можно найти поправки dX, dY и dZ к приближенным прямо­ угольным геоцентрическим координатам X , У и Z. Путем дифферен­ цирования выражений

X = N cos В cos

L ; Y = N cos В sin L

Z=N

(1 — e 2 ) s i n ß

найдем

е

/ 2 « z

dX=

 

^у^-dB — ydL;

dY

=

-1^-dB+xdL;

dZ^^yfdB.

Заменив здесь dB и dL через £ит] [см. формулы (361)], получим

Vi

где

 

 

d = l/"X 2 + F 2 ; 8 =

l + e'2.

 

Оценка точности измеренных и уравненных величин производится

по общим правилам способа наименьших квадратов.

 

Изложенный способ позволяет совместно уравнивать

разнородные

величины — стороны, направления и их разности.

 

При определении рассмотренным способом положения только

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

будут иметь следующий вид:

all + bti\ + lt

= vl

(364)

(последние два члена уравнений (362)). В этом случае будет только два нормальных уравнения, из которых -и определятся поправки g ит). Рассмотренный способ может применяться и при отсутствии избы­ точных измерений, когда для определения положения одной точки измерено только две величины и, следовательно, уравнений будет только два. Напишем эти уравнения

aiÊ + M + * i = 0;

(365)

327

Из уравнений (365) сразу же найдем поправки £ и и . Изложенный способ целесообразен во всех случаях, когда заранее известны при­ ближенные координаты определяемой точки (например, с карты или по линиям положения).

Уравнивание радиогеодезических измерений в любых геодезиче­ ских построениях можно проводить также по методу условных изме­ рений. При работе на настольных вычислительных машинах этот метод во многих случаях будет проще. Однако при вычислении с по­ мощью ЭВМ метод посредственных измерений оказывается выгоднее из-за стандартности и общности вычислительных операций.

Рассмотренными выше способами далеко не исчерпываются суще­ ствующие, весьма разнообразные, способы решения специальных задач по данным радиогеодезических измерений, наиболее подробное изложение которых дается в специальной монографии В. А. Поле­ вого [30]. /

В процессе уравнивания и окончательных вычислений может воз­ никнуть необходимость решения треугольников. Если обработка производится в проекции Гаусса, то измеренные элементы (азимуты, углы или расстояния) предварительно редуцируют на плоскость. Решение треугольников производится по формулам плоской тригоно­ метрии. Если известны углы, то стороны находят по теореме синусов. Если даны стороны, то вычисление производят по формуле

или по формуле

где

p = j(a + b + c).

При вычислениях на сфере или на эллипсоиде решение треуголь­ ников выполняют также по формулам плоской тригонометрии (на основании теоремы Лежандра). При этом, если известны сферические (сфероидические) углы, то для вычисления сторон из каждого угла предварительно вычитают одну треть сферического избытка. Если же известны длины сторон, то, найдя плоские углы, прибавляют к каж­ дому углу одну треть сферического избытка. Сферический избыток е находят по приближенной формуле-

 

absin С „

, „ „ „ ,

g =

2Д2

Р •

(367)

С вполне достаточной для радиогеодезических измерений точно­ стью решение треугольников по способу Лежандра можно выполнять для треугольников с длинами сторон до 500 км. При сторонах, пре­ вышающих 500 км, треугольники решают по формулам сферической тригонометрии, исправляя углы в необходимых случаях сфероидическими поправками.

328

§ 46. О Л И Н И Я Х П О Л О Ж Е Н И Я

ПР И РАДИОГЕОДЕЗИЧЕСКИХ РАБОТАХ

Как известно, для определения точки на поверхности эллипсоида необходимо и достаточно измерить две, а в пространстве — три неза­ висимые величины. Таковыми могут быть как расстояния или напра­ вления, так и их простейшие функции, например разности. Если измерить только одну из перечисленных выше величин, то можно построить линию, для всех точек которой значение измеренной вели­ чины будет одно и то же. Такую линию называют линией положения

или линией уровня. Например, линией положения измеренного рас­ стояния является окружность с радиусом, равным этому расстоянию. Искомая точка поверхности определится только при наличии двух таких линий-, как точка их пересечения. Построив заранее сетку из двух семейств линий положения, можно определить местоположение точки и ее координаты графически, в пересечении соответствующих линий.

На рис. 178 изображена азимутальная сетка, рассчитанная на определение точек по измеренным азимутам (дирекционным углам). Линиями положения этой сетки являются лучи, проведенные от двух исходных точек через 10°. На рис. 179 показана круговая сетка, пред­ назначенная для определения точек по расстояниям до двух исход­ ных точек. Линии положения этой сетки представлены концентриче­ скими окружностями с радиусами, изменяющимися через 2,5 мм. Для лучшей читаемости каждое семейство линий положения на сетке вычерчивают различными цветами (на рис. 178 и 179 линии одного семейства показаны сплошными, а другого — пунктирными линиями).

Область определения точек на поверхности эллипсоида, или на какой-либо проекции эллипсоида, можно представить как плоское скалярное поле, каждой точке которого соответствует некоторая ска­ лярная величина и (под и можно понимать любую измеренную

329

Соседние файлы в папке книги из ГПНТБ