Файл: Химмельблау Д. Анализ процессов статистическими методами.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 09.04.2024

Просмотров: 602

Скачиваний: 2

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

64

Глава 2

r ( 4 - ) = J £ ,

интеграл в соотношении (г) равен

так что

Если получено достаточно много экспериментальных данных, то, прежде чем предполагать, что они описываются нормальным распределением вероятности, желательно: 1) исследовать распреде­ ление их относительных частот, используя критерии согласия, как описано в разд. 3.7.7, 2) построить график накопленной суммы частот на нормальной вероятностной бумаге которая линеари­ зует график Р (х) благодаря использованию специальной шкалы, или 3) провести какие-либо другие проверки, описанные в гл. 3. Во многих случаях экспериментальные данные хорошо описывают­ ся нормальным распределением вероятности, однако часто это распределение удобно использовать для непрерывных случайных величин, не распределенных по нормальному закону, поскольку

1. Величину можно преобразовать таким образом, чтобы пре­

образованная

величина

имела нормальное распределение.

2. Распределение суммы

случайных величин, распределенных

не по нормальному закону, приближается к нормальному

распре­

делению при стремлении

к

бесконечности

объема

выборки.

3. Ошибка,

связанная

с

применением

статистических

крите­

риев, основанных на предположении о нормальном

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

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

Пример 2.3.3. Графическая проверка нормальности распределения экспериментальных данных2 )

В табл. П.2.3.3 приведены диаметры 200 частиц, задержан­ ных фильтром нефтепровода. Количество частиц, попавших в выбранные интервалы диаметров, определяет частоту появления частиц в каждом интервале. Группировка данных сглаживает случайные колебания, свойственные небольшому количеству дан-

х ) 'Применение специальным образом разграфленной бумаги, линеари ­

зующей

графики

нормального

и д р у г и х

распределений, рассматривается

в книге

[3] и

статье [4].

 

 

2 ) Данные

и

графики этого

примера

взяты из работы [5].


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

вероятности

и выборочная

статистика

65

 

 

 

 

 

 

 

Таблица

Л.2.3.3

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

200 частиц, задержанных фильтром

нефтепровода,

 

 

 

по их

размерам

 

 

 

 

Диаметр, мкм

Частота

Накоп­

Скачки

Среднее

m

Проценты

появления

ленная

значение

п + 1

 

частота

 

скачка m

 

Меньше Ѳ,30

 

2

2

1 - 2

 

1,5

1,5/200

0,75

0,31—0,40

 

33

35

3 - 3 5

 

19

19/200

9,50

0,41 - 0,5 0

 

67

102

36—102

 

69

69/200

34,50

0,51—1,00

 

5

107

103—107

 

105

105/200

52,50

1,01-2,00

 

63

170

108—170

 

139

139/200

69,50

2,01—4,00

 

5

175

171—175

 

173

173/200

86,50

4,01 - 6,0 0

 

11

186

176-186

 

181

181/200

90,50

6,01—8,01

 

1

187

187 •

 

187

187/200

93,50

8,01—10,00

 

11

198

188-198

 

193

193/200

96,50

10,01—20,00

 

1

199

199

 

199

199/200

99,50

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

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

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

частота. Каждому

наблюдаемому

значению

х приписывается

скачок,

или ранг,

минимальное значение которого

равно 1. Если частота

появления

некоторого значения х больше 1, то каждому

значению

приписы­

вается определенный ранг

(т. е. трем наблюдениям приписывается

три последовательных скачка). Дл я каждого

значения

х вычис­

ляется среднее

значение

ранга

по формуле г )

 

 

 

 

 

"V

рангов

 

 

 

 

то =

— = -

^

.

 

 

 

 

 

наблюдаемая

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

 

 

г ) К а к видно

из табл. П . 2 . 3 . 3, среднее значение ранга равно

полусумме

его г р а н и ц . — Прим.

ред.

 

 

 

 

 


66

 

Глава

2

 

По

оси ординат

при построении

графика откладывается

величинах )

 

 

 

 

Р = —^—г

(п — полный объем выборки),

которая

представляет

собой зависимую

переменную.

Так

как график накопленной

вероятности определяет вероят­

ность того, что переменная «равна или меньше, чем», то значение зависимой переменной следует откладывать над верхней границей каждого интервала. График, построенный по данным табл. П . 2 . 3 . 3 , дает сильно асимметричное распределение, на что указывает большая кривизна кривой на фиг. П . 2 . 3 . 3а . Эта кривая, а также физические соображения, что все измеряемые диаметры должны быть положительными и достигать н у л я (размеры частиц не могут быть отрицательными), наводят на мысль построить график

на

логарифмически-нормальной

бумаге, что и сделано

на фиг.

П . 2 . 3 . 36 .

 

 

 

График на фиг. П . 2 . 3 . 36 при

значениях диаметра,

больших

0,5

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

0,5

мкм он претерпевает внезапный излом и затем резко

спадает

до

нуля . Обычно это указывает

на существование каких-либо

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

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

5% проверяемых частиц имели диаметр меньше 0

,1

мкм,

что

можно было

предсказать,

продолжая

сплошную

 

линию

на

фиг. П . 2 . 3 . 36

вплоть до 0,1

мкм. Можно

ожидать,

что частиц

с диаметром больше 50 мкм будет очень

мало (менее

одной

на

тысячу).

 

 

 

 

 

 

Многомерная плотность нормального распределения вероят­ ности (2.3.7) является простым обобщением одномерной плотности

распределения. Д л я компактности запишем ее в матричной форме (приложение Б) 2 )

 

 

р (х )

=

ке-ч^,

(2.3.7)

где положительное

число

 

 

 

 

,

(det f - i ) 1

/ 2

_

1

Согласно

работе

[6].

 

 

 

2 ) Оставшуюся часть разд . 2.3.1, если это необходимо, можно изучать

вместе с гл . 5.

где даны матричные

обозначения .


99,9

О

5

Ю

15

20

 

Диаметр

частиц,

мкм

 

Ф и г . П . 2 . 3 . 36 . График данных табл .

П . 2 . 3 . 3 о размерах частиц на л о г а ­

рифмически - нормальной

вероятностной бумаге .

68

Глава 2

определяемое из условия

называется нормировочным множителем, а

 

Gin

f = Соѵ {Хг Х,} =

— матрица размерности п х п,

l _ Ö n l

. . . Опп J

I f I — детерминант матрицы f.

Величины [X; и Оц — некоторые постоянные, \it является средним значением соответствующей переменной Xt, a сг^- — дисперсия или ковариация величин (Xt, Xj). (Заметим, что tfu =

Пример 2.3.4. Двумерная плотность нормального распределения

Запишем формулу для двумерной плотности нормального рас­ пределения вероятности, полагая в выражении (2.3.7) п = 2. Это распределение применяется при рассмотрении турбулентных полей скоростей, картографировании и плановых съемках, а также при построении эмпирических моделей.

(x — n f =

[(xi — ц 4 ) , (x2

| i 2 ) l i

 

* '—Г о н

cr1 2 " »

 

 

_°*21

G 22.

 

det î = I f Icriiez —a1 2 o-2 i = cr*o\— o2i2

(так как a 1 2 = a2 1 ),

p o2 2

— o"ia -î

 

 

d e t f

d e t f

 

 

Q=[(Xl V-l), (^2-^-1*2)1

z z £ s l

_ £ L

J

l > 2 — И г )

L

| f |

|f|

 

(Ді Hi) 2

°1 2 (жі — |ii) (ж2 — (i2 ) ai2 + (s2 Иг)2 в?


Распределения вероятности и выборочная статистика 69

Теперь пусть

так что

Р = Рі2

=

O-jO-jj

'

 

 

 

СХ12 =

0х02р,

 

 

 

 

 

 

 

 

а\а\ -

of2

-

ojajps

а\ (1 _ рЯ) '

 

 

 

ff^Kl-p2)

 

ffi<r2(l-pa)

'

 

 

 

о?

 

1

 

 

 

Тогда

I f I

° i ( l - P 2

) '

 

 

 

 

 

 

 

 

 

Ѵ~^—) ~2 р

1~5Г~) І -

^ - ) +

\ - 0

— )

 

 

 

1—р2

 

 

 

где

 

 

 

 

 

 

 

 

 

 

і = 1,

2.

 

 

Существуют таблицы

двумерного

нормального

распределе­

ния [7].

 

 

 

когда величины Хі

и Х2 стати­

Д л я важного частного случая,

стически независимы,

р =

0 и

 

 

 

 

 

 

 

2.3.2. Распределение

%2

Распределение %2 находит широкое теоретическое и практиче­

ское

применение;

некоторые

примеры

его применения будут

описаны в гл. 3

*), а

именно:

 

1)

проверка согласия экспериментальных наблюдений с пред­

полагаемыми распределениями

вероятности;

2)

получение доверительных интервалов для дисперсии и сред­

него

квадратического

отклонения;

 

3)проверка независимости переменных;

4)получение выборочного распределения для среднего квадра­

тического отклонения,

ковариации, относительного

отклонения

и т. д.

 

 

Пусть Хи Х2, . . ., Хѵ

— набор ѵ случайных величин, распре­

деленных по нормальному закону с соответствующими

параметра-

*) Методы использования ^ - р а с п р е д е л е н и я хорошо изложены в книге [8] .