Курсовая работа: Теория математической обработки геодезических измерений

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Пункты A и B исходные, а пункты C и D определяемые. В геодезическом треугольнике требуется выполнить: уравнительные вычисления измеренных углов в геодезическом четырехугольнике коррелатным способом; провести оценку точности угловых измерений по результатам уравнивания; произвести оценку точности стороны AD.

Таблица 4 - Исходные данные

№ угла

Измеренный угол

1

4803'41,6”

2

4522'49,7”

3

4227'08,4”

4

4406'22,5”

5

4926'17,3”

6

4400'06,8”

7

4030'27,4”

8

4603'05,1”

Ход решения

Перед математической обработкой геодезических измерении необходимо выполнить расчет числа независимых условных уравнений поправок r в геодезической сети, которое будет равно разнице между числом всех измерений n и числом необходимых измерений k.

r = n - k

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

При условии фигур, в геодезическом четырехугольнике возникает несколько совместно зависимых уравнений:

Их зависимость выражена следующим образом:

Нужно выбрать максимальное число независимых условных уравнений, где нельзя было бы выразить одно уравнение как линейную комбинацию других. Однако любые три комбинации уравнений независимы:

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

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

Из уравнения полюса получаем справедливое равенство , из которого выразим каждую из функций:

Отсюда имеем одно независимое уравнение полюса из пяти приведенных выражений (17).

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

k = n - r = 8 - 4 = 4

Величину r еще называют числом избыточных измерений.

Составление условных уравнений и вычисление невязок

В геодезическом четырехугольнике рассмотрим три треугольника ABC, BCD, CDA. Как известно, сумма углов в плоском треугольнике равна 180°, следовательно, сумма измеренных углов в данном треугольнике будет отклоняться от 180°. Величина отклонения называется невязкой треугольника. Условные уравнения фигур имеют линейный вид. Следовательно:

Затем составим условные уравнения поправок:

Далее, в нашем случае, запишем условное уравнение полюса, например для т. А:

Прологарифмируем обе стороны выражения (21) для приведения уравнения полюса к линейному виду:

Уравнение полюса преобразуем в уравнение поправок:

где I - изменение lg sin при измерении угла на 1 секунду.

Выражения (20) и (21) являются независимыми условными уравнениями.

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

Как указывалось, задачу нахождения минимума функции независимых уравнений [vv]=min решают способом Лагранжа, вводя вспомогательные множители независимых условных уравнений.

где n - число измерений, r - число избыточных измерений (число независимых условных уравнений).

Как видно, из системы уравнений (23) в правой части стоят нули, но для величин Xi получены результаты измерений x1,…, xn. Так как значения xi отягощены ошибками измерений, и при подстановке их в левые части условных уравнений (23), то в правых их частях получаются не нули, а невязки W.

Следовательно, уравнения (23) примут вид:

Чтобы невязки были равны нулю, необходимо их устранить за счет поправок v. В результате получаем:

Поскольку задачу на словный экстремум решают по правилам Лагранжа при помощи неопределенных множителей условных уравнений, тогда введем функцию Ф, зависящую от измеренных значений xn и неопределенных множителей . В результате получаем:

(22)

Искомые значения поправок v должны удовлетворять равенствам вида:

(i = 1, 2, …, n)

(j = 1, 2, …, r)

Далее составим систему уравнений:

В результате получаем суммы n + r уравнений и n(v)+r() неизвестных.

Если условные уравнения будут иметь нелинейный вид, то задача станет практически неразрешимой. Для этого, пренебрегая нелинейными членами разложения функций в ряд Тейлора, можно выражения (25) записать следующим образом:

Заменяя

и учитывая равенства (27), получим:

(j = 1, 2, …, r)

Это равенство (29) называют условным уравнением поправок.

В случае если переменная v связаны между собой уравнениями (25) и подставим их вместо функции в выражения (22). Для удобства вычислений множители Лагранжа обозначим так: . Множители k1, k2, …, kr - коррелаты.

Далее возьмем частные производные по аргументам vi и приравняем их нулю, получим: уравнивание коррелатный геодезический сеть

откуда общий вид уравнения поправок равноточных измерений:

где - коэффициенты условных уравнений поправок, - неопределенные множители (коррелаты). При соблюдении условия [vv]=min, получим уравнения коррелат.

Вычисление коэффициентов уравнивания фигур

Для вычисления коэффициентов уравнения фигур необходимо взять частные производные выражений (20) по аргументам . После чего получаем:

для 1-го уравнения:

для 2-го уравнения:

для 3-го уравнения:

Значения коэффициентов уравнения фигур записывают в таблицу 6 (столбцы ai, bi, ci).

Вычисление коэффициентов полюсного уравнения и его невязки W4

Таблица 5 - Вычисление невязки полюсного уравнения

Углы

Значение

Углы

Значение

4

44є06'22,5''

9,8426037

21,72

2+3

87є49'58,1''

9,99968925

0,8

6+7

84є30'34,2''

9,9980029

2,02

5

49є26'17,3''

9,8806446

18,00

2

45є22'49,7''

9,8523499

20,76

7

40є30'27,4''

9,8126120

24,62

[ ]

29,6929565

[ ]

29,6929458

Величина рассчитывается по формуле:

=

Невязка уравнения полюса вычисляется как разница между суммами правым и левым столбцом таблицы 6.

Вычисление коэффициентов уравнения полюса:

Вычисление коэффициентов условных уравнений поправок

Таблица 7 - Вычисление коэффициентов условных уравнений поправок

№ угла

Сумма S`

1

1

0

0

0,00

0,00

1,00

2

1

0

0

-2,00

2,08

1,08

3

1

1

0

0,08

0,00

2,08

4

1

1

0

-2,17

0,00

-0,17

5

0

1

1

1,80

0,00

3,80

6

0

1

1

-0,20

0,00

1,80

7

0

0

1

2,26

-2,46

0,80

8

0

0

1

0,00

0,00

1,00

[ ]

4

4

4

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

Поскольку полученная система из четырех линейных уравнений не решаема, преобразуем ее на основе способа наименьших квадратов ([vv]=min) в четыре нормальных уравнения коррелат:

Таблица 8 - Вычисление коэффициентов нормальных уравнений коррелат

a]

b]

c]

d]

F]

S`]

Невязка Wi

S=S`+Wi

[a

4

2

0

-4,09

2,08

3,99

2,20

6,19

[b

2

4

2

-0,49

0,00

7,51

-5,00

2,51

[c

0

2

4

3,86

-2,46

7,40

-3,40

4,00

[d

-4,09

-0,49

3,86

17,11

-9,72

6,66

-10,67

-4,00

[F

2,08

0,00

-2,46

-9,715

10,38

0,28

0,00

0,28

[S`

3,99

7,51

7,40

6,66

0,28

25,83

Заполненные таблицы 8 осуществляется по данным таблицы 7. Значения величин, выделенные полужирным курсивом и другим типом шрифта, не записываются. В нашем случае, значения величин нужны для контрольных вычислений.

Решение нормальных уравнений коррелат

Решение нормальных уравнений коррелат было выполнено по схеме Гаусса, ход которого приведен в таблице 9.

Значения коррелат определены, начиная с 4-го коэффициента, по формулам (37). Значения дробных соотношений берутся из таблицы 9.

Таблица 9 - Решение нормальных уравнений коррелат по схеме Гаусса

№ строк

Наименование уравнения

К1

К2

К3

К4

F

W

S

Контроль

1

N1

4

2

0

-4,09

2,08

2,20

6,19

2

E1

-1

-0,50

0,000

1,023

-0,519

-0,550

-1,547

-1,547

3

N2

4

2

-0,49

0,00

-5,00

2,51

4

(1)

-1

0

2,05

-1,04

-1,1

-3,09

5

a (3+4)

3,00

2,000

1,551

-1,039

-6,100

-0,588

-0,588

6

E2

-1

-0,667

-0,517

0,346

2,033

0,196

7

N3

4

3,86

-2,46

-3,40

4,00

8

(1)

0

0

0

0

0

9

(2)

-1,33

-1,034

0,693

4,067

0,392

10

b (7+8+9)

2,67

2,82

-1,77

0,67

4,39

11

E3

-1,0

-1,059

0,664

-0,250

-1,645

-1,645

12

N4

17,11

-9,72

-10,67

-4,004

13

(1)

-4,18

2,12

2,25

6,33

14

(2)

-0,80

0,54

3,15

0,30

15

(3)

-2,99

1,87

-0,71

-4,65

16

c (12+13+14+15)

9,13

-5,18

-5,97

-2,019

-2,019

17

E4

-1

0,567

0,654

0,221

0,221

18

N5

10,38

0

0,28

19

(1)

-1,079

-1,143

-3,214

20

(2)

-0,360

-2,112

-0,204

21

(3)

-1,174

0,442

2,912

22

(4)

-2,938

-3,387

-1,146

23

d (18+19+20+21+22)

4,828

-6,200

-1,372

-1,372

Вычисление поправок в измеренные углы

Вычисление поправок осуществляется по уравниванию поправок и их значения записывают в таблице 10.

Таблица 10 - Вычисление поправок

вi

aiK1

biK2

ciK3

diK4

vi

vi2

1

-1,043

0

0

0

-1,043

1,088

2

-1,043

0

0

-1,307

-2,350

5,521

3

-1,043

2,324

0

0,052

1,333

1,776

4

-1,043

2,324

0

-1,420

-0,140

0,020

5

0

2,324

-0,942

1,177

2,558

6,546

6

0

2,324

-0,942

-0,132

1,249

1,560

7

0

0

-0,942

1,478

0,535

0,286

8

0

0

-0,942

0

-0,942

0,888

17,684

Контроль осуществляется следующими равенствами:

=17,684

Кроме того, правильность вычислений можно проконтролировать, подставляя в условные уравнения соответствующие поправки и свободные члены. Контрольные вычисления приводятся в таблице 11.

Таблица 11 - Уравненные значения углов

№ угла

Измеренные значения xi

Попрaвкиvi

Уравненные значения Хi

1

4803'41,6”

-1,0

4803'40,6”

2

4522'49,7”

-2,4

4522'47,3”

3

4227'08,4”

+1,3

4227'09,7”

4

4406'22,5”

-0,1

4406'22,4”

5

4926'17,3”

+2,6

4926'19,9”

6

4400'06,8”

+1,2

4400'08,0”

7

4030'27,4”

+0,5

4030'27,9”

8

4603'05,1”

-0,9

4603'04,2”

Оценка точности по результатам уравнивания

Источник: https://otherreferats.allbest.ru/download/1172147/