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

v7 = 80,0538 80,0540 = − 0,2мм .

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

§160. Оценка точности уравненных элементов и их функций

160.1.Общие положения

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

Среднюю квадратическую погрешность любой величины можно определить по формуле

M i

=

μ

 

 

 

 

 

,

(16.225)

 

pi

 

 

 

где μ - средняя квадратическая погрешность единицы веса; p – вес определяемой величины.

Обе величины, входящие в правую часть формулы (16.225) обычно неизвестны, поэтому по материалам уравнивания находят как значение средней квадратической погрешности единицы веса, так и вес оцениваемой (уравниваемой) величины. Здесь следует иметь в виду, что вес измеренной и вес той же, но уравненной величины – не одно и то же. Тем более и веса функций уравненных величин, зависящих от входящих в нее аргументов со своими весами. Отношение весов Р уравненных значений измеренных величин к весам р измеренных величин примерно равно отношению общего

числа измерений n к числу необходимых измерений k, т.е.

 

 

P

n

 

(16.226)

 

p

k

 

 

 

При определении погрешности единицы веса можно использовать формулу (3.44) из теории погрешностей, в которую вместо истинных погрешностей или уклонений от среднего значения подставляют значения полученных невязок W:

 

 

 

 

 

 

μ =

[pW 2 ]

,

(16.227)

n

 

 

 

 

 

где n – число невязок однородной величины, равных числу условных уравнений. Часто для оценки μ используют все возникающие условные уравнения (N > r).

Пользуясь материалами уравнивания, погрешность единицы веса легко можно найти по формуле

 

 

 

 

 

 

 

μ =

[

2 ]

,

(16.228)

n −

k

 

 

 

 

 

где v – значения поправок к измеренным величинам, имеющим вес pi; n – число использованных поправок, т.е. число измеренных однородных величин

495

[рν 2 ]

(углов, расстояний, превышений, пролетов и т.п); k – число необходимых измерений (n – k = r – число избыточных измерений).

При уравнивании коррелатным способом величину можно получить несколькими методами:

-по значениям поправок, полученных по результатам уравнивания (т.н. способ при помощи таблиц коэффициентов);

-по формуле

[рν 2 ] = [Wk ] ,

(16.229)

где W – невязки; k – значения коррелат;

- по схеме решения нормальных уравнений: к нормальным уравнениям коррелат (16.89) добавляют еще одно уравнение:

[qa1a1 ]k1 + [qa1a2 ]k2 + ... + [qa1ar ]kr + W1 = 0

…………………………………………

[qa j a1 ]k1 + [qa j a2 ]k2 + ... + [qa j ar ]kr + W j = 0 ,

(16.230)

………………………………………….

 

[qar a1 ]k1 + [qar a2 ]k2 + ... + [qar ar ]kr + Wr =

0

W k + W k

2

+ ... + W k

r

+ [pν 2 ] = 0

 

1 1

2

r

 

 

в котором значение [рν 2 ]

играет роль неизвестного. Т.е. получается система

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

При уравнивании параметрическим способом значение [рν 2 ] тоже можно определить несколькими путями:

-по значениям поправок (с использованием таблицы коэффициентов);

-по формуле

[pν 2 ] = [ pa1 l]τ 1 + [ pa2 l]τ 2 + ...+ [ pak l]τ k + [ pll] ;

 

(16.231)

- в схеме решения нормальных уравнений, присоединив уравнение

(16.231) к системе нормальных уравнений (16.172):

 

 

[ pa1a1 ]τ 1 + [ pa1a2 ]τ 2 + ... + [ pa1ak ]τ k + [ pa1l] =

0

……………………………………………

(16.232)

[ pai a1 ]τ 1 + [ pai a2 ]τ 2 + ... + [ pai ak ]τ k + [ pai l] = 0

 

……………………………………………

 

[ pan a1 ]τ 1 + [ pan a2 ]τ 2 + ... + [ pan ak ]τ k + [ panl]

=

0

[ pa1l]τ 1 + [ pa2l]τ 2 + ... + [ pak l]τ k + [ pll] [рν 2 ]

=

0

Данная система имеет (k+1) уравнение со столькими же неизвестными, т.е. решается исключением полученных значений поправок τj к параметрам tj. Уравнения решают по разработанному алгоритму (пример решения таких уравнений представлен далее).

160.2. Оценка точности при уравнивании коррелатным способом

Чтобы выполнить оценку точности любой величины, последнюю не-

обходимо представить в виде функции результатов измерений:

 

F = Ф(x1 , x2 ,..., xn ) ,

(16.233)

496

а после этого можно применить формулу теории погрешностей измерений для определения обратного веса функции:

 

 

 

1

=

1

é æ

Ф ö 2

1

ù

,

(16.234)

 

 

 

 

 

 

 

è

х ø

 

ú

 

 

 

 

P

P

 

р

 

 

 

 

 

 

= ê ç

 

 

 

 

 

 

 

 

F

 

Ф

ê

 

 

 

ú

 

 

 

 

 

 

 

 

 

ë

 

 

 

û

 

 

где

æ

Ф ö

- частные производные

 

функции по аргументам

xi; pi – вес

ç

 

 

 

è

х ø i

 

 

 

 

 

 

 

 

 

 

 

аргумента (измеренной величины).

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

[qa1a1 ]ρ 1 + [qa1a2 ]ρ 2 + ... + [qa1ar ]ρ r + [qa1 F] = 0

……………………………………………

 

[qar a1 ]ρ 1 + [qar a2 ]ρ 2 + ... + [qar ar ]ρ r + [qar F ] = 0

(16.235)

 

 

 

[qa1 F ]ρ 1 + [qa 2 F ]ρ 2 + ... + [qa r F ]ρ r + [qFF ] =

1

,

 

 

 

P

где ρj – некоторый неизвестный параметр;

F

 

 

[qa1 F ] =

q1a11 F1

+

q1a12 F1 + ... =

L1

 

 

[qa2 F ] =

q2 a21 F1

+

q2 a22 F2 + ... =

L2

 

 

………………………………..

(здесь L – cвободные члены новых уравнений);

[qFF ] = [q1 F1 F1 ] + [q2 F2 F2 ] + ... = L(r + 1) - свободный член последнего уравнения. Из решения системы уравнений (16.235), например, по алгоритму

Гаусса, получают

N((rr+)1)( r +1)

=

1

 

(16.236)

PF

 

 

 

Для каждой из функций

соответственно

получают свое значение

N ((rr+)1)( r +1) , т.е. необходимо выполнить столько решений уравнений (16.235),

сколько в них содержится определяемых функций F.

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

Выберем функции, оценку точности которых мы будем выполнять в данном примере:

F1 = H3 - H1 = h4- h7

F2

= H1

= H P10

+ h1

F3 = H2

= HP10 + h2

(16.237)

F4

= H3

= H P 20

- h9

F5

= H4

= HP30

- h8

497

Последние четыре функции содержат только одно значение h, т.е. их средние квадратические погрешности равны средней квадратической погрешности соответствующего превышения, поскольку HP10, HP20 и HP30 – исходные высоты (по условию задачи являются безошибочными).

Составим для обработки результатов измерений табл. 16.67, в которую занесем значения коэффициентов и обратных весов из примера 153.2, а также значения частных производных функций:

 

æ

F1

ö

= + 1

æ

F1

ö

 

æ

F2

ö

= + 1

æ

F3

ö

= + 1

æ

F4

ö

= - 1

æ

F5

ö

= - 1

 

ç

 

ç

= - 1

ç

 

ç

 

ç

 

ç

 

 

ç

h

 

 

 

;ç

h

 

 

;ç

h

 

 

;ç

h

 

 

;ç

h

 

 

;ç

h

 

 

 

è

4

ø

 

 

è

7

ø

 

è

1

ø

 

è

2

ø

 

è

9

ø

 

è

8

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(16.238)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 16.67

№№ qi

a1

 

a2

 

a3

 

a4

 

a5

 

F1

 

F2

 

F3

 

F4

 

F5

 

изм.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0,42

+1

 

 

 

 

 

 

 

 

 

+1

 

 

 

+1

 

 

 

 

 

 

 

3

2

0.68

-1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+1

 

 

 

 

 

0

3

1,08

+1

 

-1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

4

0,39

 

 

 

+1

 

 

 

 

 

 

+1

 

+1

 

 

 

 

 

 

 

 

 

3

5

1,32

 

 

 

+1

 

+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

6

1,03

 

 

 

 

 

 

+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

7

1,51

 

 

 

 

 

 

+1

 

+1

 

 

 

-1

 

 

 

 

 

 

 

 

 

1

8

1,72

 

 

 

 

 

 

 

 

+1

 

+1

 

 

 

 

 

 

 

 

 

-1

 

1

9

1,19

 

 

 

 

 

 

 

 

-1

 

 

 

 

 

 

 

 

 

-1

 

 

 

-2

(10)

W

-7

 

+18

 

-16

 

+6

 

+17

 

 

 

 

 

 

 

 

 

 

 

 

(11)

[qa1

2,18

 

-1,08

0

 

 

0

 

0,42

 

0

 

0,42

 

-0,68

 

0

 

0

 

1,26

(12)

[qa2

(1,08)

2,79

1,32

 

 

0

 

0,39

0,39

 

0

 

0

 

0

 

0

 

3,81

(13)

[qa3

(0)

 

(1,32)

3,86

 

1,51

 

0

-1,51

 

0

 

0

 

0

 

0

 

5,18

(14)

[qa4

(0)

 

(0)

(1,51)

 

4,42

 

1,72

-1,51

 

0

 

0

 

1,19

 

-1,72

 

7,33

(15)

[qa5

(0,42)

(0,39)

(0)

(1,72)

 

2,53

0,39

 

0,42

 

0

 

0

 

-1,72

 

5,87

(16)

[qF

 

 

 

 

 

 

 

 

 

 

 

 

1,90

 

0,42

 

0,68

 

1,19

 

1,72

 

 

Втаблице в скобках записаны значения коэффициентов, находящихся слева от диагональных. Строки, не относящиеся к номеру измерения (10 – 16), указаны в скобках.

Встолбцах F для значений [qa и [qF вычисления производят по формулам (16.235):

- столбец F1 : [qa1 F1 ] =

q1a11 F1 + q2 a12 F1

=

0 ; [qa2 F1 ] = 0,39 и т.д.;

- столбец F2 : [qa1 F2 ] =

0,42 ; [qa2 F2 ] =

0

и т.д.

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

В табл. 16.68 заносят диагональные коэффициенты в строки N, а у всех остальных коэффициентов меняют знак на противоположный. На противоположный знак следует изменить и значения ∑.

498