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

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

.pdf
Скачиваний:
11
Добавлен:
30.10.2023
Размер:
19.04 Mб
Скачать

9 0

ВСПОМОГАТЕЛЬНЫЕ СВЕДЕНИЯ ИЗ ТЕОРИИ■ВЕРОЯТНОСТЕЙ

[ГЛ. Н

ческой устойчивости результатов эксперимента, лежит в основе ма­ тематической статистики. Число Р называется вероятностью собы­ тия Л.

Очевидно,

что

 

н

так

как,

согласно

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

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

приблизительно

равна

частоте

некоторого

события,

то ясно,

что

 

0 < Р < 1 .

 

 

 

 

 

(2.1)

 

 

 

 

 

 

 

 

Если

А — невозможное событие,

то

Р = 0.

Если

же А — досто­

верное

событие, то

Р = 1 .

С

другой стороны,

если

известно,

что для

некоторого

события

А

мы имеем

Р — 0,

то А не обяза­

тельно будет невозможным событием. Действительно,

равенство Р — 0

означает лишь, что частота события А при большом числе N экспе­

риментов приблизительно равна нулю, так что при большом N

событие

А хотя и произойдет, но лишь

в очень малой доле всех

случаев.

Аналогичное

рассуждение

можно

было бы

провести для

случая Р — 1.

3. Функции распределения вероятности

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

лагать,

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

неко­

торое

значение

случайной

величины

причем вероятность

того,

что £ =

равна

Р.

 

 

 

Предположим,

например,

что наш

эксперимент состоит в

изме­

рении величины ошибки какого-либо прибора. В результате первого

эксперимента мы получим значение ошибки, равное

а в резуль­

тате N - го эксперимента мы получим значение величины ошибки,

равное

 

 

 

 

 

Будем теперь

откладывать вдоль оси абсцисс величину ошибки £,

а вдоль

оси

ординат

вероятность того, что ошибка

не превышает

значения

\ =

х

(эта

вероятность приближенно может

быть опреде­

лена при помощи суммирования частоты всех наблюденных значений ошибки, меньших или равных х). Эту вероятность будем обозначать

Я (5 О ) .

Полученная таким способом кривая F (х) обычно называется кривой интегрального распределения вероятности для величины ошибки. Очевидно, что интегральная функция распределения F (ос) =

= Я ( £ ^ а') есть неубывающая

функция х, обладающая свойствами:

0 < Д ( * ) < 1 ,

F ( — оо) = 0,

F (-)- со) = 1.

Действительно, вероятность

значения ошибки, равной — оо,

равна нулю; вероятность того,

что значение ошибки заключено

4]

ДВА

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

91

в пределах (— оо,

+ о о ),

равна единице, и

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

что

значение ошибки

заключено в пределах (— оо, х),

при возрастании х,

может только

возрастать.

 

 

 

 

Интегральной

функции

распределения вероятностей F (х) можно

дать простую

физическую

интерпретацию,

если

представить

себе

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

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

х,

т. е. находящихся слева от точки х, равнялось F (х).

 

Производная

 

 

dF {х)

w (д:)

(2.3)

dx

 

 

в точке х, если она существует, называется плотностью вероятностей или функцией распределения вероятностей случайной величины X и

может рассматриваться

как

плотность

единичной

массы в рассмат­

риваемой точке.

 

 

 

 

 

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

что

величина

X принимает

значение, принад­

лежащее интервалу д: < X<

х +

для

малых Дх, равна w ( x )d x ,

или

 

 

 

 

 

Р (х < X< х -(- Дд:) = w (л:) dx.

Этот дифференциал называется элементом вероятности данного

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

 

точка

разрыва функции F (л:)

со скачком, рав­

 

Если

х 0

есть

ным р0,

то масса р0 будет

сосредоточена в точке

х0, которая назы­

вается точкой

сосредоточения

массы

в данном распределении. С дру­

гой

стороны,

если

х 0

есть

точка

непрерывности,

то

количество

массы,

находящееся

в

интервале

(д:,

дг +

Дд^),

будет

стремиться

к нулю

при Дд: —>0. Если

существует

производная

F' (х) — w (х),

то

средняя

плотность

стремится

к

w (х)

при

Дд:-*0 и, таким

образом,

w (x)

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

плотность массы

в точке х.

 

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

ция, интеграл который по

бесконечному интервалу (— оо, оо)

равен

единице.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4. Два типа распределений вероятностей

 

 

 

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

бывают

одного из

следующих

двух типов:

диск­

ретного

и непрерывного.

 

 

 

 

 

 

 

 

 

 

 

 

Случайная

величина

X называется

 

величиной

дискретного

типа

или обладающей распределением этого типа, если вся распределенная масса содержится в точках сосредоточения массы и если, кроме того, всякий конечный интервал содержит конечное число таких точек. Обозначим точки сосредоточения массы через дг,. х 2.......... Д'дт, а соответствующие массы через р {, р2. . . . . р,\. Распределение случай­ ной величины X будет полностью определено, если указать для любого

F >Х)

92 ВСПОМОГАТЕЛЬНЫЕ СВЕДЕНИЯ ИЗ ТЕОРИИ ВЕРОЯТНОСТЕЙ [ГЛ. II

v(v =

1, 2.........

N) вероятность р„ того, что £

принимает значение х ч.

Так

как вся

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

то

 

 

N

 

 

 

 

2 a = i .

 

 

 

 

V= 1

 

 

Интегральная функция распределения F (х) задается при этом

соотношением

 

 

 

 

F(x) = Я ( ? < * ) = 2

Р„

(2.4)

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

которых

л г„^ х .

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

F (х) есть

 

ступенчатая

функция,

 

 

Fix)

 

 

----------

 

 

Р*

 

 

 

 

1

 

 

 

 

 

 

 

 

}

г

!

 

 

1

3

 

 

 

|

 

г

а

 

 

 

 

2

1

 

 

 

 

 

,

-----1

 

1

 

|

 

 

I

 

Ps

 

 

1:

1

 

1

 

1 -

. I

. :------ —I------ 1—!

 

xi

х,

хг

о х3

xt x s

хг

 

0 х3

x s

 

 

Рис. 2.2.

 

 

 

Рис. 2.3.

 

 

равная постоянной величине на любом интервале, не содержащем точки л\, и имеющая в каждой точке х ч скачок p v.

Распределение дискретного типа может быть графически пред­ ставлено при помощи ступенчатой функции F (х) (рис. 2.2) или графика, на котором каждой точке х ч соответствует ордината высоты р ч

(рис. 2.3).

Величина ? называется непрерыв­ ной случайной величиной, если функ­ ция распределения F (х) всюду не­ прерывна (рис. 2.4) и плотность вероятности w(x) — F/ (x) сущест­ вует и непрерывна для всех значе­ ний х. В этом случае

 

Л*

F(x) =

P ( k < х ) = J® (S)d5.

 

(2.5)

Это распределение не имеет точек

сосредоточения массы и,

следовательно, вероятность того, что величина \ принимает частное значение равна нулю при любом jcv:

■ / > ( $ = * , ) = 0 .

41

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

9 3

Вероятность того, что величина £ принимает значение, принадле­ жащее конечному интервалу (хх, х 2), равна

P ( x 1 < | < x 2) = F ( ^ 2) — / 7 (x 1) = f w (H )d l

( 2 . 6 )

l

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

СО

Jw $ ) d \ = 1.

СО

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

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

графиком функции

распределения F (х)

(рис.

2.4)

или плотности вероятности w (x )

(рис. 2.5). Кривая

y =

w(x)

на­

зывается

кривой

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

данного

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

Во многих задачах*) удобно

 

 

 

 

 

характеризовать случайную пе-

 

 

 

 

 

ременную не функцией распре­

 

 

 

 

 

деления, а характеристической

 

 

 

 

 

функцией.

 

 

 

 

 

 

 

 

 

Характеристическая

функ­

 

 

 

 

 

ция (у’ш)

случайной перемен­

 

 

 

 

 

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

 

 

 

 

 

 

00

 

 

 

 

 

 

 

W (у'ш) =

J

е~!шх w (х) dx,

 

 

 

 

 

 

 

 

 

(2.7)

 

 

 

 

 

где w (х) — плотность

вероят-

 

Рис- 4.5.

 

 

 

ности.

 

 

 

 

 

 

 

 

 

Отсюда следует, что характеристическая функция может рас­

сматриваться

как

преобразование

Фурье

плотности

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

и, следовательно, w(x) можно определить, зная Ф (у'ш), по обычным формулам преобразования Фурье (1.34):

СО

w (х) =

Исходя из свойств плотности вероятности w(x), легко показать2),* что характеристическая функция имеет следующие свойства:

1)

ЧЧО) =

1 ,

2)

1^(7“) | <

1. — оо < ш < оо,

х) См., например, Райс С., Теория флуктуационных шумов, в сб. «Теория передачи электрических сигналов при наличии помех», ИЛ, 1952.

2) См. работы, цитированные на стр. 88,

9 4

• ВСПОМОГАТЕЛЬНЫЕ СВЕДЕНИЯ ИЗ ТЕОРИИ ВЕРОЯТНОСТЕЙ

[ГЛ. И

3) характеристическая функция Ч; (уш) суммы независимых слу­ чайных переменных х = х 1- \ - х 2-\- . . . - \ - х п равна произведению характеристических функций этих переменных, т. е.

'Р'(у'ш) == Ф, (у'ш) ВДсо) . . . 4J„ (усо),

где vF(y'io) означает характеристическую функцию случайной пере­ менной х.

б. Средние значения и моменты

Определим среднее значение случайной величины с дискретным

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

,

 

Пусть с экспериментом

А связаны "события

£а, . . . , & и пусть

p it р2..........р$ — вероятности этих событий, тогда

 

N

v = 1

Средним значением или математическим ожиданием величины Е называется сумма')

N

(2. 8)

представляющая собой взвешенное среднее значение величины £, причем весами служат вероятности в точках сосредоточения массы.

Точно так же, если g (£) представляет собой функцию случайной величины С то ее средним значением или математическим ожиданием называется сумма

 

N

 

g =

M [g ® ] = '2,p,gCb).

(2.9)

 

V-1

 

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

видим, что для экспериментального

опреде­

ления среднего значения нужно взвесить наблюденные значения случайной величины по их частотам, т. е. взять сумму произведений наблюденных значений £2, . . . . на их частоты.

Рассмотрим теперь способ определения среднего значения слу­ чайной величины \ с непрерывной функцией распределения.

Так как

частота

значений £ в интервале

от

х до х - \ - А х при­

ближенно равна w(x)A x,

а значение величины $

в этом

интервале

приближенно

равно

х, то

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

приближенно

равно

х = М (х) ^ ^ xw (х) Ах,

Ад*

*) Здесь и далее математическое ожидание случайной величины обозна­ чается значком /V,

5] СРЕДНИЕ ЗНАЧЕНИЯ И МОМЕНТЫ 95

где суммирование

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

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

меньшего

x min до

наибольшего

х,пах

значения

величины

х. В пре­

деле при

Дх —>•0

мы

 

получим

точную

формулу

для среднего зна­

чения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Л ш а х

 

 

 

 

 

 

 

х =

М (х) =

 

 

1" xw (х) dx.

 

 

(2.10)

 

 

 

 

 

 

 

 

'v m i n

 

 

 

 

Обычно эту формулу

записывают в виде

 

 

 

 

 

 

х =

М(х')~-

j ' x w ( x ) d x ,

 

 

(2.11)

 

 

 

 

 

 

 

 

— СО

 

 

 

 

не уточняя интервала распределения величины £.

 

 

 

Точно

так

же

среднее

значение

величины .§ ■ ( £ ) , представляющей

собой функцию от %,

определяется

выражением

 

 

 

 

 

 

g =

 

M { g ) =

 

J

g ( x ) w { x ) d x .

 

(2.12)

Учитывая,

что

 

 

 

 

 

— СО

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w (х) dx =

dF (х).

 

 

(2.13)

вместо (2.12)

можно

написать:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

СО

 

 

 

 

 

 

 

g = M { g ) =

 

J

g (х) dF (х).

 

 

(2.14)

 

 

 

 

 

 

 

 

— СО

 

 

 

 

 

Если

для

некоторого

целого

положительного числа v функ-

ция xv интегрируема на (— со,

 

оо)

относительно

F (х), то среднее

значение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cv =

M (xv) =

J x vd/?(x)

 

 

(2.15)

 

 

 

 

 

 

 

 

 

— СО

 

 

 

 

называется моментом

порядка

v.

Очевидно, что

момент

нулевого

порядка

 

 

 

 

 

 

 

СО

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а0 =

f

dF{x)

 

 

(2.16)

 

 

 

 

 

 

— СО

 

 

 

 

 

всегда существует и равен единице.

Из механики известно, что момент первого порядка, выраженный

интегралом

 

т0 = а1— J x d F (х)

(2.17)

96

ВСПОМОГАТЕЛЬНЫЕ СВЕДЕНИЯ ИЗ ТЕОРИИ ВЕРОЯТНОСТЕЙ

[ГЛ. II

и равный среднему значению самой случайной величины, предста­ вляет собой абсциссу центра тяжести массы, имеющей распределе­ ние F (х). Точно так же момент второго порядка

аг — J х 2 dF (х)

(2.18)

— СО

представляет собой момент инерции массы с распределением F (х) относительно перпендикулярной оси, проходящей через точку х = 0.

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

Действительно, дифференцируя (2.7) по to и полагая затем со = 0,

получаем:

00

dnW (уо>) I

 

 

f

 

 

[e_/mv]u>=odx =

 

du>n Ju)=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ( - j y

1' x n w (x) dx

= ( ~ j ) n V

 

откуда

 

 

 

 

dnW (Ja) I

 

 

 

 

 

(-У Г " [

 

 

Если x0— постоянная,

 

rf<on

Jm=0‘

 

то величины

 

 

 

 

 

СО

 

 

 

 

 

 

 

£ [(S — jc0)v1=

f

(x — x oy d F ( x )

(v = 0,

1, 2....... A7) (2.19)

называются моментами относительно точки х0.

 

 

Моменты относительно

центра

тяжести массы, имеющей распре­

деление F(x), т. е.

относительно

точки

х 0 =

т0, называются цен­

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

 

 

 

СО

 

 

 

 

 

 

 

 

 

 

 

 

Р, =

£[(* — /но)’] =

 

f

т0У dF (х).

(2.20)

 

 

 

 

 

— ОО

 

 

 

Раскрывая выражение

(х — т0У,

легко найти соотношение

между

моментами и центральными моментами:

 

 

 

й р =

1 .

 

 

 

 

 

 

 

Р. =

0.

 

 

 

 

 

 

 

Рг — °-2«о.

 

 

 

 

 

(2 .21)

Рз — ct3Ътйа

-|- 2/Ио.

 

 

 

Р4 = a t — 4от0а3 + 6т1а2— З т * .

6] ХАРАКТЕРИСТИКИ РАСПРЕДЕЛЕНИЯ 97

6. Характеристики распределения

На практике часто бывает важно описать распределение хотя бы в общих чертах при помощи немногих простых параметров. Рассмот­

рим здесь некоторые из наиболее

важных параметров такого рода.

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

 

00

 

 

 

 

 

 

Д4 ($) =m0== J х dF (х)

(2.22)

 

 

 

— ОО

 

 

— первый момент

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

есть

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

самой слу­

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

 

 

 

 

 

Среднее значение т0 имеет,

как мы видели, простую

физическую

интерпретацию: это — абсцисса

центра

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

распределе­

ния F(x) единичной массы.

 

 

 

 

Среднее арифметическое значение

а

 

 

 

 

ОО

 

 

 

а = М ( | $ | ) =

f j x j d F ( x )

(2.22а)

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

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

ожидание абсолютного значения

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

 

 

 

 

Срединным (или вероятным) значением Ъ случайной

величины ?

называется ее значение на границах интервала (— Ъ< %<

Ь), вероят­

ность нахождения внутри которого равна 0,5. Итак,

величина b

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

 

 

 

ь

ь

 

 

 

 

j " d F ( x ) =

J w ( x ) d x = 0,5.

(2.23)

 

 

 

 

Среднее значение квадрата случайной величины $

 

 

 

СО

 

 

 

 

М ( ? ) =

J *2 dF(x) = <x2

(2.24)

представляет собой момент второго порядка случайной величины. Отклонением Д; случайной величины называется разность

Д; = S — т0.

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

7 Зак. 1083. В. В. Солодовников

98 ВСПОМОГАТЕЛЬНЫЕ СВЕДЕНИЯ ИЗ ТЕОРИИ ВЕРОЯТНОСТЕЙ [ГЛ. II

В качестве

меры рассеяния

относительно

среднего

значения т0

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

центральный

момент:

 

 

 

00

 

 

^

= ^ [(5 — ОТо)2]= f (х - т0)* dF (х).

(2.25)

 

 

— СО

 

 

Он называется дисперсией случайной величины и представляет момент

инерции распределенной массы

относительно перпендикулярной

оси,

проходящей через ее центр тяжести.

 

 

Если среднее

значение т 0=

0, то дисперия сводится к среднему

от квадрата случайной величины.

 

 

Чтобы иметь

характеристику

рассеяния той же

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

что

и сама случайная

величина £, часто предпочитают

рассматривать

не­

отрицательное значение квадратного корня из (32, которое называется

стандартным или средним квадратическим

отклонением и обо­

значается через с:

 

 

 

c =

=

.

(2.25а)

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

Нормальная интегральная функция распределения имеет вид

F(x) = ^

f e ~ * d l

(2.26)

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

1

W (х) = F' (х) : (2.27)

Y 2л

Графики этих функций приведены на рис. 2.6 и 2.7, Среднее зиа-

W(Xj

Fix}

Рис. 2.6,

Рис. 2.7.

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

7 НОРМАЛЬНОЕ РАСПРЕДЕЛЕНИЕ 99

]

Действительно,

J

J

х dF (х)

2 йх = 0;

(2.28)

х2dF(x)

■f.

(2.29)

2 dx — 1 .

При х = 0 кривая w(x) имеет максимум, равный

w (0 )

1

(2.30)

2^ '

/

 

Случайная величина ij называется нормально распределенной с пара­ метрами т0, с, если ее интегральная функция распределения есть

F ^ х — >я0j и р

 

определяется

формулой

(2.26).

При этом

плот­

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

равна

 

 

 

 

 

 

 

 

 

 

 

„ W

1

,

i

T i

)

.

1

 

(х-то?

 

(2.31)

 

= i

f ' (

=

- j71= e

-

 

В этом случае

согласно

формулам

 

(2.17), (2.21) будем иметь:

 

1

'п

( х - т 0У

 

 

 

 

си

 

 

 

 

 

 

 

 

 

 

 

 

 

2 dx — mQ, (2.32)

I хе

 

dx =

у4=-

/

(«о -1- «

) е

ai — c y i s

 

— СО

 

 

 

'

 

 

— СО

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

СО

 

 

 

 

 

 

 

 

 

 

 

 

с

j

(х — т0)2 е

2^

dx =

4 =

/

х^е

2 dx =

с2, (2.33)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и, следовательно, т0 и

с обозначают, как обычно,

среднее

значение

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

или

среднее

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

отклонение.

 

 

Изменение значения

т0 вызывает только

смещение кривой

те/(х)

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

Чем меньше с, тем большая часть массы сосредоточена в окрестно­ сти точки х = т0.

Из определения среднего арифметического а, среднеквадратиче­ ского с и срединного Ъ значений случайной величины £ [см. фор­ мулы (2.22а), (2.25а), (2.23)] ясно, что задание одной из них доста­ точно для определения двух остальных (если известна функция рас­ пределения).

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