Геодезические сети

Функция вида U = l1 + l2

Определить СКП U, где l1 и l2 – независимые слагаемые со случайными погрешностями ∆l1 и ∆l2. Тогда сумма U будет содержать погрешность:

∆U = ∆l1 + ∆l2.

Если каждую величину слагаемого измерить n раз, то можно представить:

∆U1 = ∆l1' + ∆l2' – 1-е измерение,

∆U2 = ∆l1" + ∆l2" – 2-е измерение,

…………………

∆Un = ∆l1(n) + ∆l2(n) – n-е измерение.

После возведения в квадрат  обеих частей каждого равенства  почленно их сложим и разделим на n:

∑∆U2 / n = (∑∆l12)/n + 2×(∑∆l1×∆l2)/n + (∑∆l22)/n.

Так как в удвоенном  произведении ∆l1 и ∆l2 имеют разные знаки, они компенсируются и делим на бесконечно большое число n, то можно пренебречь удвоенным произведением.

mU2 = ml12 + ml22;

mU = √( ml12 + ml22 ).

СКП суммы двух измеренных величин равна корню квадратному из суммы квадратов СКП слагаемых.

Если слагаемые имеют  одинаковую СКП, то:

ml1 = ml2 = m;

mU = √(m2 + m2) = √2m2 = m√2.

В общем случае:

mU = m√n,

где n – количество аргументов l.

Функция вида U = l1 - l2

mU = m√n;

mU = √( ml12 + ml22).

СКП разности двух измерений  величин равна корню квадратному из суммы квадратов СКП уменьшаемого и вычитаемого.

Функция вида U = l1 - l2 + l3

mU = √( ml12 + ml22 + ml32…)

СКП суммы n измеренных величин равна корню квадратному из суммы квадратов СКП всех слагаемых.

Линейная функция вида U = k1l1 + k2l2 + … + knln

mU = √[ (k1ml1)2 + (k2ml2)2 + … + (knmln)2],

т.е. СКП алгебраической суммы  произведений постоянной величины на аргумент равна корню квадратному из суммы квадратов произведений постоянной величины на СКП соответствующего аргумента.

Функция общего вида U = ƒ( l1, l2, …, ln)

Это наиболее общий случай математической зависимости, включающий все рассматриваемые выше функции, являющиеся частным случаем. Это  значит, что аргументы l1, l2, …, ln могут быть заданы любыми уравнениями. Для определения СКП такой сложной функции необходимо проделать следующее:

1. Найти полный дифференциал функции:

dU = (dƒ/dl1)×dl1 + (dƒ/dl2)×dl2 + … + (dƒ/dln)×dln,

где (dƒ/dl1), (dƒ/dl2), …,(dƒ/dln) – частные производные функции по каждому из аргументов.

2. Заменить дифференциалы  квадратами соответствующих СКП,  вводя в квадрат коэффициенты  при этих дифференциалах:

mU2 = (dƒ/dl1)2×ml12 + (dƒ/dl2)2×ml22 + … +(dƒ/dln)2×mln2.

3. Вычислить значения  частных производных по значениям аргументов:

(dƒ/dl1), (dƒ/dl2), …,(dƒ/dln).

И тогда mU = √[ (dƒ/dl1)2× ml12 + (dƒ/dl2)2×ml22 + … +(dƒ/dln)2×mln2].

СКП функции общего вида равна корню квадратному из суммы квадратов произведений частных производных по каждому аргументу на СКП соответствующего аргумента.

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

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

mlср. = ½ √∑d2/n

где d – разности в каждой паре; n – количество разностей.

Формула Бесселя:

mlср = ½ √∑d2/n-1

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

μ=√∑ [f2 /n]/N,

где - СКП одного угла;

f – невязка в полигоне;

N – количество полигонов;

n – количество углов в полигоне.

4. Определение дополнительных пунктов

4.1 Цель и методы  определения дополнительных пунктов

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

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

4.2 Передача координат  с вершины знака на землю. (Решение примера)

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

Исходные данные: пункт  A с координатами X_A, Y_A; пункты геодезической сети B (X_B, Y_B) и C (X_C, Y_C).

Полевые измерения: линейные измерения выбранных базисов  b₁ и b'₁; измерения горизонтальных углов ß₁ , ß'₁ , ß₂ , ß'₂ ; б , б'.

Требуется найти координаты точки P – X_P, Y_P.

Решение задачи разделяется  на следующие этапы:

Решение числового  примера

Исходные данные

Обозначе- ния	А ХА, YА	B ХB, YB	C ХC, YC	β1 β2	β2 β2`	β1 β1`	б б‘
Численные значения	6327,46	8961,24	5604,18	266,12	38o26'00"	70o08'54"	138o33'49"
	27351,48	25777,06	22125,76	198,38	42˚26'36"	87˚28'00"	71˚55'02"

 

Вычисление расстояния D_АР

Обозначе- ния	B1 B2	sinβ2 sinβ‘2	sin(β1+β2 ) sin(β‘1+β‘2)	B1 sinβ2 B2 sinβ‘2	D1 D2	D1 -D2 2D/T	Dср
Численные значения	266,12	0,62160	0,94788	165,420	174,52	0,00	174,52
	198,38	0,67482	0,76705	133,871	174,52

 

 

Решение обратных задач

Обозначения	YB YА	ХB ХА	YC YА	ХC ХА	tgαAB αAB	tgαAC αAC	sinα AB sinα AC cos αAB cosαAC	S AB S AC
Численные значения	10777,06	8961,24	7125,76	5605,08	-0,5977	7,23421	-0,51309 -0,99058 0,85833 -0,13693	3068,48
	12351,48	6327,46	12351,48	6327,46	329˚07'55"	262o07'51"		5275,51

 

Вычисление дирекционных углов α_{АР =} α_D

Обозна- чения	D	sinб sinб'	S AB S AC	sin ψ sin ψ'	ψ ψ'	φ φ'	αAB αAC	αD α'D	αD-α'D õmß
Численные значения	174,52	0,66179	3068,48	0,03950	2o15'50"	39o10'41"	329o07'55"	8o18'36"	∆α=1'30"
		0,95061	5275,51	0,03292	1o53'13"	106o11'46"	262o07'51"	8o18'37"

 

 

sin ψ = D×sinб/ S _AB; sin =174,52×0,66179/3068,48=0,03950;

sin ψ' = D×sinб'/ S _AС; sin `=174,52×0,95061/5275,51=0,03292;

ψ = arcsin 0,03950 =2 ^o15` 50``;

ψ'= arcsin 0,03292=1 ^o53` 13``;

φ = 180 ^o – (б+ ψ) = 180 ^o – (138^o33` 49``+2 <sup>o</sup>15` 50``) = 39<sup>o</sup>10` 41``

φ`= 180 <sup>o</sup> – (б`+ ψ` ) = 180 <sup>o</sup> – (71<sup>o</sup>55` 02``+1 <sup>o</sup>53` 13``) = 106 ^o11` 46``

α_D = α_AB ± φ =329^o07` 55``+ 39<sup>o</sup>10` 41``= 8<sup>o</sup>18` 36``

α_D`= α<sub>AC</sub> ± φ`=262^o07` 51``+ 106 <sup>o</sup>11` 46``= 8<sup>o</sup>18` 37``

Контроль:

(α_D – α'_D) õm_β;

где m_β –СКП измерения горизонтальных углов.

Знак «+» или «-» в формулах вычисления дирекционного угла берется в зависимости от взаимного расположения пунктов А, Р, В и С.

(8^o18` 36``-8<sup>o</sup>18` 37``) ≤ 30``

0^o00` 01`` ≤ 30``

Решение прямых задач (вычисление координат т.Р)

Обозначения	αD αD'	sinαD sinαD'	cosαD cosαD'	DcosαD DcosαD'	DsinαD Dsinα'D	∆Х - ∆Х' ∆Y - ∆Y'	ХА YА	Хp = ХА+ ∆Х Х'p = ХА+ ∆Х' Yp = YА+ ∆Y Y'p = YА+ ∆Y'
Численные значения	8o18'36"	0,14453	0,98950	172,69	25,22	∆=00,00 ∆=00,00 ∆доп=25см	6327,46	6500,15
	8o18'37"	0,14454	0,98950	172,69	25,22		12351,48	12376,70

 

Х_p = Х_А+ ∆Х,Y_p = Y_А+ ∆Y,

Х'_p = Х_А+ ∆Х',Y'_p = Y_А+ ∆Y'.

∆Х= Dcosα_D,∆Y= Dsinα_D,

∆Х'= Dcosα'_D,∆Y'=Dsinα'_D.

Расхождение координат не должно превышать величины õm_ß×p, где p=206265", m_ß – средняя квадратическая погрешность измерения угла.

Оценка точности определения  положения пункта P.

Средняя квадратическая погрешность определения отдельного пункта вычисляется по формуле:

M²_p = m²_X +m²_Y,M²_p = m²_D +(D×m_α / P)²

где m_D- определяется точностью линейных измерений, а m _α – точностью угловых измерений.

Пример: m_D =2см, m_α= 5``, тогда

M_p =√ [(0,02) ²+(170×5/2×10⁵)²] ≈ 2×10^-2 = 0,02м.

4.3 Решение прямой  и обратной засечки (по варианту задания)

Определение координат пункта прямой засечкой (формулы Юнга).

Для однократной засечки  необходимо иметь два твёрдых  пункта. Контроль определения осуществляется вторичной засечкой с третьего твёрдого пункта.

Исходные данные: твердые  пункты А(Х_АY_А); B(Х_BY_B); С(Х_СY_С).

Полевые измерения: горизонтальные углы β1, β ₂, β`<sub>1</sub>, β`_2.

Определяется пункт P.

Формулы для решения задачи:

Х_p -Х_А=((Х_B-Х_А) ctg β ₁+(Y_B-Y_А))/ (ctg β ₁+ ctg β ₂);

Х_p= Х_А+∆Х_А;

Y_p -Y_А=((Y_B-Y_А) ctg β ₁+(Х_B-Х_А))/ (ctg β ₁+ ctg β ₂); Y_p= Y_А+∆Y_А;

Оценка точности определения  пункта P.

Вычисление СКП из 1-го и 2-го определения:

M₁ =(m_β×√(S₁²+ S₂²))/p×sinγ₁;

M₂ =(m_β×√(S₁²+ S₂²))/p×sinγ₂;

 

Значения величин, входящих в приведённые формулы следующие:

m_β =5``, p=206265``; γ=73˚15,9`; γ=62˚55,7`; S₁=1686,77 м; S₂=1639,80 м; S₃=2096,62 м.

Стороны засечки найдены из решения обратных задач.

M₁ = (5``×√2,86+2,69)/(2×10⁵×0,958)=0,06м.

M₂ = (5``×√2,69+4,41)/(2×10⁵×0,890)=0,07м.

M_r = √ (M₁² +M₂²); M_r =√ [(0,06) ²+(0,07) ²]=0,09м.

Расхождение между координатами из двух определений