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

книги из ГПНТБ / Разумов, О. С. Пространственная геодезическая векторная сеть

.pdf
Скачиваний:
5
Добавлен:
19.10.2023
Размер:
5.69 Mб
Скачать

Г л а в а 6

ГЕОДЕЗИЧЕСКАЯ ВЕКТОРНАЯ СЕТЬ

§ 16. РАСПРЕДЕЛЕНИЕ ОШИБОК В ЭЛЕМЕНТАРНЫХ ПОСТРОЕНИЯХ

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

Таким образом, в пространственной прямоугольной системе координат OXyZ основным элементом сети будет являться вектор, соединяющий наблюдательные станции. Приращения координат между ними будут равны

(2.203)

Отвечающие этим уравнениям дифференциальные зависимости имеют вид

 

cos ф cos Л

— L sin ф cos Л

dL =s

cos ф sin А

L sin ф sin Л

 

sin ф

L cos Ф

(2.204)

Переходя по известным правилам от соотношений (2.204) к тен­ зору ошибок вектора L, получим

(2.205)

где Л42 ^ \ — тензорошибок непосредственно измеренных величин. При независимых результатах измерений составляющие эле­

менты тензора (2.205) равны

т2хх

cos'2 ф cos2 Am2 + L2sin2 ф cos2 Arri^ -f L2 cos2 ф sin2 Апгл2,

(2.206)

m2y =

cos2 \|)sin2 Аm'l 4- L2 sin2 ф sin2 A/n2, + L2 cos2 ф cos2 Лm?A,

(2.207)

 

m\, —- sin2 ф/?г2 -J- Lr cos2 Фm2, •

(2.208)

 

mry — cos2 ф sin Л cos Am2 -|- L2 sin2 ф sin Л cos A/n2 —

 

 

L2 cos2 ф sin А cos A/n2v

(2.209)

130

іп\г — cos Ф sin фcos Ami — cos 'Фs'n Ф cos Am^,

(2.210)

m2z = cos фsin ф sin Anil L2cos ф sin ф sin Am^.

(2.211)

Если от некоторого начального пункта прокладывается век­ торный ход, состоящий из нескольких хорд, то погрешности на­ капливаются, и тензор ошибок положения конечного пункта хода найдется как

 

 

 

' K c t f K J I ' K j n

 

 

 

 

 

 

M l = к г;-( К Л " К , ] ; K i r

 

(2 . 2 1 2 )

 

 

 

 

 

 

 

 

 

Г №

J

 

 

 

При

этом СОВОІ супиое положенне всех

пунктов

 

такого

хода

может оцениваться блочной матрнцей

 

 

 

 

 

/

МІ

/И І.

Ml

 

 

 

« г,

^

1

 

Ll

L,

LX

 

 

 

 

 

 

 

 

 

К]:

 

 

 

м

 

 

 

 

 

 

 

M l

 

 

 

 

 

 

 

М* =

A ll

1

[ « ? ] ;

n :

 

 

 

п :

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

\

М І.

Iм?j:

 

 

 

 

т

 

 

 

 

 

 

 

 

где блоки Ш 77І

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

(2.212).

 

 

 

Приняв для оценки точности положения точки формулу

 

 

 

 

т2 — in2

-f т2 4- tri1,

 

 

(2.214)

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

т\ = [тЦ + (L2m^] -f [Z.2 cos2t|w^ ] .

(2.215)

Обозначив

 

 

/n,|, = L/щ;

тл L соэфѵпд,

(2.216)

уравнение (2.215) приведем к виду

 

ml [тЦ +

I + W a \ ■ .

(2.217)

5* 131

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

т\ = (ml +ігі* т',д п.

(2.218)

Из уравнений (2.215—2.217) следует, что

в векторных ходах

ошибки положения пунктов уменьшаются по

мере уменьшения

длин хорд — сторон хода.

 

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

Рассматривая г-товын пункт векторного хода как узловую точ­

ку системы двух ходов, направленных к ней от твердых

пунктов,

можно составить для

точки

два

значения

тензора

ошибок

по

( 2.212)

 

 

 

 

 

 

 

 

 

 

М?. = [МГ2 }1

и

 

М%= [/Иг]!,,.

 

(2.219)

Величины, обратные тензорам

ошибок,

 

называют

весовыми

тензорами положения

пунктов в нашем

случаеони

будут равны

Рп =

Ш Ь ) - 1

 

И

Рі2 =

 

(МІ2Г ' .

 

(2.220)

И так как при сложении распределений

 

веса складываются,

итоговый вес уравненного положения пункта

будет

найден

так­

же в виде тензора

 

 

 

 

 

 

 

 

 

 

 

Рі0 =

Р п +Рі2.

 

 

 

 

(2.221)

Выполнив обратный переход от весового тензора к тензору

ошибок, получим

 

 

 

 

 

 

 

 

 

 

 

МЬ>

= РТо1.

 

 

 

 

(2.222)

В том случае, когда твердые

 

пункты

уравновешиваемой

си­

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

и их

положениям

свойственны некоторые тензоры

ошибок,

то

 

перед

вычислением

весовых тензоров положения определяемого пункта элементы тен­ зоров ошибок исходных пунктов нужно сложить с соответствую­ щими элементами тензоров М2п и М]2.

 

Для предварительной

оценки точности взаимного

положения

произвольных пунктов ходов

і и j

воспользуемся

теоремой

161],

о весе уравновешенного

значения

измеренного

элемента

сети

как

сумме весов двух величин — веса непосредственного измере­

ния

и веса того же элемента,

найденного косвенным

путем по ре­

зультатам остальных измерений в системе. Здесь рекомендуется следующий порядок вычислений.

132

1. Вычисляют тензор ошибок положения пункта / относительно пункта і по результатам непосредственных измерений в ходе

МІ = [мЩ .

(2.223)

2. Опустив использованные измерения, тензор ошибок взаимного положения пунктов вычисляют вторично, как сумму тензоров ошибок конечных пунктов свободных ходов, направленных к пунк­ там і и / от исходных пунктов

 

М;-= Ml+ М2=[М г]! + [л*г]"+1..

(2.224)

3.

Вычисляют весовые

тензоры

взаимного положения

пунктов

і

и j

 

 

2 X—1

 

 

 

Рц = «

р'ч = (м ;;)-1

(2.225)

 

)

и итоговый весовой

тензор

 

 

 

 

Poij =

Рч +

Р'ч-

(2.226)

 

4. Определяют тензор ошибоквзаимного

 

положения пунктов хода в принятой си­

 

стеме координат

 

 

 

 

 

 

М \ц =

Р~оі)

(2.227)

 

5. Длянахождения погрешностей состав­ ляющих элементов вектора L{L, -ф, Л), соединяющего рассматриваемые пункты хода, тензор (2.227) умножают на матри­ цу (1.99) направляющих косинусов гради­ ентов составляющих элментов этого векто­ ра

Рис. 49. Система ходов с узловым пунктом

оііУ- пм02ііпт.

(2.228)

После этого ошибки направляющих углов и модуля

вектора L

вычисляют по формулам (1.103).

 

Расчет погрешностей положения точек в простейших построе­ ниях векторных ходов с одной или несколькими узловыми точка­ ми может быть выполнен приемами метода эквивалентной замены с использованием тензоров ошибок и весовых тензоров положе­ ния пунктов.

Например, точность уравненного положения пункта В (рис. 49)

может быть найдена из таких расчетов:

М2А

и М2Аа

положения

1) подсчитывают тензоры ошибок

конечной точки А в ходах 1 и 2;

по (2.220)

и (2.221)

находится

весовой тензор эквивалентного хода 1, 2

 

 

 

 

Ра і >2~ ( М а^

1+ (М~аУ)

 

 

и тензоп ошибок точки А этого хода'М2.

2

=

Дт1 ;

 

*

 

 

-г"11о

 

133

2) вычисляют

тензоры

ошибок

положения точки

В по ходам

3, 1, 2 + 4, причем

 

 

 

f2.

 

 

 

Mi I ,2 + 4 =

 

 

 

 

^~ А І ,2 +

 

 

3) вычисляют

весовой

тензор

и тензор

ошибок

уравненного

положения точки В

 

 

 

 

 

 

Рз0— (MbJ

1+

(^Я| ,2+4')

*’

 

мъ л = Рв!.

Обратим внимание читателя еще иа одно обстоятельство. Предположим, что при производстве непосредственных изме­

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

т, Lm^ = L cos ?тА = lil

(2.229)

или с учетом (2.216)

mL = % = ,пА =

(2.230)

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

Ml = [pi]" Е33,

где E3Z— единичная матрица третьего порядка.

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

II

тf*w=”

 

то

 

 

Р‘- { и :

Е33

Е33,

+ № .

1

и

,2іі г,,21«

 

 

 

 

M l

ЕЯз-

 

№ ) 1

 

(2.231)

ходе, более

(2.232)

(2.233)

(2.234)

То есть,

погрешности

при любой геометрической

струк­

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

туре хода,

а если ход равносторонний,

то

 

 

2

2 і (п І)

г.

(2.235)

 

^ІО — \^L

-£88.

134

На основании соотношении (2.234), (2.235) и (2.227) опреде­ ляется и стандартный тензор ошибок совокупного положения всех вершин векторного хода, проложенного между твердыми пунк­ тами (рис. 50), с соблюдением условия (2.230)

АР = р|-/(\

(2.236)

где составляющие элементы квадратной матрицы АР («-го порядка) равны

/?І7. - К

і (п—і)

і(п —/)

ть

(/ > 0;

Если ход неравносторон­ ний, то

Рис. 50. Схема Е е к т о р и о г о хода

 

2 )/i

 

П

mf. =

[M i [VLli+l

тч

[PLlflPZl Ж

 

п

 

 

 

1

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

менты этого тензора, приходим к выводу, что точность

уравно­

вешенных значений

каждого из измеренных элементов

вектора

1 в таком

ходе

при постоянстве рг,

оказывается

одинако­

вой. Например, делая расчет по схеме (1.98), получим

 

mt

-а ,

=

— К "“

О Н - + 1)[Я — (і + 1)] -

 

 

 

-

2і [п -

+ 1)]) =

(- ^ -

(2.237)

 

 

 

 

 

п

 

Подмеченные

особенности векторного

хода сохраняются и в

сетях из векторных ходов произвольной формы, если в них по каж­ дой стороне хода соблюдается принцип равного влияния (2.230).

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

Кроме того, если в векторной сети ошибки положения исход­ ных пунктов описываются сферами, то такая форма сохраняется и у ошибок положения всех определяемых пунктов сети (при со­ блюдении ранее поставленных условий). А если в этой сети исход­

135

ные пункты расположены равномерно, то их ошибки распределя­

ются между элементами

заполняющей сети также равномерно,

в соответствии с весами

измеренных элементов.

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

ции и сопоставим эти данные с точностью

взаимного

положения

наблюдательных станций, соединенных вектором L

(табл. 7).

 

Приведенные

результаты показывают

некоторые

преимуще­

ства векторной

сети по сравнению с другими

системами

по­

строения.

 

 

 

 

 

 

 

 

Таблица

7-

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

 

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

 

 

Геометрические

Точность непо­

Ошибка

по­

Системы

средственных

ложении

данные

 

измерений

пункта

в м

 

 

Трилатерация...........................

(метод засечек ПВЦ) . . . .

Триангуляция ...........................

(метод засечек ПВЦ) . . . .

Триангуляция ...........................

(засечки направлений хорд) .

Векторная сеть ...........................

L = 1000 км ер = 60°

г' = 1000 км

L = 1000 км ß = 60°

L = 1000 км

тг' = ± 5 м

и =5-10 о

11 сл

о

1

а

ц = 5-10-°

15-16

14,5

9

8,5

Расчет погрешностей положения пунктов в обширных вектор­

ных сетях, учитывая их уникальный характер,

нужно

выполнять

при помощи обратной

матрицы коэффициентов нормальных урав­

нений. Приближенный расчет здесь не приводит к большой

эко­

номии времени и не

обладает

достаточной

надежностью.

На

рис. 26 показаны итоги расчета ожидаемых средних

квадратиче­

ских ошибок координат пунктов и отдельных

направлений

сто­

рон глобальной векторной сети,

построенной

по геометрической

схеме проекта космической триангуляции И. Д. Жонголовича [20]. При выполнении этого расчета не соблюдался в полной мере прин­ цип равного влияния (2.230), а были приняты следующие зна­ чения средних квадратических ошибок непосредственно измерен­ ных величин

mL = + 20 м; іп^ — + 0,5"; т"к — + 0,5" (L = 6400 км).

Расчет точности был выполнен методом обращения матрицы коэффициентов нормальных уравнений.

Данные, приведенные на рис. 26, а также ранее выполненный автором анализ родственных азимутальных построений [65], по­ зволяют сделать следующие выводы:

136

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

ошибки

положения

пунктов в векторной сети уменьшаются

в

среднем

на

22—30%,

ошибки уравновешенных значений длин сторон на 27%

(антипод­

ной стороны 9— 10 на 45%) и ошибки

уравненных

значений на­

правляющих углов ф и Л хорд в среднем на

15%);

 

 

абсолютные значения ошибок положения пунктов в сплош­ ной векторной сети в среднем уменьшаются пропорционально кор­ ню квадратному из числа векторов L, сходящихся в данном пунк­ те, независимо от углов между ними;

в свободной векторной сети ошибки положения точек воз­ растают по мере удаления от исходного пункта;

точность векторной сети повышается по мере уменьшения

длин ее сторон.

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

1)проложение отдельных траверсов через труднодоступные районы для усиления существующих астрономо-геодезических се­ тей (АГС);

2)построение жесткого каркаса для континентальных АГС,

свободного от влияния уклонений отвесных линий и имеющего абсолютную ориентировку в пространстве (при достижении точно­ сти составляющих элементов 2-10-6 такой каркас со сторонами 100—200 км мог бы служить сетью так называемого нулевого класса);

3) объединение отдельными траверсами континентальных се­ тей в единую мировую геодезическую сеть.

§17. УРАВНОВЕШИВАНИЕ ВЕКТОРНЫХ СЕТЕЙ

Взамкнутой фигуре из векторных ходов сумма составляющих

векторов должна равняться

нулю (замкнутой фигурой называть

и ту, которая опирается на

два исходных пункта, если вектору,

соединяющему эти пункты,

придать вес, равный бесконечности).

Поэтому условные уравнения поправок в векторных ходах, полу­ чаемые на основе соотношений (2.204), имеют вид

[cos ф cos Л ül ]------

[L sin ф cos Л Уф] —

 

 

Р

 

 

---- —[L cos ф sin Лол] + / с — 0

 

Р

 

 

 

[cos ф sin Aul ] -----—[L sin ф sin Лг\|,] +

(2.238)

 

P

 

 

-f- — [Lcos ф cos Лил] + / ѵ=

0

 

P

 

 

 

[sin ф vL] -f ~

[L cos фУф] + fz =

0

 

137

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-----^ [Az cos Atu]-----—[А(/иЛ] +

f K=

o

 

 

—^Vi

P

 

 

P

 

 

0

(2.239)

 

---- - [Az sin Луц,] +

— [Аллі 4- fy =

 

'

Ді/

1 j

 

 

 

 

 

 

 

 

 

L

L

LJi

P

 

 

P

 

 

 

 

 

 

 

+ — [Az ctg г[)У,|,] + f z = 0

 

 

 

 

 

 

 

P

 

 

 

 

 

 

 

где fx, fy,

fz— невязки, определяемые суммой векторов

той пли

иной уравниваемой фигуры.

 

 

 

 

 

 

 

Придав результатам непосредственных измерений веса

 

 

 

 

С

 

С

 

 

С

 

(2.240)

 

 

 

 

•) 1

/ > * = —

Р л = —

 

 

 

 

 

 

 

 

 

 

 

 

 

т Ф

 

ГП'

 

 

 

 

 

 

 

 

 

 

 

 

система условных

уравнении

поправок

решается под

условием

 

 

 

\Pl vI\ + І Ѵ ф1 + а ѵа\ = nlin-

 

 

Если

результаты

непосредственных

измерений

зависимы

между

собой,

то решение уравнений

 

(2.238)

осуществляется

по так назы­

ваемому обобщенному принципу наименьших квадратов. В этом

способе матрица коэффициентов нормальных уравнений

корре-

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

как произведение АМ\ Ат, где А1®— недиаго­

нальная корреляционная матрица непосредственных

измерений,

А — матрица коэффициентов условных уравнений поправок.

Системе уравнений (2.238) эквивалентна система

 

 

[cos ф cos A.vl ] — [sin ф cos Ауф1— [sin Аѵл[ -f- /д. —. 0

 

 

[cos ф sin Adl ] — [sin ф sin A uA[ +

[cos Aua] + fy = 0

■>

(2.241)

 

[эіпфУі.] + [соэфУф]

+ f x = 0

 

 

решаемая под условием

 

1

 

----------1

 

t о

с * Г " м

1

4

<

П11_

 

 

 

 

о

п

 

L <

J

 

п 1 А

где

 

 

 

 

 

4 = L

 

;

 

ѵл = L cos ф— .

 

P"

 

 

P

(2.242)

(2.243)

Последняя форма записи (2.241) будет удобна. В ней веса измерений имеют одинаковые размерности и потому могут быть сопоставимы. Наложение принципа наименьших квадратов при

138

решении системы (2.241)

приводит к решению системы

нормаль­

ных уравнений

 

 

 

 

 

 

AKt +

DIU +

ЕКЯ+ /, = ОI

 

 

DKi +

ВКі + FK3+ fy = 0

,

(2,244)

 

EK1 + FK-2 + CK3 + fz = oj

 

где

 

 

 

 

 

 

А =

[m2 cos2гр cos2 Л] -f-

sin2ij)cos2л |

+

[тд sin2 a J

 

В =

[ т 2 cos2 гр sin2 Л] -f- [/?/(“ sin2 гр sin2 л ]

+

[т л cos2 л]

 

С =

[ т - sin2ip] + [т ,|, соэ2ф]

 

 

 

 

D =

[ т 2 cos2 гр sin Л cos Aj +

sin2 гр sin Л cos л] —

(2.245)

[тл sin Л cos л ]

Е[ т 2 sin гр cos гр cos Л] — [m(j, sin гр cos гр cos л ]

F = [ml sin op cos гр sin Л] — [m^ sin op cos op sin л]

Если результаты непосредственных измерений в сети отвечают принципу равного влияния (2.230), то

ß = C = [pi2]; D=-.E = F = 0.

(2.246)

И тогда

 

Кг

 

I x

I

K 2

 

і

К з

Г 2 і

 

 

 

 

 

 

 

 

 

 

 

LmiI

 

[h i

 

Kl,

Uv co s

 

 

 

 

s in АІ + h s in гр;)

П /. = —

c o s Л / + /у c o s

Wl\

 

 

 

 

 

 

 

'°Ч,- =

У'

Kl,.

 

 

 

sin гр sin At —/г cos гр;)

“7“ ■у—öV Ux sin гр/ cos Л; + /

1

Lc

I ІЧ]

 

 

 

 

 

 

ѵА. =

 

 

IЧ,.

 

 

 

COS A j )

 

і-і cosap;

ГГ (/л: sin A j —

/ у

 

 

[у Ц

 

 

 

 

 

 

А так как подставляя

(2.248)

и (2.204) можно получить

 

 

А,

 

 

 

9

 

9

Пд.ѵ. —

 

 

 

 

14;

 

„ 11Ь,-

Ar

|.^2іГ >’

°Л-Ѵ/; = ~ ty

 

и Лг .

/ г

 

 

x l

\

 

 

 

 

 

[и!

(2.247)

(2.248)

(2.249)

139

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