Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
учебники / Клюшин Е.Б., Михелев Д.Ш., Барков Д.П. - Практикум по прикладной геодезии - 1993.pdf
Скачиваний:
379
Добавлен:
26.04.2015
Размер:
13.91 Mб
Скачать

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

створов

и

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

общего

и

последовательных створов.

6 . Исходя из совокупности величин рассмотренных ошибок, на­ иболее точной будет схема последовательных створов, наименее точной— схема частных створов.

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

§ 41. Статистический анализ результатов геодезических измерений при наблюдениях

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

выполняемых

в

сложных (горных,

внутрицеховых, температурных

и др.) условиях

среды.

 

Поэтому,

прежде чем приступить

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

геодезических измерений данного цикла по известным формулам, основанным на законе нормального распределения ошибок измерений, необходимо установить: 1) в какой степени фактическое распределение ошибок измерений отличается от этого закона; 2 ) будут ли измерения

внутри одного цикла и между

циклами равноточными;

3) имеется

ли

взаимное

смещение центров

распределения

ошибок

измерений,

т. е.

остается

ли постоянной

систематическая

ошибка

измерений

внутри цикла и между циклами.

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

ивычисленных по ним величин деформаций инженерного сооружения,

атакже разработать практические рекомендации по методике измере-

176

ний, направленные на ослабление влияния систематических ошибок

в последующих циклах

наблюдений.

Решение этих задач целесообразно выполнять на основе стати­

стического анализа рядов разностей

« И Л Ь - ( / . К

(317)

двойных измерений (/q^

и (/0)2- Такой выбор исходной информации

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

Исходя из решаемых на основе статистического анализа задач, целесообразно использовать ряды разностей (317) таких пар измере­

ний, которые отличаются большей независимостью, определены

с

некоторым

разрывом во времени, при разных установках прибора

и

визирных

целей.

 

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

следующую

совокупность статистических критериев.

Проверка нормальности

распределения

 

 

 

 

d

1. Критерий

максимальной

разности:

максимальная

разность

должна

подчиняться

условию

 

 

 

 

 

m rd x \d \^ tm d,

 

 

 

 

 

 

 

(318)

где / = 3

при « > 2 0 ;

г = 4 при

1 0 < « < 2 0 ;

t = 5 при « < 1 0 .

 

 

 

2. Критерий Шарлье: число / разностей d, превышающих предел

 

max|^|>Z«?d,

 

 

 

 

 

 

 

(319)

не

должно

превышать: /^1

при «^50;

1^2 при 30<«<50;

/^3

при

1 0

< «< 3 0 ,

/^ 4 при « < 10. Коэффициент Z выбирается из специальных

таблиц

(Смирнов

Н. В.,

Белугин Д. А. Теория вероятности

и мате­

матическая

статистика

в

приложении

к

геодезии.— М.,

1969)

по

вычисленной вероятности

 

 

 

 

 

 

 

Ф (г )= п - \/2 п .

 

 

 

 

 

 

(320)

 

Указанные два критерия (318) и (319) позволяют выявить грубые

измерения.

 

 

 

 

 

 

 

 

 

 

 

3. Характеристика

эксцесса

распределения

 

 

 

/ =

 

0,8

 

1,8

 

 

 

 

(321)

 

' -

 

 

 

 

 

 

 

 

 

 

 

ч / ^ Т ’

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где величина

ш

П:

(322)

 

 

 

есть вероятная разность i-й группы (ряда). Генеральная совокупность разностей d одного цикла названа здесь рядом, а выборка из этой совокупности, составленная по признаку постоянства условий измере­

ний,— группой.

Число

групп

ряда

должно

быть не

менее

трех,

а объем каждой группы целесообразно

иметь более

1

0 .

 

 

4. Критерий

Романовского

 

 

 

 

 

 

 

 

 

 

 

Л( = 0,707 |х?-1|<3,

 

 

 

 

 

 

 

 

 

 

(323)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xiJiL!Z 2 );

 

 

 

 

 

 

 

 

 

 

 

 

(324)

 

я,-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

г,— число

положительных разностей

в

i-й

группе

(ряде).

 

 

5. Критерий

Ястремского

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,707|Ех?-5|

 

 

 

 

 

 

 

 

 

 

 

 

 

/ = -----= =-- <3,

 

 

 

 

 

 

 

 

 

 

 

(325)

где S — число

 

групп объема п( в ряде

N-ro

объема.

 

 

 

 

Критерии (323), (325) характеризуют асимметрию распределения.

Проверка

равноточности измерений

 

 

 

 

 

 

 

 

 

6 . Критерий

Бартлета

 

 

 

 

 

 

 

 

 

 

 

Хб = £ & < 5 с!.

 

 

 

 

 

 

 

 

 

 

 

 

(326)

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ft= - (n ,- l)In F w;

 

 

 

 

 

 

 

 

 

 

 

(327)

F,„ =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(328)

mdN

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т<и— средняя

квадратическая

разность

i-й группы

данного

ряда;

тdN

средняя

квадратическая

разность

ряда;

%1— критическая об­

ласть, выбираемая

из

таблиц

по

числу

степеней

свободы

K = S — l

и заданному

уровню

значимости

q.

 

 

 

 

 

 

 

 

 

В случае, если отсутствуют таблицы

натуральных

логарифмов,

но имеются таблицы десятичных логарифмов, значения

можно

вычислить по

формуле

 

 

 

 

 

 

 

 

 

 

 

lnF,.N = 2,301g/^.

 

 

 

 

 

 

 

 

 

 

 

(329)

7. Критерий

Кочрена

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'Lmii

Критическая область Gg выбирается из таблиц по числу степеней

свободы

K l =ni

и

K2 = S —1

и заданному уровню значимости q.

Критерий

(330)

применяется

при одинаковом объеме групп ряда,

т. е. nl =n2 = ... = ni

и nS=N .

 

8 . F-критерий

 

 

 

 

 

 

 

(331)

где i и j — номера групп (или рядов) объемами nt, nj (Ni, Nj) соответственно; Fq— критическая область, выбираемая из таблиц по числу степеней свободы ^ = ^ — 1 и Kj = nj—1 и заданному уровню значимости q.

Проверка несмещенности центров распределения разностей d

9. Критерий Снедокора — Фишера

 

 

 

 

 

 

 

(332)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(333)

ср

СР

 

 

 

 

 

(334)

 

 

 

 

 

 

m du — общерядовая

средняя квадратическая разность,

определяемая

по формуле

(250);

mdS— межгрупповая

средняя

квадратическая раз­

ность; 5SJy — разность между средней i

группы

dn.

и

средней

всего

ряда dNc; Fq— критическая область, выбираемая

из

таблиц по

числу

степенейРсвободы

Ks = (S—1) и KN = (N—S).

 

 

 

 

10. Критерий Романовского

 

 

 

 

 

 

 

 

 

 

 

 

(335)

где

 

 

 

 

 

 

 

величина FSN определяется по формуле (332). 1 1 .г(-критерий

где / и j — номера групп (рядов) объемами

nj (Ni, Nj) соответственно;

tq— критическая область,

выбираемая из

таблиц по числу степеней

свободы K=(ni+ nj — '2) и

заданному уровню значимости q.

1 2 . /^-критерий

 

 

SsV4(iV - 2)

(338)

где величина bs определяется по формуле (334); mdN— по выражению (333); i— номер группы объемом п\ N — объем всего ряда; tq— критическая область, выбираемая из таблиц по числу степеней свободы К =(N—2) и заданному уровню обеспеченности q.

Для применения названных критериев нет необходимости в систе­ матизации исходной информации, кроме обычного группирования ее по исследуемому признаку. Для пяти из перечисленных критериев не требуется применять таблицы, для использования остальных критериев необходимы таблицы х 2-, F-, G-распределений, которые приводятся во всех руководствах и справочниках по математической статистике.

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

Все наблюдения выполнены одним исполнителем, одним и тем же прибором за два дня. При этом первый день был жарким, солнечным, второй— прохладным, облачным. Кроме того, условия наблюдений в первый и второй половинах дня также различались между собой. С учетом этого все превышения разобьем на четыре группы: первый день, утро— 1 0 превышений (I); первый день, вечер— 10 превышений (II); второй день, утро— 12 превышений (III); второй день, вечер— 13 превышений (IV).

Рп

Рис. 73. Схема

высокоточ­

ной нивелирной

сети

Разности двойных

измерений

 

 

<i=hпр +^ р .

 

 

 

 

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

группам

(табл. 40).

 

Т а б л и ц а

40. Разности

dx двойных

измерений, мм

 

№№

 

 

 

Группа

 

 

 

 

 

 

 

I

II

 

III

IV

1

—0,12/0,014

-0,20/0,040

10,20/0,040

-0,13/0,017

2

+ 0,12/0,014

+ 0,07/0,005

—0,11/0,012

-0,14/0,020

3

-0,08/0,006

-0,40/0,160

+ 0,34/0,116

-0,22/0,048

4

-0,04/0,002

-0,15/0,022

-0,07/0,005

+ 0,29/0,084

5

+ 0,11/0,012

+ 0,13/0,017

-0,15/0,022

+ 0,31/0,096

6

—0,22/0,048

+ 0,05/0,022

-0,16/0,26

-0,04/0,002

7

-0,04/0,002

-0,27/0,073

-0,18/0,032

+ 0,16/0,026

8

+ 0,01/0

+ 0,17/0,029

+ 0,01/0

+ 0,10/0,010

9

+ 0,08/0,006

-0,33/0,109

-0,27/0,073

+ 0,03/0,001

10

-0,19/0,036

-0,01/0

 

-0,22/0,048

+ 0,10/0,010

11

 

-0,12/0,014

-0,10/0,010

12

 

+ 0,09/0,008

-0,24/0,058

13

 

-0,05/0,002

Примечание. В числителе приведены значения d, в знаменателе

d 2.

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

Т а б л и ц а 41.

Вычисление критериев нормальности

 

 

 

Блок

 

 

 

Критерий

в группе

 

 

 

Стро­

Обозна­

 

 

 

 

(4) +(5) +

Критерий

 

ка

чение

I

II

III

IV

+(6)+ (7)

ряда

 

 

 

 

 

1

2

3

4

5

6

7

8

9

 

1

+ £

+ 0,32

+ 0,42

+ 0,64

+ 0,99

+ 2,37

_

 

2

- I

-0,69

-1,36

-1,28

-0,92

-4,25

 

3

£|</|

1,01

1,78

1,92

1,91

6,62

 

4

Y.d

-0,37

-0,94

-0,64

+ 0,07

-1,88

1

5

п

10

10

12

13

5

 

6

^ср

-0,037

-0,094

-0,053

+ 0,005

-0,042

 

7

+ 0,13

+ 0,16

+ 0,12

+ 0,12

±0,063

 

 

 

8

TLd2

0,140

0,457

0,396

0,384

1,377

 

9

™d

0,118

0,214

0,182

0,172

0,183

 

10

max I d I

0,22

0,40

0,34

0,31

_

0,040

 

 

11

t m d

0,35

0,64

0,55

0,52

0,55

 

12

0 ( Z )

0,450

0,450

0,458

0,461

0,489

 

13

Z

1,7

1,7

1,7

1,8

2,3

 

 

 

 

Критерий

в группе

 

 

 

 

Блок

Стро­

Обозна­

 

 

 

 

 

(4) +(5) +

Критерий

 

ка

чение

I

II

III

IV

+(6)+ (7)

ряда

 

 

 

 

 

1

2

3

4

5

6

7

8

9

 

14

Z m d

0,20

0,36

0,31

0,31

_

0,42

 

 

15

1

1/5

1/5

1/4

1/4

1/3

2

16

vd

0,101

0,178

0,160

0,147

0,147

 

17

0,06

0,03

0,08

0,05

0

 

f

 

18

 

0,60

0,60

0,54

0,52

0,27

 

19

 

4

4

4

6

 

18

 

20

г

0,40

0,40

1,33

0,08

2,21

1,8

 

21

Л,

0,42

0,42

0,23

0,65

 

0,56

 

22

/

0,75

 

23

FiN

0,42

1,37

0,99

0,88

_

_

 

3

24

\n F iN

-0,86

+ 0,31

-0,01

-0,13

 

25

Qi

+ 7,74

-2,79

+ 0,12

+ 1,69

+ 6,76

 

26

 

+ 0,005

-0,052

-0,011

+ 0,047

_

0,138

4

27

&SNn

0

0,027

0,001

0,029

0,057

 

28

0,8

 

RSN

Здесь в блоке

1 выполнены предварительные

вычисления:

суммы

+Z

и — Е

положительных

и отрицательных

разностей d\

сумма

X|d| абсолютных разностей \d\9 сумма 2,d разностей d; среднее значение dhc разностей d, вычисляемое по формуле

*СП= L

( 3 3 9 )

средняя квадратическая разность md, определяемая по формуле Гаусса

если d <Ad, и по формуле Бесселя

тЛ=

если dcp>Ad. Величина Ad вычисляется по формуле

Ad = 3тнуД /^ /п ,

( 3 4 0 )

( 3 4 1 )

где mh— предполагаемая средняя квадратическая ошибка измеренных превышений, образующих разности d (в нашем примере mh = 0 , 1 мм).

Из сравнения строк 6 и 7 заключаем, что во всех группах величину md следует вычислять по формуле (340). В графе 8 приведены

суммы

по

соответствующим строкам,

а в графе 9— значения dNcp

и A dN,

вычисленные по

выражениям

(339), (342) при замене п на

N = 4 5 ,

и

mdv, которая

находится по

формуле

В блоке 2 по пяти критериям проводится проверка нормальности распределения разностей d: наличие в них грубых ошибок (строки

10— 15),

характеристики

эксцесса

(строки

16— 18)

и

асимметрии

(строки

19— 22). В блоке 3 приведены

данные

для

вычисления

критерия

равноточности

разностей

d в группах, а

в блоке 4— для

критериев несмещенности центров их распределения (критерии си­ стематических ошибок).

Перейдем к анализу вычислений, выполненных в блоке 2 табл. 41.

Сравнив данные строк 1 0 и 1 1 , можно

сделать вывод,

что критерий

(318)

в нашем примере выполняется и

для групп ряда (0,22 <0,35;

0,40 <0,64; 0,34 <0,55;

0,31 <0,52; 0,40 <0,55). При этом

принято /= 3,

хотя

в группах «<20

и, следовательно, можно было бы

принять /= 4.

В

строке 15 приведено фактическое

(в числителе)

и допустимое

(в знаменателе) число / разностей d, не превышающих критических значений, вычисленных в строке 14. Следовательно, критерий (319) в нашем примере также выполняется.

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

Сравниваем данные строк 17 и 18 и убеждаемся, что эксцессы эмпирических кривых распределений разностей d для групп и для ряда несущественны (критерий (320)). Анализ данных строки 21 (критерий (323)) свидетельствует, что и асимметрия распределения также несущественна. Незначительность асимметрии распределения

кривой

для

ряда подтверждает и критерий (325),

приведенный

в

 

строке

 

2 2 .

 

 

 

 

 

Таким

образом, на основании вычислений, приведенных в блоке

2

,

можно

утверждать, что

распределение разностей

d в группах

и

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

 

 

 

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

измерений

с

использованием

статистических таблиц

распределений

X 2, F, t

(табл. 42).

 

 

Т а б л и ц а 42. Вычисление критериев равноточности измерений и систематических ошибок

 

 

 

Критерий при

уровне

 

 

 

 

значимости

(в %)

 

Блок

Строка

Обозначение

 

 

 

Фактический

 

 

 

10

5

1

критерий

 

 

 

 

1

2

3

4

5

6

7

1

1

х 2

6,3

7,8

и ,з

6,76

 

2

F.J

2,3

3,2

5,4

3,29

 

 

 

 

 

 

 

3

FSN

2,3

2,8

4,3

0,57

2

4

1,7

2,1

2,8

1,18

 

5

1,7

2,0

2,8

1,56

 

UN

Примечание. В блоке 1 приведены критерии равноточности, в блоке

2 — несмещен­

ности центров распределения разностей

двойных

измерений.

 

Прежде чем перейти к разбору выполненного в табл. 42 анализа, напомним об одном известном положении: для того чтобы уменьшить риск отвергнуть правильную и принять неверную гипотезу, критичес­ кие области х 2>Fq, (я выбирают для трех разных уровней значимости <7=10, 5, и 1%. Таким образом эмпирическое значение критерия сравнивается с тремя его теоретическими значениями. Ответ считается положительным, если эмпирический критерий не превышает хотя бы одного из трех теоретических значений. Эта рекомендация вытекает

из

того положения, что если выдвинутая гипотеза

Н 0 принимается

при

данном уровне значимости q0, то противоположная гипотеза

Нк

отвергается с тем же уровнем значимости qK,

т. е. q0 = qK- Но,

если гипотеза Н0 отвергается при уровне значимости q'0, то вероят­

ность

q'K того,

что верна

противоположная

гипотеза, очень мала,

т. е. q'K<q'0.

 

 

 

Из

данных

строки

1 следует, что

эмпирическое значение

Хб (критерий (326)) не превышает своего теоретического предела

только

при уровнях

значимости

q = 5%

и

q= 1% (6,75 <7,8;

6,76 <11,3)*

В

строке 2

вычислен критерий

(331)

для

групп I, II,

имеющих,

как

это видно

из строки 9

(см. табл. 41),

максимальное

и минимальное значения md. Эмпирическое значение этого критерия (F2l = 3,29) практически находится в пределах критической области при тех же уровнях значимости q = 5% и q = \ % (3,29«3,2; 3,29<5,4). Таким образом, на основании критериев (326) и (331) следует признать, что измерения в группах равноточны, т. е. изменения условий измерений между группами не оказали существенного влияния на точность результатов. В противном случае было бы необходимо уравнивание данного цикла наблюдений выполнять с учетом соотношения весов

групп (например P, = U / > « - ( ^ ) ; '>>“ ( ^ ) • / ’‘ " ( ^ ) '

Критерии (335) и (332), приведенные в строках 28 (см. табл. 41) и 3 (табл. 42) соответственно, показывают, что центры распределения разностей d в группах не смещены относительно друг друга (0,57 <2,3 и 0,83). Критерий (337) (строка 4 табл. 42) подтверждает этот вывод, показывая, что максимально отличающиеся друг от друга средние

разности

dcp групп II и IV

(см. строку 6 табл. 41)

принадле­

жат одной и той же

генеральной совокупности

(1,18 < 1,72). То же

самое утверждение справедливо и в отношении

средних

dcp груп­

пы II и всего ряда,

о чем свидетельствует критерий (339) строка

5 табл. 42

(1,56, 1,7).

Таким

образом, заключаем, что

в измере­

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

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

Соответствующий анализ можно проводить для каждого цикла

наблюдений. Затем следует

проверить, равноточны

ли измерения

в двух смежных циклах I,

II и постоянна ли их

систематическая

ошибка. Этот анализ выполняется с использованием критериев (330), (331), (337), (338) относительно средних квадратических отклонений mds циклов I и II и их средних dNcp. Если измерения в двух смежных циклах окажутся неравноточными, то уравнивание результатов послед­ него цикла наблюдений необходимо выполнять с учетом соотношения

 

т 1

весов, приняв, например,

Если средние разности

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

Ad= dcp- d cptr

(344)

доверительный интервал которой

A d - t pM hn< A d < A d + t pM l и,

(345)

185

Соседние файлы в папке учебники