Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

2965

.pdf
Скачиваний:
5
Добавлен:
15.11.2022
Размер:
2.69 Mб
Скачать

 

Поскольку угловые невязки f

 

в полигонах вычисляются по тому же

принципу, что и истинные погрешности

, то в формуле (58) можно заменить

на

f , а число погрешностей n на число невязок (полигонов) N. Тогда

 

 

 

 

 

 

 

 

 

 

pf 2

 

 

 

 

 

 

 

.

(76)

 

 

N

 

 

 

 

 

 

 

 

 

 

Примем вес измерения i-го угла (i = 1, 2, …, n)

j - го полигона (j = 1, 2, …,

N) равным единице, т.е. p ij 1 p .

Согласно выражению (66) вес суммы углов j - го полигона

 

 

p j

 

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

nj

 

 

 

 

Тогда формулу (76) можно записать в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(77)

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

Величина в формуле (77) есть средняя квадратическая погрешность

измерения одного угла, так как

p

2

 

 

1 , т.е.

m .

m2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3. В сети N полигонов геометрического нивелирования получены невязки f h1 , fh2 , …, fhN в суммах превышений полигонов (ходов), периметры кото-

рых L1 , L2 , …, LN (в километрах).

Невязки fhi в полигонах представляет собой погрешности соответст-

вующих сумм превышений. Принимая вес превышения в полигоне длиной 1 км равным единице, имеем

pкм 2 1 , т.е. mкм .

mкм2

Число таких приведенных на 1 км превышений в j-ом полигоне равно Lj , а вес суммы превышений в этом полигоне

 

 

 

 

p j

1

.

 

 

 

 

 

 

 

 

 

Lj

 

 

Как и в предыдущем примере, заменим в формуле (58) на fh , n на N

и

p j

 

1

. Тогда средняя квадратическая погрешность единицы веса опреде-

 

 

 

 

Lj

литься по формуле

41

 

 

f 2

 

 

 

 

 

 

 

h

 

 

 

 

 

 

L

 

 

 

 

 

 

,

(78)

 

N

 

 

 

 

 

 

 

 

где L – длина хода, км.

В случае тригонометрического нивелирования превышения вычисляют

как

h d tgv ,

где d – горизонтальное проложение линии; – угол наклона линии. При малых углах наклона можно принять

h d vрад ,

где рад – угол наклона линии в радианной мере.

В соответствии с формулой (65) обратный вес превышения

 

1

 

 

1

 

v2рад

 

1

d 2 .

(79)

 

 

 

 

 

 

рvрад

 

ph

pd

 

 

 

 

 

 

 

 

Пренебрегая первым членом выражения (79) ввиду его малости, имеем

1

 

 

1

 

 

d 2 .

 

(80)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ph

 

рvрад

 

 

Считая измерения углов наклона равноточными и приняв

рvрад 1, полу-

чим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

1

.

 

(81)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

d 2

 

 

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

 

 

f

2

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

2

 

 

 

 

D

 

 

,

(82)

 

N

 

 

 

 

 

 

 

 

где D – периметр полигона, км.

§16. Оценка точности вычислений с приближенными числами

Равномерное распределение случайных погрешностей измерений. В

некоторых случаях законы распределения элементарных погрешностей могут

42

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

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x)

 

при

 

 

x

 

 

 

 

2

,

(83)

 

 

 

 

 

 

 

 

 

f (x) 0

при

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где – максимальное значение погрешности.

Как следует из графика функции (83), показанного на рис. 2, вероятность погрешности на всем интервале ; одинакова, т.е. f x c . Поскольку площадь прямоугольника под графиком будет равна

Рис.2. График равномерного распределения случайных погрешностей

f (x)dx 1,

 

 

 

 

 

 

 

то c 1,

откуда

c

1

и

 

 

 

 

 

 

2

 

 

 

 

 

 

 

f x

1

.

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

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

Математическое ожидание

M X

1

 

xdx

x

2

 

 

1

2 2 0 .

 

 

|

 

 

 

 

 

2

 

4

 

4

 

Дисперсия может быть получена как

D X M X 2 M (X ) 2 .

 

 

 

2

 

 

1

 

2

1 1

3

1 3

 

 

3

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

X

 

 

 

2

x dx

2

 

x |

6

 

 

 

 

 

.

 

 

 

 

 

 

 

3

 

 

 

 

 

 

3

Тогда стандарт распределения, т.е. средняя квадратическая погрешность округления

 

 

 

 

 

,

 

m (x)

D(x)

(84)

 

 

 

0

 

 

3

 

 

 

 

 

 

 

 

 

где – абсолютное значение предельной погрешности округления.

Средняя квадратическая погрешность округления. В процессе вычис-

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

43

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

Процесс округления приближенного числа до некоторого разряда состоит в отбрасывании всех цифр, стоящих правее этого разряда. В вычислительной практике принято округлять числа по правилам Гаусса. Чтобы округлить число до n значащих цифр, отбрасывают все цифры, стоящие справа от n -ой значащей цифры. При этом следует соблюдать следующие требования:

если первая из отброшенных цифр меньше 5, то оставшиеся десятичные знаки сохраняются без изменения (например, 72,4≈72);

если первая из отброшенных цифр больше 5, то к последней оставшейся цифре прибавляется единица (84,56≈84,6);

если первая из отброшенных цифр равна 5 и среди остальных отброшенных цифр имеются ненулевые, то последняя оставшаяся цифра увеличивается на единицу (64,51≈65);

если первая из отброшенных цифр равна 5 и все остальные отброшенные цифры являются нулями, то последняя оставшаяся цифра

сохраняется неизменной, если она четная, и увеличивается на единицу, если она нечетная (например, 52,50≈52, а 53,50≈54).

Случайные погрешности округления обладают следующими свойствами: а) предельная погрешность округления 0,5 единицы последнего зна-

ка;

б) положительные и отрицательные погрешности округлений равновозможны;

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

Случайное значение числа a может находиться с равной вероятностью в пределах от a 0,5t до a 0,5t , где t 2 − значение единицы последнего со-

храненного разряда. В единицах округляемого разряда возможны следующие погрешности округления:

−0,5; −0,4; −0,3; −0,2; −0,1; 0,0; +0,1; +0,2; +0,3; +0,4; +0,5.

Если их рассматривать как истинные погрешности, тогда средняя квадратическая погрешность округления будет

 

 

 

2

 

 

 

m0

 

 

 

 

0,3 .

n

 

 

 

 

 

 

 

 

Принимая mo 0,3 и 0,5 , можно записать соотношения между ними:

1,7mo mo 3

или, как уже было доказано ранее,

m

 

.

(85)

 

 

 

 

 

 

o

3

 

44

В общем случае выражение

(85)

 

можно записать как

 

mo

 

 

 

 

t

 

 

 

0,3t .

(86)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

2

 

3

 

 

 

 

 

 

 

Оценка точности арифметических действий с приближенными чис-

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

Погрешность округления алгебраической суммы приближенных чисел

определится как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

my mo

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n ,

 

 

(87)

где n − число членов суммы или разности.

 

 

 

 

 

 

 

 

 

 

 

 

С учетом выражения (85)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

my

 

 

 

n

.

 

 

 

(88)

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Предельная погрешность

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M у пред 3my

3n .

(89)

Относительная погрешность алгебраической суммы чисел

 

 

M у пред

 

 

 

 

 

 

 

 

 

 

 

 

 

3n

.

 

(90)

 

 

у

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

Оценка умножения и деления двух приближенных чисел в относительной

форме выражается соотношением

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

my

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 .

 

 

 

(91)

 

 

y

x

 

 

 

 

 

 

 

 

 

 

o

 

 

 

 

 

 

 

 

 

 

 

 

Для произведения нескольких чисел предельная погрешность

 

 

M y пред

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3n .

(92)

 

 

у

 

х

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При n 2 выражение (92) можно записать как

My пред 2 .

ух

Относительная предельная погрешность произведения n сомножителей

45

M y пред

 

 

 

 

.

(93)

 

у

 

х

 

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

M y пред

n

 

;

M y пред

 

1

 

 

.

(94)

у

х

у

n

х

 

 

 

 

 

 

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

§17. Понятие о прямой и обратной задачах теории погрешностей измерений. Принцип равных влияний

При практическом использовании теории погрешностей измерений часто решают две принципиально разные задачи.

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

С помощью обратной задачи можно установить средние погрешности

отдельных измерений для обеспечения заданной точности некоторой функции измеренных величин y f x1, x2 , , xn .

Воснову решения обратной задачи положено выражение

M y F m1 , m2 ,

, mn ,

(95)

где M y − погрешность функции (допустимое значение конечного результата), m1, m2 , , mn − искомые средние квадратические погрешности измерений.

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

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

 

 

2

 

df

2

2

 

 

df

2

2

 

 

 

 

df

2

 

2

 

 

M y

 

 

 

mx1

 

 

 

 

mx2

 

 

 

 

 

 

mxn .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx1

 

 

 

 

dx2

 

 

 

 

 

dxn

 

 

 

 

Приняв, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

df

 

2

2

 

 

df

2

2

 

 

 

 

df

 

2

 

2

 

M y2

 

 

 

 

 

 

mx1

 

 

 

 

mx2

 

 

 

 

 

mxn

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

dx1

 

 

 

 

 

 

dx2

 

 

 

 

 

 

dxn

 

 

 

 

46

получим

df

mx

 

df

mx

 

 

df

mx

 

 

M y

 

.

(96)

dx1

dx2

 

dxn

 

 

 

 

 

 

 

n

1

 

 

2

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Пример. Превышение между точками, определенное методом тригонометрического нивелирования, рассчитано по формуле h D sin v . В результате измерений получено наклонное расстояние D 200 м и угол наклона v 5 . Величина превышения h должна быть получена со средней квадратической погрешностью M h 0,10 м. С какой точностью должны быть измерены наклонное

расстояние D и угол наклона v ? Погрешность функции h D sin v

2

 

dh 2

2

dh 2

2

 

2

2

 

2

 

2

 

mv2

 

M h

 

 

 

mv

 

 

 

mD

sin

 

v mD

D

 

cos

 

v

 

 

.

 

 

 

 

 

 

2

 

dD

 

dv

 

 

 

 

 

 

 

 

 

 

При условии, что влияние погрешностей измерений расстояния и угла наклона на величину средней квадратической погрешности превышения одина-

ково, можно записать

 

 

 

sin2 v mD2

D2

cos 2 v

mv2

 

 

 

M h2

,

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M h2

2 sin2 v mD2 ;

M h2

2D2

cos 2 v

mv2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D cos v

mv

.

M

h

 

 

2 sin v m ;

M

h

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mD

 

 

 

M h

;

mv

 

 

 

M h

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 sin v

 

 

 

2 D cos v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя в полученные выражения численные значения известных величин, получим:

mD

 

 

Mh

 

 

 

 

 

0,10 м

0,81

м;

 

 

 

 

 

 

1,41

0,087

 

2 sin 5

 

 

 

 

 

 

 

 

mD

 

0,81 м

 

 

 

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

D

200 м

 

 

250

 

 

 

 

47

 

 

 

Mh

 

 

0,10 м 206265

 

mv

 

 

 

 

 

 

1,41 200 м 0,996

73 .

2 D cos 5

 

 

 

 

 

Следовательно, в рассматриваемом примере для обеспечения требуемой точности определения превышения ( M h 0,10 м) измерение угла наклона мож-

но выполнять теодолитом типа Т30, а расстояние измерять с помощью нитяного дальномера.

48

Раздел второй Уравнивание геодезических сетей сгущения и съемочных сетей

Глава 5. ОБЩИЕ СВЕДЕНИЯ ОБ УРАВНИВАНИИ ГЕОДЕЗИЧЕСКИХ

СЕТЕЙ

§ 18. Сущность уравнительных вычислений. Метод наименьших квадратов

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

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

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

Избыточные измерения позволяют получить определяемые величины несколькими путями, при этом значения одной и той же величины будут различны между собой, т.е. возникает неоднозначность решения задачи. В результате уравнивания каждая определяемая величина должна иметь одно и то же значение независимо от пути ее вычисления. Таким образом, уравнивание сети есть вычислительный процесс по распределению невязок, т.е. нахождение системы поправок, которая доставляет наилучшие значения измеренных и искомых величин, а также оценка их точности. Наилучшими значениями будем считать величины, наиболее близкие к действительным (истинным) их значениям.

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

 

 

 

 

 

 

 

2

 

(97)

 

 

 

 

 

 

p

 

min,

где

 

2

2

2

2

; pi – вес i-го измерения;

i – вероятнейшая

p

 

p1 1

p2 2

pn n

поправка i-го измерения;

i=1, 2, …, n.

 

 

 

Данный принцип был предложен в 1806 г. французским математиком А.

Лежандром для обработки результатов наблюдений и получил развитие в работах К.Гаусса, русского математика А.А. Маркова и др.

49

Уравнительные вычисления с применением принципа наименьших квад-

ратов называют уравниванием по методу наименьших квадратов. Теоретически доказано, что соблюдение условия (97) при уравнивании приводит к наилучшим оценкам определяемых величин. В этом случае уравнительные вычисления называют строгими.

Строгие способы уравнивания обеспечивают соблюдение трех требо-

ваний:

а) выполнение всех геометрических условий в сети; б) нахождение наиболее надежных значений измеренных величин под

условием p 2 min ;

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

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

Для уравнивания геодезических построений используют параметрический и коррелатный способы. При коррелатном способе поправки отыскивают непосредственно к измеренным величинам, а при параметрическом – к их функциям. Оба этих способа уравнивания являются строгими и приводят к равноценным результатам. Выбор того или иного способа уравнивания определяется в основном числом исходных данных и избыточных измерений, а также с учетом использования оптимального алгоритма решения задачи. В настоящее время уравнивание геодезических сетей выполняют на ЭВМ по параметрическим алгоритмам и программам, основанным на итерационных методах, позволяющих решать системы уравнений с некоторой заранее заданной точностью с помощью сходящихся бесконечных процессов, т.е. последовательными приближениями.

Наряду со строгими методами существуют приближенные (упрощенные) способы уравнивания, в которых соблюдаются геометрические условия в сети, а другие требования выполняются частично.

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

50

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]