книги из ГПНТБ / Разумов, О. С. Пространственная геодезическая векторная сеть
.pdfГ л а в а 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