Файл: Химмельблау Д. Анализ процессов статистическими методами.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] .