Файл: Разумов, О. С. Пространственная геодезическая векторная сеть.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 = dR2— dRlt |
|
|
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