Файл: Разумов, О. С. Пространственная геодезическая векторная сеть.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 29.10.2024

Просмотров: 56

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

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

В соответствии с принципом наименьших квадратов наиболее надежное значение вектора X определяется соотношением

X = — (ATPA)-'ATPL = — QÄrPL,

(1.74)

а точность искомых неизвестных характеризуется матрицей

 

М2= р2 (АТРА)-' = p2Q,

(1.75)

где ц — ошибка единицы веса, Q — матрица весовых коэффициен­ тов (обратная матрица коэффициентов нормальных уравнений по­ правок неизвестных).

Уравнение (1.75) имеет тензорную форму и операции с вели­ чинами М2 подчиняются правилам тензорной алгебры. А так как

при определении

положения точек в пространстве искомыми неиз­

вестными являются

их координаты,

то,

следуя

 

предложению

Ю. А. Гордеева [15], матрицу М2 мы будем

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

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

точек. Для

отдельно взятого пункта

 

 

 

 

 

 

 

 

 

(mU

m

12

m

13

 

 

m*

m

xz

 

 

 

 

2

xy

 

 

Ml =

m

 

 

 

 

 

m2

m2

(1.76)

m22

m*

=

m~

m*

mz

 

 

xy

yy

yz

 

i m.

m

23

 

 

m

32

2

 

,nz2

m2

m2 ,

 

 

 

 

 

ззу

 

' XZ

yz

22/ к

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

Тензор ошибок М2 непосредственно связан с матрицей Р коэф­ фициентов корреляции соотношением

М2= triPrn,

(1.77)

где in — диагональная матрица, составленная из средних квадра­ тических ошибок искомых неизвестных

mл

 

 

m,

 

 

 

m =

 

 

 

 

(1.78)

 

\ о

 

 

m„

 

 

1

Pl2 •

• Pln\

 

Р =

Рзі

1

• P2/1

(1.79)

 

 

 

 

 

vPnl Рл2

• •

• 1 /

 

2* 35


В свою очередь

 

 

 

 

9

_

Qu

 

 

 

 

 

Р

тч

 

 

(1.80)

 

tnnnijj

VQnQjj

 

 

 

 

 

 

 

 

 

В тех случаях,

когда

для определения положения

отдельно

взятой точки измеряется

лишь

необходимое

число

элементов

и,

матрица частных

производных

от

координат

искомого пункта

по

результатам измерений

 

дх

-

дх

дх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

д«і

ди.

дч3

 

 

 

 

 

В =

 

ду

 

дц

ду

 

 

 

(1.81)

 

 

dui

 

дп.

ди3

 

 

 

 

 

 

 

 

 

 

 

 

\

dz

 

dz

dz

 

 

 

 

 

 

 

 

"ІнГ

ди3

 

 

 

 

взаимообратима с матрицей

А

частных

п

водных

от

прибли-

женных значений

измеренных

величин

по

координатам этого

пункта

 

 

дщ_

 

дих

дііі

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/

дх

~дІ

dz

 

 

 

 

 

 

А =

 

ди. ди. du,г

 

 

(1.82)

 

 

дх

~ду

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

\

ÈL*.

 

ди3

du2

 

 

 

 

т. е.

 

\

дх

 

ду

dz

 

 

 

 

 

 

 

В — Л-1.

 

 

 

 

(1.83)

 

 

 

 

 

 

 

Тогда в силу известного соотношения

 

 

 

 

 

 

(.АТРА)~1=

А~1Р~Х(АТ)~

 

 

(1.84)

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

 

Ml =

вмівт,

(1.85)

где Ml — средний

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

ошибок (корреляцион­

ная матрица) непосредственных измерений.

 

Формула (1.85)

позволяет

делать предварительные расчеты

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

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

36


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

Основные операции с тензорами ошибок следующие.

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

тіа = аМі а — т~хх cos2a -j- тууcos2ß -f- m~zzcos2у +

2m^ cos а cos ß +

-f 2mv2cos а cos у + 2nfyzcos ß cos y,

(1.86)

где

 

а—■ (cos а cos ß cos y ) .

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

М \ - а 0 = Х а 0,

(1.87)

где X— главные значения тензора.

Уравнение (1.87) равносильно системе трех однородных урав­ нений

т2хаох + т% аоу + inxz aoz = Хаох

 

 

1ПХуй0Х"T" Чіцу Иду

 

myZ0Q — Хйду

I

 

(1.88)

mxzaox + rrCyzа0у + m2za0z = XaoZ

 

 

 

которая имеет решение,

отличное

от нуля только

в том

случае,

когда ее определитель равен нулю,

т. е.

 

 

 

 

(т 2хх X )

 

 

 

 

 

9

 

 

 

 

т ху

 

 

Шхг

 

 

 

Мху

(Щу -- а.)

 

/72.

= 0.

(1.89)

 

yz

 

 

 

 

2

 

/72.

 

 

(mL — X)

 

 

 

fflxz

 

yz

 

 

 

 

 

 

 

Или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m uu

m .■yz

m i x

9

 

X 3 — X 2 (n fxx +

niyy -f m : z) +

 

 

т - х г

 

X

 

mx

2

m l z

+

 

 

 

 

 

 

 

ГПгг

m 2zz

 

 

2

2

 

m xx

m .■xy mlz

 

 

 

+

Ш х х

т Ху

+

 

2

 

2

 

 

(1.90)

2

2

tTL■vuxy

myy

rjlyz

0 .

 

 

т ху

г п уу

)

 

2

2

2

 

 

 

 

 

 

 

tTlxz

fflyz

tTLzz

 

 

 

Поскольку тензор ошибок всегда имеет симметричный вид, то уравнению (1.90) отвечают три вещественных положительных кор-

37


ня — главные значения тензора. В этих новых составляющих тен­ зор ошибок записывается в виде

 

 

( К

0\

 

 

(ЛІ“) '=

К

.

 

(1-91)

 

 

\0

V

 

 

Контролем вычислений здесь может служить равенство

 

К + 'Ч +

=

тІх +

тІи+

mlz-

(1.92)

Тензору (1.91) и сопоставляется поверхность эллипсоида

с полу­

осями

 

 

 

 

 

а = Ѵ К і

ь = Ѵ \ \

с =

|/Л Г ,

(1-93)

направления которых соответствуют главным направлениям тен­ зора. Эти направления определяются из решения системы одно­ родных уравнений (1.88) после подстановки в них соответствую­ щего корня л,- кубического уравнения (1.90).

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

тіа =

cos cpj + Xocos cp., -|- X.j cos cp3)

(1.94)

где cp — углы, составленные заданным направлением с направле­ ниями полуосей эллипсоида ошибок.

3. Преобразование тензора (1.76) при переходе к другой ор гональной координатной системе осуществляют по правилу

М І = иМ 2к ПТ,

(1.95)

где П — матрица направляющих косинусов координатных осей но­ вой системы относительно старых осей.

Для тензора ошибок (1.75) совокупного положения нескольких точек справедливо равенство

/ П

 

\

,п г т

 

°

\

/ П

0

/ Пг

 

м\ =

 

м%

 

 

 

\ о

’ п

1

Ѵ о

т

1

 

пт

где блоки П составлены из матриц направляющих косинусов (1.95),

иих число равно количеству трансформируемых пунктов.

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

38


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

 

 

 

 

2

 

 

 

 

 

 

 

Піг

m^lZl

 

 

 

 

 

г а .

•л-'іУ

 

 

 

 

 

 

2

 

О

°

 

 

 

 

щЧ,/,1/,

 

 

tn\

 

 

 

 

 

9

4hV

 

 

 

ІПг,.ѵ,

9

/пЦ

.2

 

М 12

т 2,1/,

9

tnг,1/2

 

тх„

' ф ,

mJA

m . v „2,І / 2

 

 

 

/n.vaA.a

 

 

 

т.ц9*„х. Щп!1\

,n!hZi

/wjU.

inІ/2І/2

 

 

 

 

2

2

«l :

m2-1/2

 

 

 

тZ..V, /?г2»//,

mZnZ

Под

совокупной погреш-

 

 

 

ностыо

положения

одного

 

, Z

 

пункта

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

друго­

 

 

 

 

 

го будем

понимать

ошибку

 

 

 

вектора Li2,

 

их соединяющего

 

 

к

(рис. 14).

 

 

 

 

у ч

 

Тогда, дифференцируя век­

 

 

 

 

 

торное уравнение, получим

 

'

,1

dLn = dR2dRlt

 

 

X /1

 

или

 

 

 

 

 

а

dLx = dXa — dXu

■.v,z2

\

2

\

mi/122

 

2

 

ffZz.z.'

(1.97)

2

m .v.2z2

 

9 "

mZz-

\ и \ г ^

 

І V

_

!

 

-4)'

dLy = dK2— dYx,

dLz = dZ2— dZy.

Рис. 14. Положение вектора L в пространстве

Принимая во внимание, что в рассматриваемом случае

Lv = Х, — Х ъ Ly —-Y* — Yu Lz = Z%— Zx.

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

динат)

получим,

взяв частные производные функции L12 по коор­

динатам точек, т. е.

 

 

 

 

 

 

 

 

 

 

 

dLx

.

dLx

— o-

dLx

=

0;

dLx

= 1;

dLx

П-

dLx

=

0;

дХ1

'

avx

dZx

 

 

dX2

dY3

— u,

dZ2

dL>

0.

dLy

д1У

0;

dLy

П-

dLy

 

dLy

=

0;

ÖAT

 

dYx

dZx

- =

dX,

— u,

öYn ■= l;

dZ2

dLZ

Q.

dLz

o-

dL*

-

l;

dLz

- = 0;

dLz

= 0;

dLz

 

1

öXi

 

dYx

dZx

 

 

 

dX2

dYo

 

dZ2

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( ~

l

 

0

0 1 0 0\

 

 

 

 

 

 

 

 

T = \

0

—1 0 0 1 0 ,

 

 

 

 

 

 

 

 

 

V

0

 

0 - -1 0

0 1/

 

 

 

 

 

39