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

книги из ГПНТБ / Химмельблау Д. Анализ процессов статистическими методами

.pdf
Скачиваний:
34
Добавлен:
27.10.2023
Размер:
31.74 Mб
Скачать

92

 

 

 

 

 

 

 

Глава

2

 

 

 

 

 

 

 

 

 

 

А = 0,81 м2 с неопределенностью 3sA =

0,045 м2 . Поэтому

 

 

 

 

 

 

sA = ^

= 0,015 м2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

a i «

4 = 0,000225 м4 .

 

 

 

 

 

 

сгдТдфф: Температура

пара,

который

предполагался

насыщенным,

 

определялась

по

его

давлению,

измеряемому

ртутным

 

манометром. Давление поддерживалось вблизи установ­

 

ленного

уровня 33 360 Н/м в

пределах

 

ошибок

± 2 5 мм

 

рт.

ст.,

т. е.

3336

Н/м. Следовательно, колебания тем­

 

пературы

пара не

превышали

+0, 9 °С. Таким

образом,

 

 

 

 

 

 

S 2

T s =

0,09(°C)2 .

 

 

 

 

 

 

 

 

Температура воды измерялась с помощью двух термопар,

 

каждая из которых обладала максимальной

погрешностью

 

± 0 , 3 °С:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тю — ~2 ІХWl -f-

TWi),

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно, так как Л^эфф = Ts

Tw,

 

 

 

 

 

оІтэфф

«

4 т э ф

ф =

0,09 +

0,005 =

0,095

(°C)2 .

 

adT/dt

'• Производная

от температуры

по времени в момент

t, при

 

котором

АГЭфф = 30 °С, определялась

с

помощью

каса­

 

тельной, проведенной на графике зависимости температуры

 

воды от времени. После ряда построений, при которых

 

использовались различные способы проведения плавной

 

кривой по отдельным точкам и касательных к этим кривым,

 

была получена разумная оценка производной dTldt,

рав­

 

ная

1,5

 

град/мин,

с

дисперсией

 

 

 

 

 

 

 

 

 

 

 

 

oh/dt « éridt « 0,015 (град/мин)2 .

 

 

 

Среднее значение общего эффективного коэффициента тепло­

передачи

£7Эфф в

результате

расчета

оказалось

равным

 

 

 

 

 

fj

80-4187.1,5

n/r

тэ

, г

-град,

 

 

 

 

 

 

 

U m

= о,81-30-60

= 3

4 5

Вт/м2

 

 

 

а оценка

дисперсии

Г0 '0 1 2 8

 

 

 

 

 

 

 

 

 

 

 

-2

а

_

/д/ел2

0,000225

+

0,095 , 0,015-] _

 

 

°иэфф

« 8иафф

-

(МО)

L -gö§

г

(o,81)2

-gor +

T ^ F j

-

 

 

 

= (345)2

[2. Ю - 6 + 3,45.10"4

+ 1,05 • 10"4

+ 0,00667] =

 

 

 

 

 

 

 

 

 

 

 

 

=

118500-0,007122 = 844.

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

Оценка

стандартного отклонения С/Эфф

 

ôu

« Su

=

29.

 

 

эфф

эфф

 

Эта

оценка °"иЭфф является

лишь

приближенной, так как

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

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

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

ний. Как отмечалось во введении

к этой главе, в полную ошибку

может

дать вклад некоторый постоянный сдвиг, или смещение,

и этот

вклад нельзя исключить

повторением измерений. Если

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

всегда удается работать с недавно калиброванными

приборами;

по этой причине не следует

забывать о

возможных

смещениях.

2.4.5.

Выборочный

коэффициент

корреляции

Выборочный

коэффициент

корреляции

р х у представляет собой

оценку коэффициента корреляции Рху> который определяется выражением (2.2.14). Выборочная ковариация случайной пере­ менной выражается следующим образом:

^хг = 1 ^ Т 2 V i - * ) (Yt-Y)

m

(2.4.27)

i

 

 

и может быть

также подсчитана

по формуле

s*r =

-n~T (S ntXtYt—L

2 ntXt

2 n t Y t ) , (2.4.27a)

 

i

i

i

94 Глава 2

если учесть, что

 

 

^(ХІХ)

(Yi —

Y)ni=2lniXiYi—nXY.

г

 

г

Следовательно, выборочный

коэффициент корреляции равен

 

i

 

— K p i r < + 1 .

(2.4.28)

Если предложена некоторая эмпирическая модель и выполнен ряд соответствующих экспериментов с целью определения связи

4,

,

-і,о

о

 

Ar

 

Ф и г . 2.4.9. Плотность

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

выборочного коэффи­

 

циента к о р р е л я ц и и PXY-

 

между двумя переменными и в измерениях одной или обеих пере­ менных содержится ошибка, то, как описано в гл. 4 и 5, можно применить регрессионный анализ. С другой стороны, если просто измерять или наблюдать две переменные в некоторой случайной выборке, можно вычислить меру линейной связи между этими переменными, а именно выборочный коэффициент корреляции. При этом несущественно, являются ли эти переменные независимы­ ми или нет. Если по какой-нибудь причине значения одной из переменных, которая в генеральной совокупности является слу­ чайной величиной, ограничены некоторыми пределами или иссле­ дуются некоторые предварительно отобранные значения пере­ менной, то выборочный коэффициент корреляции дает искажен­ ную оценку коэффициента корреляции [13, 14].

Распределение выборочного коэффициента корреляции весьма сложное. Оно симметрично только при р Х у = 0 и резко асимме-

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

вероятности

и выборочная

статистика

95

трично, если абсолютное значение | pXY слишком большое (фиг. 2.4.9). Фишер

I велико, а число п не нашел преобразование

1 1 - Р х у

Z* = Arc th PXY = y l n — -PXY

такое, что величина Z* распределена приблизительно по нормаль­ ному закону при любых значениях pXY и умеренных значениях п.

• ••

а

о

Ф и г . 2.4.10. Диаграммы рассеяния для гипотетических

данных пр и п р и б л и ­

зительно нулевой к о р р е л я ц и и .

 

Выборочный коэффициент корреляции дает оценку pXY; про­ верки, которые производились 1 ) , основывались на предположении о том, что совместное распределение для обеих переменных X и Y

а

*

б

X

Ф и г . 2.4.11. Диаграммы

рассеяния данных

пр и сильной

положительной

 

к о р р е л я ц и и .

 

 

является нормальным. Отклонение от нормального закона может привести к сильно смещенным оценкам и как следствие к ошибоч­ ным заключениям.

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

г ) Различные проверки, которые

можно выполнить для pXY с помощью

таблиц и графиков, обсуждаются в

[15—17].

96 Глава 2

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

чить, что между двумя переменными может существовать

нелинейная

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

использую­

щим в качестве меры этой связи выборочный коэффициент корре­ ляции. Фиг. 2.4.11 иллюстрирует необходимость использования однородных данных для того, чтобы избежать ложной корреляции, которая возникает, если при вычислении выборочного коэффи­ циента корреляции объединяются две неоднородные группы дан­ ных. Наконец, важно помнить, что значительная корреляция еще не доказывает, что между двумя переменными существует

причинная связь.

Пример 2.4.8. Выборочный коэффициент корреляции

Были взяты восемь проб полимера и измерены две его характе­ ристики: 1) скорость седиментации и 2) степень кристалличности. Чему равен выборочный коэффициент корреляции между этими двумя переменными?

Скорость седиментации

15

И

8

8

6

4

3

1

Степень к р и с т а л л и ч н о с т и

8

8

7

5

4

3

2

1

Решение

Выборочный коэффициент корреляции можно подсчитать по формуле (2.4.28), как показано в табл. П.2.4.8.

 

 

 

 

 

 

Таблица

П.2.4.8

Степень

Скорость

 

 

 

 

 

 

кристаллич­

седиментации

г

- X)

( Г ; - Y )

І - X)*

(YT - Y)2 (XT-X)(YT-Y)

ности Х^

Y i

 

 

 

 

 

 

1

1

 

- 4

—6

16

36

24

2

3

 

- 3

- 4

9

16

12

3

4

 

—2

—3

4

9

6

4

6

 

—1

—1

1

1

1

5

8

.

0

1

0

1

0

7

8

 

2

1

4

1

2

8

11

 

3

4

9

16

12

8

15

 

3

8

9

64

24

40

56

 

б

0

52

144

81

X = 5

У = 7

^(Zt-X)(Xt-Y)

81

PXY-

=0,937

V^iXi-X^Wi-YT-

~ V52 - 14 4

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

выборочная

статистика

97

Вообще говоря, значение pXY

0,937

является

высоким;

чтобы провести соответствующую

проверку,

следует

обратиться

к литературе, цитированной выше.

Задачи

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

сительные частоты

с предсказанными вероятностями.

 

а)

Распределение монет по размерам с округлением до 1 см.

б)

Распределение монет по размерам с округлением до 0,5

мм.

в)

Распределение монет по размерам с округлением до 0,01

мм.

Используйте микрометр.

 

г)

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

орлов при однократном подбрасывании

каждой монеты; распределение решек; распределение попаданий

на ребро. Возможны ли

другие

исходы?

д)

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

по

годам выпуска

монет.

е)

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

по последним

цифрам

года выпуска; распре­

деление по первым

цифрам.

 

 

2.2. Постройте график распределения вероятности и распре­ деления накопленной вероятности дискретной случайной вели­ чины, для которой задано распределение накопленной вероят­

ности

 

 

 

 

Р{Х<х}

 

 

 

 

Чему равно значение Р {X = 3}?

 

 

2.3. При условии

что

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

непрерывной случайной

величины

равно

 

 

 

0,

ж < 0 ,

Р{х) =

Р{Х<Сх}=<

—,п

О-^Сх^Сп,

 

 

 

1,

х>п,

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

2.4. Плотность рэлеевского распределения вероятности равна

.

P(r)=^e-ryz°\

r > 0 ,

98

 

 

Глава

2

 

где

а 2

— постоянная.

Найдите рэлеевское распределение накоп­

ленной

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

(г) и постройте оба графика р (г) и

Р (г)

при

нескольких значениях а 2 .

(г) описывает вероятность

того,

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

деления,

лежит

внутри

круга

радиуса г.]

 

 

 

 

 

 

2.5.

По аналогии

с термодинамикой

«энтропию»

для

дискрет­

ного

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

вероятности

можно

определить

как

 

 

 

 

 

 

 

 

 

 

71

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я ( л ) = — 2

 

 

P{xh)\nP(xk),

 

 

 

 

 

 

 

 

 

 

 

 

 

fc=i

 

 

 

 

 

 

 

 

 

 

 

где

P

(xh)

= P

{X

=

xh}.

В

каком

случае

имеет

место

равен­

ство H = 0 и как интерпретировать ответ? Когда имеет место

наибольшая неопределенность и чему тогда равно значение P

(xk)?

Какова тогда энтропия H (п)1 В каком случае энтропия мини­

мальная? Чему тогда равна H (п)?

 

 

 

 

 

 

 

 

 

2.6.

Плотность

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

вероятности

Максвелла

равна.

 

 

 

 

 

р (х) =

-Щ=х2е-х2'2аг

 

 

U

(х).

 

 

 

 

 

 

 

 

 

 

 

 

а3

у

л

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

оо

 

 

 

 

 

 

 

 

 

 

Покажите,

что р

(х) >

0

и что

j " p

(x)

dx

=

1.

U

(х)

единич-

 

 

 

 

 

 

 

 

 

 

— оо

 

 

 

 

 

 

 

 

 

ная

ступенчатая

функция,

 

а

а

— некоторая

постоянная.

2.7. Совместная вероятность X и Y дана в таблице внизу.

Покажите,

что

X

и

Y

независимы.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

3

 

 

4

 

 

 

 

 

 

 

 

 

 

О

 

1

 

 

1

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

24

 

 

12

 

12

 

24

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

1 .

 

1

 

 

 

 

 

 

 

 

 

 

 

 

12

 

 

6

 

6

 

12

 

 

 

 

 

 

 

 

 

 

2

 

1

 

 

1

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

24

 

 

12

 

12

 

24

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.8.

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

вероятности

X

равна

 

 

 

 

 

 

у2пОх

 

 

 

L

 

2 о х

Л

 

 

 

 

 

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

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

вероятности и выборочная

статистика

99

то какова плотность распределения вероятности X и У, р (х (] у) —

— Р (х> У)і

если X

и

У" — независимые случайные величины?

(Термин П

объяснен

в приложении А.)

 

 

2.9. Какой смысл «совместного распределения накопленной вероятности»? Приведите пример дискретного и непрерывного совместного распределения. Нарисуйте график и обозначьте оси.

2.10.Можно ли интегрировать и дифференцировать случайные величины?

2.11.Найдите математическое ожидание каждой из следую­ щих величин (У" — случайная переменная, у — детерминирован­

ная переменная, а — постоянная):

б ) ^ •

°' dfi '

в) / ( У ) = j(6Y-r-5)e-*£tt;

п

г) f(Y) = S iY\.

i=l

2.12.Найдите математическое ожидание зависимой перемен­

ной в

каждом из следующих дифференциальных

уравнений (X

и У — случайные

величины;

а и Ъ — постоянные;

у — детерми­

нированная

величина):

 

 

 

 

. dW .

dY

v ...

 

 

 

 

...

dY .

dY

, . „

,

 

 

 

6 ) - w + a - d x - = b ( ^ - y ) -

 

 

 

2.13.

Найдите

среднее значение

(математическое ожидание);

а)

рэлеевской

случайной

величины (см. задачу

2.4);

б)

случайной величины X , равномерно распределенной в интер­

вале

a

X

b

и

равной

нулю вне его.

 

2.14. Дл я какой

из

следующих

плотностей

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

вероятности

случайная

переменная

X строго стационарна?

а) р (х)=—у=—ехр

 

— -

 

(нормальное распределение);

б) р (х) = е~м ^j—

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

Пуассона).

 

2.15. Найдите среднее значение зависимой переменной в сле­ дующих моделях процессов (случайные величины обозначены про­ писными буквами):

100

 

 

Глава 2

 

 

 

а)

перенос тепла

 

 

 

 

 

б)

массообмен

d x 2 - U

Т(х2)

=

Т20;

 

С (0,

х) =

0,

 

 

 

 

H-a散

дх*

C(t,0)

=

c0,

 

dt

l i m С (t,

х)

= 0.

 

 

 

2.16. Найдите среднее значение и дисперсию случайной вели­ чины X с равномерной плотностью распределения

=

-

при — 2 - < Z < - 2 - ,

р (ж) =

0

в остальной

области.

2.17. При условии

что

дисперсия

случайной величины X

равна

0,75,

найдите дисперсию следующих случайных величин:

а)

5Х;

 

б)

,

 

в)

X +

7;

X

х - з

 

г )

- 2 ~ -

2.18.В каком случае % {X} % {Y} = g {ХУ}?

2.19.Определите для каждого случая, является ли случайная переменная X стационарной (в широком смысле) или нет, и пояс­ ните, почему.

а)

X

(t)

=

cos (at

+

Г)

(Г — случайная величина);

б)

X

(t)

=

A

cos

(dt +

В sin (ùt (А

и 5 — случайные величины);

в)

X

(t)

=

aY

+

bt

(Y

— случайная

величина);

г)

X

=

aY

+

b

(Y

— случайная

величина).

2.20.

Д л я

заданных ниже случайных

переменных Y (t) и соот­

ветствующих плотностей распределения вероятности вычислите

автокорреляционные

функции.

 

 

 

 

Переменные:

 

 

 

 

 

 

 

а) Y

(t) =

Аеш

 

— случайная

величина, і = У

—1);

б) Y (t) = Ai cos

cot + A2

sin (ot

(Au

A2

— независимые слу­

чайные

величины) ;

 

 

 

 

 

 

в) Y

= Ах

+

b

— не

зависящая

от

времени

случайная

величина).

 

 

 

 

 

 

 

 

Плотности

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

 

 

 

 

 

а ) - Р ( а ) = 4"'

- М І ^ Т '

 

 

 

 

 

р(а) = 0,

\ А \ > \ ;

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

вероятности и выборочная

статистика

101

б) р (а) = се~аіІь,

где с — постоянная,

которую необходимо

определить.

2.21.Найдите автоковариацию для случайных переменных задачи 2.20.

2.22.В таблице записаны значения совместной плотности рас­ пределения случайных величин X и Y. Являются ли они незави­ симыми?

X

1

2

3

1

1

1

4

8

8

1

1

1

8

16

16

1

1

1

8

16

16

2.23.При условии что X — случайная величина со средним

значением

ц х , найдите % {X},

% {2Х},

 

% {X

+

1},

% {2Х

+ 1},

% {X2} и

Ш {X

- ііх}.

Заметим,

что

Ш {X2}

Ф

%2

{X}.

Найдите

соответствующие дисперсии, т. е. Var {X},

Ѵаг {2Х}

и т. д.

2.24. Совместная плотность распределения вероятности двух

случайных

величин

X

и Y равна

 

 

 

 

 

 

 

 

 

 

 

р (х,

у)

= X +

у

при

0 ^

X

^

1,

 

 

 

 

 

 

 

 

 

 

 

 

0 < У

< 1 ,

 

 

 

 

 

 

р (х, у) = 0

 

 

в остальной

области.

 

 

Найдите

коэффициент

корреляции

между

X

и

Y.

 

 

2.25.

Используя

указанные

совместные

плотности

распределе­

ния, вычислите взаимную корреляционную функцию гХу

ІЧі

случайных

функций

X

и

Y.

 

 

 

 

 

 

 

у

 

а) X

и

Y — независимые

случайные

функции,

равномерно

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

соответственно

в интервалах

(0, а)

и (0, Ъ).

б) X

и

Y — случайные

функции

с

совместной

нормальной

плотностью

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

 

 

 

 

 

 

 

 

 

 

р

(х,

у) = А ехр [— (ах2

+ Ьху + су2 + dx +

еу)],

 

где ах2

+ Ьху + су2 + dx -J-

еу^>0.

 

 

 

 

 

 

 

2.26.

Некоторый случайный процесс с сосредоточенными

пара­

метрами

можно

описать следующей

моделью:

 

 

 

 

 

 

 

 

• £ Ж + 2У(*) =

Х(«),

 

 

 

 

 

Y(0) = 0,

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