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

книги из ГПНТБ / Болдырев, В. С. Методы математической статистики в гидрографии и кораблевождении учеб. пособие

.pdf
Скачиваний:
7
Добавлен:
19.10.2023
Размер:
6.33 Mб
Скачать

Таким образом,

данный интервал

с

вероятность» ß

на­

крывает точку,

обозначающую истинное значение параметра

(рис. 12).

 

 

 

 

 

-<г

 

* е

 

I

у///,//////\//7ЛГ/\_________ х

 

о

' а

а

!

 

 

u

 

^

 

 

Рис.12

 

 

 

Если ошибки измерений распределяются по нормальному

закону, то оценка математического

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

будет

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

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

причем

параметрами этого закона будут

т

и Q l .

Пользуясь

 

*

п

 

этим обстоятельством, можно вычислить вероятность того, что оценка математического ожидания не выйдет за задан­ ные пределы:

Р ( \^ х ~ т х \ ^ ^ ) ^ Р 1

(1.67)

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

Для упрощения вычислений при нахождении величины интервала существуют специальные таблицы обратной функции Лапласа

37

В большинстве руководств вычислено

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

 

 

(1.69)

где

(э~

- точечная оценка среднего квадратическо-

го отклонения.

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

( м о ;

Здесь D * - точечная оценка дисперсии.

Если пренебречь членами, в знаменателе которых стоит

п 1 и п 3 , то получим приближение

или

для нормального закона распределения

38

д= < Г 2 ;

Л- 3 6 '*-

Подставив эти величины в (1.70), получим

или приближенно

 

 

 

 

 

 

s

[

d } ^ ( 3

‘' ■

 

Отсюда среднее квадратическое еткленение дисперсии

 

 

б [ Ъ \

 

( I .7 I )

Среднее квадратические еткленение стандарта найдем,

рассматривая функции

 

 

 

 

f

 

£де

X = S

 

Применяя

теоремы о числовых характеристиках, будем

Л етъ

 

 

 

 

 

 

_ ^

г

I

\

_______

 

 

 

 

 

 

2 p 1 S! 2 G f T

2 G t f ?

Так как для нормального закона

7

то

бГ

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

39

 

в і

(1.72)

 

]/2 (п -і)'

Заметим,

что надежность оценки, полученной по (1.72),

не слишком

велика к

при / 7 = 7 она меньше 0,95.

Теперь можно построить доверительный интервал для дисперсии

О Ы - ß r ^

*

(1.73)

 

(1.73) получена из ( І .7 І ) путем

замены <3 на D .

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

закон распределения оценки дисперсии подчиняется ^ 2- распределению.

С учетом особенностей - распределения построен точный метод вычисления доверительного интервала для дисперсии.

Рассмотрим несмещенную оценку дисперсии

 

Е

( * г ^

2

 

 

 

 

 

п ■1

 

 

 

 

Эта

величина подчиняется

закону

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

с гг-1

степенями свободы.

Такому

хе

2

распределе­

-

нию будет подчиняться

и безразмерная

дробь

вида

 

 

{ n - i) D

 

 

 

(1.74)

 

І А =

-

 

 

 

 

D

40

где D

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

 

Из этой дроби может быть найдено выражение для оценки

дисперсии:

 

 

 

 

п - і

(1.75)

 

 

 

Зная закон распределения величины «/• , мы можем по­

строить доверительный интервал и для дисперсии.

Так как

закон

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

несимметричен, условимся

вы­

бирать

доверительный интервал так, чтобы вероятность

невыхода величины с/-

за пределы этого интервала

по обе

стороны его была одинакова (рис.13).

 

f ( V )

 

 

 

Рис.13

 

 

 

Обозначим эти

вероятности

через Pf = cL

и 4 -

/-О І

 

 

 

~2

2

где ot-

- некоторая

вероятность. Для удобства ремения

задачи

обозначим значения аргумента,

соответствующие

этим вероятностям,

2

2

» Доверительный

через

и

интервал - через

 

 

 

 

 

 

По интервалу

 

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

интервал *Jjb . Этот интервал накрывает истинное значе­ ние дисперсии с вероятностью. ^

 

P(J><D<D2 ) = ß ■

(1 .7 6 )

Таким образом, неравенства

 

М п - і) < „

 

 

V,

д

(1.77)

Ъ % Р - > ѵ

.

 

равносильны, а

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

 

Таблица - распределения составлена так, что в каждой строке даются два значения, соответствуйте границам интервала. Входами в эту таблицу являются число степеней свободы и заданная вероятность. При больших /2 (7 ? ;> 4 о ) распределение становится практически симметричным и поэтому оба множителя одинаковы.

Задаваясь вероятностью 0,68, вычислим погрешность в величине стандарта в процентах для различных п,

(табл .1 .2).

 

 

 

 

 

 

Таблица 1.2

п

3

5

ТО

20

30

40

50

 

%*

0,72

0,38

0,25

0,16

0,13

0,11

0,10

4 2

§ 5. Резко выделяющиеся наблюдения

Из закона нормального распределения следует, что ве­

роятность

появления

ошибок измерения, больших 3 3 »равна

0,003.

Вероятность того,

что

из п

результатов измере­

ний не

будет

ни одного,

имеющего погрешность больше 3 3 ѵ

 

Р ~ ( і - 0 90 0 3 ? * * і- 0 9 0

0 3 п .

(1*79)

Т .е. вероятность того, что одно из 10 измерений

будет иметь погрешность > 3 3

,

составляет

0,03, одно

из ТОО измерений - ~

0,3 .

 

 

 

 

Для того

чтобы решить вопрос

о

закономерности по­

явления

в

ряду наблюдений отсчетов,

резко

отличающихся

по величине от соседних, необходимо составить безраз­

мерную дробь

вида

 

 

 

**■'т ах ~

(1.80)

 

 

£/ =

где

Xm a x

<*х

 

- максимальная или минимальная величина

 

 

в данном ряду наблюдений;

 

 

/77. •7 е х

- параметры распределения, вычисленные

 

с учетом "подозрительного"

результата.

 

 

Считая, что данные измерений распределяются по нор­ мальному закону, можно получить закон распределения дроби V- . Это распределение будет зависеть от объема выборки п . По объему выборки п и заданной вероят­ ности Р можно найти такое число , при котором будет соблюдаться неравенство

= Р .

(I .8 I)

43

Таким образом, для решения задачи необходимо: вычис­

лить V- I по

заданной вероятности,

пользуясь таблицей

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

, найти величину

, а затем сравнить

табличное и полученное значение.

Пример. В результате измерений получены следущие

отсчеты:

 

 

 

 

3,68

3,11

4,76

2,76

4,15

5,08

2,95

6,35

3,78

4,49

2,81

4,65

3,27

4,08

4,51

4,43

3,43

3,26

2,48

4,84

т „ ш 3,943f

<3L.= 0,943;

сл =

Ь 35г5г. Ш — = 2,552

А

 

 

*

 

0,943

 

 

Зададимся вероятностью 5#. Из таблицы распределения

V найдем по

Р = 5# и п

= 20

^

= 2,623.

 

Значение

нашей дроби меньше

,

и поэтому

отсчет

6,35 следует оставить, так как с

вероятностью

95# он

не нарушает закона распределения.

 

 

 

 

§ 6

.

Отличие реальных

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

ошибок

 

 

нормального

закона

рас­

 

измерений отпределения

 

 

 

 

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

среднего арифметического, полученного по малому числу измерений, можно применить закон распределения Стьюдента. Это распределение (Т - распределение^выведено для дроби

г=

Х 9

(1.82)

 

где X - центрированная случайная величина, распре­ деленная нормально с параметрами 0 |І ;

А- число степеней свободы ( к = п - і )

Ценным качеством распределения Стьюдента является то, что сюда не входит дисперсия. Эту дробь иногда за­ писывают в виде

 

 

 

 

(I.8S )

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

Стьюдента

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

 

 

 

 

(1.84)

Математическое

ожидание

Т = 0 (/*пг =

0 ).

При увеличении

А - распределение

Стьюдента стре­

мится к нормальному (рис.14).

Рис. 14

4 5

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

Указанные особенности могут быть объяснены наличием в реальных распределениях асимметрии и эксцесса. Так, например, эксцесс погрешности определения поправки ин­

декса

секстана - Оі Г8 ;

 

 

- ошибки измерения горизонтального направления -

О’, 27;

 

 

 

 

-

ошибки измерения высоты Солнца секстаном - 0^26;

-

ошибки в измерении наклонения

горизонта -

0', 64

(для

закона

равномерной плотности

эксцесс равен

- 1 ,2 ).

Для тоге

чтобы учесть влияние

асимметрии, рассмот­

рим обобщение формулы плотности распределения для нор­ мального закона.

Формула нормального закона для центрированной слу­

чайной величины

 

в '<£

(І.8~ч

может быть обобщена в виде

f p x 2

(1.86)

f & h C e ^

Коэффициент С найдем в результате нормировки, так как должно быть справедливо очевидное соотношение

46

Соседние файлы в папке книги из ГПНТБ