Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
geokniga-геодезия-попов-вн-чекалин-ви-2007.pdf
Скачиваний:
3614
Добавлен:
16.03.2015
Размер:
39.37 Mб
Скачать

s

 

=

s

 

sin γ 1 ;

 

 

 

 

 

1

 

 

0 sin β

α 1

=

α 0 ± 1800 + η 1

 

 

 

 

 

 

1

 

s2

=

s0

sin γ 1 sin γ 2 ;

α 2

=

α 0 ±

2 1800 + η 1 − η 2

(16.42)

 

 

 

 

 

sin β 1 sin β 2

α 3

=

α 0 ±

3 1800 + η 1 − η 2 + η 3

 

s

 

=

s

 

sin γ 1 sin γ 2 sin γ 3

 

3

0

;

 

 

 

 

 

 

 

sin β 1 sin β 2 sin β 3

 

 

 

 

 

Представим уравнения (16.41) через поправки в углы в линейной форме:

å

å

æ

ϕ

x

ö

 

ç

 

 

ν

 

 

ç

γ

 

 

 

 

è

i

ø 0

æ

ϕ

y

ö

 

ç

 

 

 

ν

ç

γ

 

 

 

 

è

i

ø 0

γi

γi

+å

+å

æ

ϕ

x

ö

 

ç

 

 

ν

 

 

ç

β

 

 

 

 

è

i

ø 0

æ

ϕ

y

ö

 

ç

 

 

 

ν

ç

β

 

 

 

 

è

i

ø 0

βi

βi

+å

+å

æ

ϕ x

ö

 

ç

 

 

ν

 

 

ç

η i

 

 

 

è

ø 0

æ

ϕ

y

ö

 

ç

 

 

 

ν

ç

η i

 

 

 

è

ø 0

η i + Wx = 0

(16.43)

η i + Wy = 0

Свободными членами Wx и Wy в уравнениях (16.43) являются приближенные значения искомых функций (16.41), вычисленные по измеренным горизонтальным углам с использованием равенств (16.42).

Из уравнений (16.41) найдем частные производные (коэффициенты условных уравнений поправок) и подставим их в уравнения (16.43). При этом поправки углов выражаются в секундах, свободные члены – в дециметрах, а разности координат – в километрах. Получим условные уравнения поправок:

- для координат х (при уравнивании углов):

å

xn xi

ctg γ iν γ i

å

xn xi

ctg β iν β i

å

yn yi

(±ν η i ) + Wx = 0;

(16.44)

 

 

20,206

20,206

 

20,206

 

 

 

 

- для координат у (при уравнивании углов):

å

yn yi

ctg γ iν γ i

å

yn yi

ctg β iν β i +

å

xn xi

(±ν η i ) + Wy = 0

(16.45)

20,206

20,206

20,206

 

 

 

 

 

 

 

Здесь: xn и yn - координаты последнего пункта (для рис. 16.6 xn = x5 , yn = y5); xi и yi - координаты текущего пункта i ходовой линии, проходящей через вершины промежуточных углов ηi ; νβi и νγi - поправки связующих углов β и γ (угол γ лежит против исходной стороны треугольника); νηi - поправки промежуточных углов ηi (записываются со знаком «плюс» для левых по ходу углов, со знаком «минус» - для правых по ходу углов).

§ 151. Методы решения систем линейных нормальных уравнений

151.1.Способ последовательной подстановки

Рассмотрим решение системы линейных уравнений способом последовательной подстановки на следующем примере. Имеется система линейных уравнений с четырьмя неизвестными х1, х2, х3, х4 :

1. 4х1 2х2 + 3х3 – 2х4 – 1 = 0

(16.46)

2. - 2х1 + 5х2 – 2х3 + х4 – 6 = 0

3.3х1 2х2 + 3х3 – 4х4 + 8 = 0

4.- 2х1 + х2 – 4х3 + 2х4 + 4 = 0

Шаг 1. Выразим в уравнении 1 (16.46) х1 через остальные неизвестные,

х1 = 0,5 х2 – 0,75 х3 + 0,5 х4 + 0,25,

(16.47)

396

и подставим его значение в уравнения 2, 3 и 4. Получим новую систему уравнений с тремя неизвестными:

2. 4 х2 – 0,5 х3 – 6,5 = 0

(16.48)

3. – 0,5 х2 + 0,75 х3 – 2,5 х4 + 8,75 = 0

4. – 2,5 х3 + х4 + 3,5 = 0

 

Шаг 2. выразим в уравнении 2 (16.48) х2 через остальные неизвестные (в

данном случае – через х3),

(15.49)

х2 = 0,125 х3 + 1,625,

и подставим его в уравнения 3 и 4 (16.48). Получим систему из двух линейных уравнений с двумя неизвестными:

3. 0,6875 х3 – 2,5 х4 + 7,9375 = 0

(16.50)

4. – 2,5 х3 + х4 + 3,5 = 0

Шаг 3. Выразим в уравнении 4 (16.50) х3 через х4 ,

(16.51)

х3 = 0,4 х4 + 1,4 ,

и подставим его значение в уравнение 3 (16.50). Получим

(16.52)

- 2,225 х4 + 8,9 = 0,

откуда х4 = - 8,9 / - 2,225 = + 4.

 

Из уравнения (16.51) х3 = 0,4 х 4 + 1,4 = + 3.

 

Из уравнения (16.49) х2 = 0,125 х 3 + 1,625 = + 2.

 

Из уравнения (16.47) х1 = 0,5 х 2 – 0,75 х 3 + 0,5 х 4 + 0,25 = + 1.

 

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

151.2.Способ матричных преобразований

Рассмотрим пример, приведенный в 151.1.

Для решения системы линейных уравнений матричным способом необходимо переписать их в виде

1.

4х1 2х2 + 3х3 – 2х4 = +1

(16.53)

2.

- 2х1 + 5х2 – 2х3 + х4 = + 6

3.3х1 2х2 + 3х3 – 4х4 = - 8

4.- 2х1 + х2 – 4х3 + 2х4 = - 4

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

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

397

 

+ 4.....

2..... +3.....

2.....+1.....

+ 4

 

 

 

 

 

 

 

 

 

 

 

2.....

+5.....

2.....

+1.....+ 6.....

+8

 

 

 

+3

2

+3

4 8

8

 

(16.54)

 

 

 

 

 

 

 

 

 

 

 

2

+1 4 + 2 4

7

 

 

 

 

 

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

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

Из математики известно, что результат решения не изменится, если:

-любую строку матрицы поменять местами с другой строкой;

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

Решение матрицы сводится к образованию т.н. треугольной матрицы вида

a .....a

 

.....a

 

.....a

 

.....

 

m

 

.....t

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

3

 

4

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0......b .....b .....b

.....

 

 

 

m

 

 

 

.....t

 

 

 

 

 

 

 

 

 

 

 

 

 

2

3

4

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

(16.55)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0......0......c3 .....c4

......

m3

 

 

 

.....t3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0......0......0......d4

.....

 

 

m4

 

 

.....t4

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом получается система линейных уравнений

a1k1 + a2k2 + a3k3 + a4k4 = m1

b2k2 + b3k3 + b4k4 = m2 (16.56) c3k3 + c4k4 = m3

d4k4 = m4

для, например, четырех линейных уравнений. Из последней строки находят значение k4 :

k4 =

m4

(16.57)

d4

 

 

398

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

k в исходные уравнения (16.53).

Проследим решение на приведенном примере.

Шаг 1. Образовать 1-й нулевой столбец в строках 2, 3 и 4 матрицы (16.54). Для этого умножим 2-ю и 4-ю строки на (+2), а 3-ю строку – на (-4/3). Получим:

 

+ 4

2

+3

2

 

 

 

 

 

+1

 

 

 

 

+ 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

+10 4

+ 2

 

 

 

+12

 

 

+16

 

 

 

 

 

 

4

+8 / 3 4

+16

 

+32 / 3

 

+32 / 3

 

(16.58)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

+ 2

8

+ 4

 

 

 

8

 

 

 

 

 

 

14

 

 

 

 

 

 

.........

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Затем последовательно сложим 2-ю, 3-ю и 4-ю строки (16.58) с первой строкой этой матрицы:

 

+ 4 2 +3 2

 

 

 

 

+1

 

 

 

+ 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

..0

+8

1

0

 

 

+13

 

 

 

 

 

+ 20

 

(16.59)

 

 

 

 

 

 

 

 

 

 

 

 

........

 

 

 

0

+ 2 / 3 1 +16

 

+35 / 3

 

+ 44 / 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

..0

0

5

+ 2

 

 

 

7

 

 

 

 

10

 

 

 

 

 

 

.........

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Шаг 2. Образовать 2-й нулевой столбец в строках 3 и 4 (16.59). При этом в примере для строки 4 нет необходимости в преобразованиях, поскольку в ней на второй позиции уже имеется ноль. В связи с этим достаточно преобразовать только 3-ю строку. Для этого умножим ее на (-12)

 

+ 4 2 +3 2

 

 

 

 

+1

 

 

+ 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(16.60)

 

..0

+8

1

0

 

 

+13

 

 

 

 

+ 20

 

 

 

........

 

 

0

8

+12

40

 

140

 

176

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

..0

0

5

+ 2

 

 

 

7

 

 

 

 

10

 

 

 

 

 

 

 

.........

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

399

и сложим полученную строку со 2-й строкой той же матрицы:

 

+ 4 2 +3 2

 

 

 

 

 

+1

 

 

 

+ 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

..0

+8

1

0

 

 

 

+13

 

 

 

 

 

+ 20

 

 

 

 

........

 

 

 

 

0

0

+11

40

 

127

 

 

156

 

(16.61)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

..0

0

5

+ 2

 

 

 

7

 

 

 

 

 

 

10

 

 

 

 

 

 

.........

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Шаг 3. Образовать нулевой 3-й столбец (16.61) в строке 4, для чего требуется умножить его на (+2,2)

 

+ 4 2 +3

2

 

 

 

 

 

+1

 

 

 

 

 

+ 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

..0

+8

1

0

 

 

 

+13

 

 

 

+ 20

 

 

 

 

0

0

+11

40

 

127

 

156

 

 

 

......

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

..0

0

11

+ 4,4

 

15,4

22

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(16.62)

и сложить со строкой 3 этой же матрицы (16.62):

 

+ 4 2 +3

2

 

 

 

 

 

+1

 

 

 

 

+ 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

..0

+8

1

0

 

 

 

+13

 

 

 

 

 

 

+ 20

 

 

 

..........

 

 

 

 

0

0

+11

40

 

127

 

 

 

156

 

 

 

 

........

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

..0

0

0

35,6

142,4

178

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(16.63) В результате система линейных уравнений (16.53) преобразуется к виду:

1. +4k1 – 2 k2

+ 3 k3 - 2 k4 = + 1

 

2.

+8 k2

- k3

= + 13

(16.64)

3.

 

+11 k3 - 40 k4

=

- 127

4.

 

– 35,6 k4

=

- 142,4

 

Из уравнения 4 (16.64) находим k4 = +4. Из уравнения 3 подстановкой в него значения k4 находим k3 = +3. Из уравнения 2 подстановкой k3 (коэффициент при k4 равен нулю) находим k2 = +2. Из уравнения 1, после подстановки

400