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

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

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

22 Глава 1

в) Дисперсия в трубе с насадкой

d?z .

dz . ,

, ,

1.3. Классифицируйте

каждое

уравнение задачи 1.1 (или 1.2)

по признакам, указанным в табл. 1.1.2.

1.4.Можно ли, выполняя измерения какой-либо величины, выяснить, является ли она случайной или детерминированной?

1.5.В табл. 3.1.5 представлены данные [6] о теплоте парооб­ разования пропана. Имеет ли какое-либо отношение среднее откло­ нение к точности или воспроизводимости представленных данных?

 

 

 

 

Таблица 3.1.5

 

 

Диапазон температур

Среднее откло­

Автор

Номер ссылки

минимальное

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

нение темпера­

по работе [6]

туры парообра­

 

 

значение

значение

зования по всем

 

 

 

 

данным

А

14

100

135

1,12

В

16

103

167

1,43

С

4

100

190

0,98

1.6. Укажите два способа, с помощью которых

детерминиро­

ванный

входной сигнал

х (t). = a cos (dt можно

преобразовать

вслучайный.

1.7.Каким еще путем, кроме способов, указанных на фиг. 1.2.3, можно ввести ошибку в модель процесса?

1.8.Является ли случайной (стохастической) ошибка, возникаю­ щая при численном интегрировании модели в форме дифферен­ циального уравнения? Является ли случайной вводимая в модель процесса.ошибка, связанная с отбрасыванием членов при аппрок­ симации дифференциального уравнения уравнением в конечных разностях?

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

1.10.Можно ли из данного графика относительной плотности распределения измерений некоторой предполагаемой случайной

величины узнать, имеется ли смещение измерений? Поясните.

i l

• •

г

m •

 

J I i _ J L _ j I I I

Введение

23

 

1.11. Является ли случайной величиной функция случайной

переменной?

Объясните.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Л И Т Е Р А Т У Р А

 

 

 

 

 

1.

Himmelblau

D . M . , Bischoff

К . В . , Process Analysis and

Simulation,

W i l e v ,

 

N . Y . ,

 

1968.

 

Technometrics,

 

 

 

 

 

 

 

 

 

2.

Youden

W .

J.,

4,

111

(1962).

 

 

 

 

3.

Youden

W .

J.,

Physics

Today,

14,

 

32

(Sept.

1961).

 

 

4. Eisenhart С ,

J.

Res.

Nat.

Bur.

Standards,

67C, 161

(1963).

 

5. Stanton B . D . , ISA

J., p.

77

(Nov.

1964).

 

 

 

 

 

6.

Helgeson N .

L . ,

Sage

В . H . , / .

Chem.

Eng.

Data,

12,

47 (1967).

 

 

 

 

 

 

 

 

ДОПОЛНИТЕЛЬНАЯ

ЛИТЕРАТУРА

 

 

 

Formby

J.,

A n

Introduction

to

the

Mathematical Formulation of Self-

organizing

Systems,

Van Nostrand, Princeton, N . J . ,

1965.

 

 

Papoulis A . , Probability,

Random

Variables

and

Stochastic Processes,

M c G r a w - H i l l ,

N . Y . ,

1965.

 

 

 

 

 

 

 

 

 

 

 

Peterson

E .

L . ,

Statistical

Analysis

and Optimization

of Systems,

W i l e y ,

N . Y . , 1961, Chs. 1—3;

есть русский перевод: Питерсон

И . Л . , Статистический

а н а л и з

и

оптимизация

систем автоматического

у п р а в л е н и я , изд-во

«Сов.

Радио»,

1964.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Глава 2

Р А С П Р Е Д Е Л Е Н ИЯ ВЕРОЯТНОСТИ И ВЫБОРОЧНАЯ СТАТИСТИКА

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

Эксперимент

Случайный исход Таблица экспериментальных р е з у л ь ­

татов Все возможные исходы

Асимптотическая относительная частота п о я в л е н и я некоторого результата (за большое время)

Таблица асимптотических относи­ тельных частот каждого исхода

На к о п л е н н а я сумма относительных частот

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

Случайная переменная Выборочное пространство

Гене р а л ь н а я совокупность Вероятность события

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

Распределение накопленной веро­ ятности (функция распределе ­ н и я вероятности)

Г

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

 

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

 

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

связать

 

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

 

цией —распределением вероятности, которое достаточно

хорошо

 

соответствует относительной частоте появления этих данных.

 

Тогда,

используя эту

функцию, он

может

делать

различные

 

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

 

эксперимента. Однако часто объем экспериментальных данных

 

бывает

недостаточен и

тогда экспериментатор

в

лучшем

случае

V

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

!

и, быть

может, дисперсии по

ансамблю случайной

величины.

 

В этой главе будут описаны некоторые наиболее часто исполь­

 

зуемые

функции распределения. Кроме того, рассматриваются

 

такие характеристики ансамбля, как среднее значение, диспер­

 

сия, ковариация, коэффициент корреляции, которые используются

 

при анализе процессов. Затем будет обсужден первый из двух

 

основных методов оценивания

средних

значений

по-ансамблю,

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

вероятности

и выборочная

статистика

25

а именно: 1) выборочных средних и 2) средних по времени. В изло­ жение включены некоторые выборочные распределения, которые окажутся полезными при последующих обсуждениях интерваль­ ных оценок и методов проверки гипотез. Средние по времени будут рассмотрены в гл. 12.

2 . 1 . П Л О Т Н О С Т Ь Р А С П Р Е Д Е Л Е Н И Я

В Е Р О Я Т Н О С Т И

И Р А С П Р Е Д Е Л Е Н И Е Н А К О П Л Е Н Н О Й В Е Р О Я Т Н О С Т И

Д л я упрощения записи распределение

накопленной

вероятности

для X (t) будет обозначаться

символом 1 )

 

 

Р{Х {t)

^х}

= Р

(х;

t),

(2.1.1)

где X — некоторое число.

Аргумент

в левой части

тождества

(2.1.1) означает: «все значения

случайной

величины X

(t) меньше

 

Время

t

Ф и г . 2.1.1.

Повторные измерения температуры жидкости в заданной точке;

хі, х2

и х3 — различные уровни

случайной переменной X (t).

детерминированной величины х или равны ей». Использование символа X, а не постоянной к связано с тем, что во многих случаях ограничивающая величина сама является детерминированной переменной. Величина Р (х; t) иногда называется распределением накопленной вероятности первого порядка, потому что это распре­ деление .содержит только одну случайную величину в некоторый момент времени.

Используя понятие частоты, можно дать физическую интерпре­ тацию Р (x;t). Пусть многократно проводились измерения темпе­ ратуры жидкости. В результате было получено семейство кривых X (t), некоторые из которых приведены на фиг. 2.1.1. Д л я каждой

1 ) В отечественной литературе под функцией распределения

и л и интег ­

ральной

функцией

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

принято

понимать

величину,

удовлет ­

в о р я ю щ у ю условию

Р {X (t) <

х} =

Р

(х; t).

 

 

 

 

 

 

 

 

Е с л и переменная X

(t) непрерывна,

это

определение

совпадает

с

(2.1.1).—

Прим..

ред.

 

 

 

 

 

 

 

26

Глава 2

 

 

кривой в момент времени t =

ty проверяется, выполнено ли условие

X (ti) ^ X. Пусть в момент

времени ty будет

всего пн показаний,

для которых X (ty ) ^ X, а полное число кривых равно N.

В пре­

 

дельном случае,

когда

 

N->

оо,

 

 

 

 

 

Ясно,

что

 

величина

 

 

 

 

Р {х; t) лежит в интер­

 

 

 

 

вале между 0 и 1. На

 

 

 

 

фиг.

2.1.2, а

показано

 

 

 

 

некоторое

распределе­

 

 

 

 

ние

накопленной

веро­

 

 

 

 

ятности

для

случая,

 

 

 

 

когда

оно

 

является

 

 

 

 

функцией времени. При­

 

 

 

 

мер

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

на­

 

 

 

 

копленной

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

 

 

 

 

не зависящего от време­

 

 

 

 

ни,

дан на фиг. 2.2.2, б.

 

 

 

 

 

Предположим, что те­

 

 

 

 

перь исследуются экспе­

 

 

 

 

риментальные

графики

 

 

 

 

как

в

момент

времени

 

 

 

 

t

=

ty,

так

и

в

момент

 

 

 

 

t

— t2.

Тогда

совместное

 

 

 

 

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

случай­

 

 

 

 

ных

величин

X

(ty) и

 

 

 

 

X (t2) можно обозначить

 

 

 

 

как

 

 

 

 

 

 

 

 

 

 

 

 

 

Р

(Ху, Х2,

ty,

t2) SI

Ф и г. 2.1.2. Распределения накопленной ве­

 

=

Р {X

(ty) <

Ху,

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

роятности .

 

 

X(t2)^x2};

 

 

 

 

(2.1.2)

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

 

 

 

 

 

ющееся некоторой функцией

времени;

б — распре­

величина

P

(ху,

хг;

ty,

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

зависящее от

 

времени.

 

t2) называется распре­

 

 

 

 

 

 

 

 

делением

 

накопленной

вероятности

второго

порядка величины X (t),

a

Ху и

х2

пред­

ставляют собой некоторые два числа. Слова «второго порядка» отражают тот факт, что изучается совместное распределение одной

и той же величины,

измеряемой в два разных момента

времени.

С частотной

точки

зрения величина P (Xl,

х2;

ty,

tz)

является

предельным

значением доли совместных

событий

{X

(ty) ^ Ху)

и {X {t2) ^ х2) в двумерном пространстве при N-+-

оо. Если вели-

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

вероятности

и

выборочная

статистика

27

чина X не изменяется со временем,

функциональная зависимость

от t может быть опущена.

 

 

 

 

Каждому распределению накопленной вероятности соответ­

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

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

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

сле­

дующим образом:

р(х; t ) =

дР{хд1t] ,

(

 

 

 

 

 

Заметим, что плотность распределения вероятности обозначается строчной буквой р, а распределение накопленной вероятности — прописной буквой Р. На фиг. 2.1.3 приведены кривые типичных процессов и соответствующие им плотности распределения вероят­ ности первого порядка. Использование термина «плотность» ста­ новится понятным, если заметить, что для того, чтобы величина Р (х; t) была безразмерной, единицы измерения р (х; t) должны

быть обратны единицам измерения х, т. е.р

(х; t)

является

вероят­

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

(Иног­

да выражение р (х; t) dx используется для обозначения

вероят­

ности того, что величина х лежит в интервале

от х

до х

+

dx.)

До сих пор рассматривались плотность распределения вероят­

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

непрерывных

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

(но не плотностью)

вероятности

для

дискретной

величины

X

(і) является

величина

Р (xh;

t) ==

= Р {X(t) =

Хъ},

а

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

накопленной

вероятности

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

собой

не

интеграл, а сумму

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

Р(х;

t)=.P{X{t)<.xK}=

2

P(xt;

t).

 

 

 

 

 

 

 

 

i =

l

 

 

 

 

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

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

величины может быть записано в виде

X

 

 

 

Р{х;

t)=

j р(х';

t)dx',

 

 

 

 

 

 

 

— 00

 

 

 

 

 

Р(х2; t) — P{xt;

t)=^p{x;

t) dx

= P { ^ < X

(t)*£x2),

 

 

 

 

 

Xi

 

 

 

 

где

штрихом

обозначена

немая

переменная.

Следовательно,

Р{Х

=

х0} = 0,

так

как

интервал

интегрирования равен нулю.

Кроме

того, по

определению,

 

 

 

о о

j p (х; t)dx = l.

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

вероятности

и выборочная

статистика

29

Аналогичные соотношения можно написать для распределения вероятности второго порядка. Соотношение между плотностями распределений первого и второго порядка имеет вид

о о

р ( я і ; ti)= j р(хи

хг; tu h)dx2.

(2.1.4)

 

— о о

 

 

Величина p (x4 ; ti)

называется маргинальной

(безусловной) плот­

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

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

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

вероятности величины X (tt) безотносительно к тому, каково зна­

чение X (£2)-

Совместное распределение накопленной вероятности двух раз­

личных случайных величин X (t)

и

Y

(t)

называется двумерным

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

и может быть

записано

в

виде

 

 

 

Р (х,

y;t)=P

{X (t) ^х;

Y (t) <

y}

=

 

 

 

 

 

 

 

 

 

 

X

y

 

 

 

 

 

 

 

 

 

 

=

j

,J

p (x',

y'; t) dy'

dx'.

(2.1.5)

 

 

 

 

 

— o o

— o o

 

 

 

 

 

 

Два

типичных,

не

зависящих

от

времени

двумерных'["рас­

пределения вероятности показаны на фиг. 2.1.4.

условная

плот­

В

последующих

главах будет

использоваться

ность

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

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

 

Условное

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

копленной вероятности случайной величины Y,

получающееся

при условии, что случайная величина X

равна значению х,

дается

выражением

 

 

 

 

 

 

 

 

 

 

 

 

P (у I X =

х) =

lim P (у I x <

X < x +

Ах),

 

Дх->-0

где вертикальная черта означает «для данного» или «при условии, что». Затем, используя соотношение для верхнего и нижнего пре­ делов непрерывной величины, аналогичное соотношению (А.8) приложения А, получим

n ,

I V

ч

P(x-\-hx,

у)— Р (х, у)

дР(х,

и)/дх

Р (у \Х = х) = hm рТ I а ; р / t =

Я р

к

а 1

'

дзс^о Р(х

+ Ах) — Р(х)

дР(х)Ідх

 

Соответствующая плотность распределения вероятности получает­

ся дифференцированием P

(у \ X = х)

по у:

 

 

 

 

р(у\Х

= х) =

Р

^ .

 

 

 

(2.1.6)

Д л я упрощения записи условную

плотность

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

ч - Ф и г . 2.1.3. Типичные

диаграммы процессов

(слева)

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

им (не з а в и с я щ и е от времени) плотности распределения

вероятности

(справа).

а — синусоидальная волна

(со случайной начальной фазой

W); б — синусоидальная

волна плюс случайный шум; в — узкополосный случайный

шум; г — широкополосный

случайный шум. (Из работы

[1],

стр.

17, 18.)

,

30

Глава

2

ятности принято

записывать

в

виде

 

Р {у \ х) =

р

{у \ Х = х).

Так как

оо

Р(У) = j Р (Х> У) D X >

— о о

асовместная плотность распределения, согласно выражению

(2.1.6),

равна

 

 

 

 

 

 

 

• р (х,

у)

=

р (у

I х)

р

{х),

можно

записать

 

 

 

 

 

 

 

 

о о

 

 

 

 

 

 

Р (У) =

j

Р

Iх)

Р

(х)

d x -

 

 

— о о

 

условия

X

— х нужно условную

Другими словами, для снятия

Ф и г. 2.1.4. Д в а

не з а в и с я щ и х

 

от времени

двумерных

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

а — двумерная

плотность

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

величин;

б - дву­

 

 

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

 

плотность

распределения умножить на

плотность

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

X и проинтегрировать по всем значениям

X.

 

 

 

Обобщая выражение (2.1.5), можно ввести и-мерное распреде­

ление накопленной вероятности. Изучение п различных

случайных

величин Хі,

Х2,

. . ., Хп

эквивалентно

рассмотрению

одного

тг-мерного

случайного вектора

X = (Xt,

Х2,

• • -,

Хп). Одномер­

ные

переменные

называются

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

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

если

1

) Статистическую независимость двух случайных величин можно

и н т е р -

" претировать

следующим образом. Вероятность получения какого - либо

значе ­

н и я

одной

из

величин не д о л ж н а

зависеть от

того, каково

фиксированное

значение другой величины . Известны непрерывные

процессы,

д л я к о т о р ы х

предыдущие

значения

переменных

действительно

в л и я ю т

на следующие

j з н а ч е н и я .

 

 

 

 

 

 

 

 

 

 

 

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

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

статистика

31

для всех возможных значении этих переменных и для всех сов­ местных распределений накопленной вероятности выполняется условие

Р (х) = Р (хі, . . ., хп; h, . . ., tn) = P (ti) . . . P (xn; tn).

(2.1.7a)

Эквивалентное соотношение для плотностей распределения веро­ ятности имеет вид

р (хи . . ., хп; ti, . . ., tn) = p (U) . . . p (xn; tn). (2.1.76)

Аналогичное выражение для независимых дискретных переменных получается обобщением выражения^(А.б) приложения А

P (xik,

. . ., xnh)

=

P {Xi (t) ^

xih,

. . . ; Xn (t) xnk}

=

 

 

 

 

 

 

 

 

 

 

— P (xlh)

• • • P

(^nft)-

Пример 2.1.1. Двумерная функция распределения

 

 

 

Пусть p (х,

у;

ti,

t2) > 0 — двумерная

плотность

распреде­

ления

вероятности

для двух

случайных

величин

X

(t) и Y (t).

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

P (at < X < а2, Ьі < Y < Ъ2) =

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

j

j Р (ж, y; h,

t2)

dy dx,

(a)

 

 

 

 

 

 

 

ai

bi

 

 

 

 

 

P ( — оо < X < X, — oo < Y < y) =

 

 

 

 

 

 

 

 

 

 

X

 

y

 

 

 

 

 

 

 

 

 

 

=

j

 

J p (ж',

y'; t u

t2)

dy' dx'

(6)

 

 

 

 

 

— OO

— o o

 

 

 

 

 

 

 

j

j p (ж, г/; tu

t2) dy dx = 1.

 

 

(в)

 

 

— o o

— o o

 

 

 

 

 

 

 

 

 

Приведем простой численный пример,"предполагая, что плот­

ность распределения p (х, у; tu

t2)

 

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

 

р (х, у; ti, t2)

= е-( З С + у >

 

 

для X

> О, Г > О,

 

 

p (х, у; tu

t2)

= 0

в

остальной области.

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

2

 

 

 

 

 

 

P ( - | - < Z < 2 ; 0 < y < 4 J = j Ç g - * * - »

dy =

 

 

 

 

О

1/2

=

(е-Ѵ2 _ e-s) (i _ e-i) = 0.4Ç2.

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