
книги из ГПНТБ / Бендат Дж. Измерение и анализ случайных процессов
.pdf78 Глава 7
ная функция x(g) есть действительная |
п-значная функция g,.t |
||||
где п — целое |
число. Тогда |
|
|
|
|
Р [g<Г S (х (fe)) < g + Agl |
|
nP [X <~ X Ik) ^ X-f- Ал-] |
|
||
|
A g |
“ |
|
A g |
|
|
пР [.V<( a (к) |
:gi x + A-v] |
A.v |
(3.12) |
|
|
Ал: |
|
|
' Ag |
|
|
|
|
|
||
Следовательно, |
при dg/dx Ф 0 |
|
|
|
|
|
p{g) = n p { x )№ |
up (X) |
(3.13) |
||
|
I dg/dx I |
||||
|
|
|
|
При использовании этой формулы1) необходимо заменять перемен ную X в правой части соответствующим значением g.
Пример3.3. Распределение гармонического колебания. Гар моническое колебание постоянной амплитуды X и постоянной частоты /0 можно рассматривать как случайную величину, если начальная фаза 0 = G(£) есть случайная величина. Рассмотрим, в частности, случай фиксированного времени ^0; тогда случайное гармоническое колебание описывается выражением
X (k)=x (Ѳ)=Х sin [2л/0/0 + Ѳ(/?)].
Пусть C(k) имеет равномерную плотность распределения
^ _1(2л)_1 при 0 ^ Ѳ^ 2я,
\ 0 при других Ѳ.
Какова в этом случае плотность распределения р(х) функции
*(*)?
Здесь прямая функция х(Ѳ) однозначна, а обратная функция Ѳ(х) двузначна. Подставив в формулу (3.13) 0 вместо х и х вме сто g, получим
2 р (Ѳ) |
при |
dx |
Р (*) = dx/dQ |
~Ж Ф 0. |
|
где |
|
|
^ Г = Х сос (2л/0/0+ Ѳ)= Х-j/1 —s.nz (2л)0/0 + 0) = у/ X*—хй.
L) При выгоде данной формулы авторы, по-видимому, предполагают, что величина dg/dx одинакова при всех п значениях х, соответствующих данному g. В общем случае это не так, и правильно записать
Р (Xj)
Р ( е > = Е - ^ - ' dg
I ~г— dxit
где Х[ = xt(g) — функции, обратные функции g(x).— Прим. ред.
Математические основы анализа случайных процессов |
79 |
^Таким образом, |
|
р(х) = \ ( я / * " —х*)~1 ПРИ \х \ < х > |
(3.14а) |
\ 0 при 1хI ^ X |
|
[см. формулу (1.24)]. Соответствующая функция распределения гармонического колебания имеет вид
0 |
при X< —X, |
X |
|
Р(х) = I р(Ю^= = 4 “ ( т _ + агсз'п ' т ) |
при — Х < х < X, (3.146) |
- X |
|
1 |
при X > X. |
Пример 3.4. Неравенство Чебышева. Пусть x(k) — произ вольная случайная величина со средним значением р.х, средним значением квадрата Ф'2 и дисперсией о2. Пусть одномерная плот ность распределения, вообще говоря неизвестная, есть р(х). Тогда
“ |
уЛр (x)dx ^ |
г |
х-р (x)dx ^ е2 |
С |
р (x)dx, |
Ф* _ I |
j |
I |
|||
—со |
|
I X I ^ 8 |
■ |
I л I ^ 8 |
|
поскольку подынтегральное выражение неотрицательно и х2 ^ е2
6ВО всей области интегрирования [|х( |
ej. |
Следовательно, |
|
|
Р [ | * ( * ) | > е ] = f |
Р{x)dx < |
-Hi- |
(3.15) |
|
\ x \ ^ |
e |
|
|
|
Заменим теперь x(k) на x{k) — р,^. Тогда Ф2 заменится на а2 и со отношение (3.15) примет вид
|
Р ||х ( А ) - а л | ^ |
е ] < - ^ . |
(3.15а) |
В частности, положив е — сах, найдем |
|
||
|
Р IIX (k) —цх I ^ |
сах] < - і - , |
(3.156) |
что |
эквивалентно |
|
|
* |
Р [ ! * ( * ) — |
> 1 ------ |
(3.15В) |
Любая из форм неравенств (3.15) известна под названием неравен ства Чебышева.
80 |
Глава 3 |
|
3.1.2. |
Совокупность двух случайных величин |
|
Рассмотрим теперь две случайные величины x(k) и y{k), где k — |
||
точки |
соответствующего выборочного пространства. |
Пусть Р(х) |
и Р(у) |
— функции распределения, соответствующие |
случайным |
величинам х{Іг) и у(/г). Совместная функция распределения Р(х, у)
определяется как вероятность, соответствующая подмножеству точек k выборочного пространства, удовлетворяющих одновремен
но неравенствам x(k) ^ х и y(k) < у. |
Это можно записать так: |
|
Р(х, у ) = Р [ к ( к ) ^ х |
и y{k)^y] . |
(3.16) |
Очевидно, совокупность всех точек k удовлетворяет неравенствам
x(k) |
^ оо и yik) ^ оо. Отсюда |
|
|
|
|
|
|
|
|
||||
Далее, очевидно, |
Р(оо, оо) = 1. |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
Р ( - о о , у)=Р(х, —оо) = 0. |
|
|
|
|
(3.17) |
|||||
Предполагая, как и ранее, что случайные величины непрерыв |
|||||||||||||
ны, |
можно . |
определить |
совместную |
плотность |
распределения |
||||||||
р(х, |
у) соотношением1) |
|
|
|
|
|
|
|
|
|
|
||
|
р(х, |
|
Р |х < х (ft) < |
X -f Дх и у<Су (k) |
< |
у + |
Ду) і |
(3.18) |
|||||
|
у) — Пт |
|
|
Дл:Ді/ |
|
|
|
|
|
|
|||
|
|
|
Д.ѵ->) |
|
|
|
|
|
|
|
|
|
|
|
|
|
Ді/->0 |
|
|
|
|
|
|
|
|
|
|
Заметим, |
что |
|
р (х, |
у) |
^ 0, |
|
|
|
|
|
|
(3.19) |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
j |
j р (х, y)dxdy= 1, |
|
|
|
|
(3.20) |
||||
|
|
|
— со |
|
|
|
|
|
|
|
|
|
|
|
|
|
У |
|
|
|
|
|
|
|
|
|
|
|
Р(х, |
у)— |
1 j p ( 6 , ^didrt, |
J - \ d- |
^ j l \ |
= |
p(x, |
у). (3.21) |
|||||
|
|
|
— СО — СО |
|
|
|
|
|
|
|
|
|
|
Одномерные |
плотности |
распределения |
находятся |
по |
формулам |
||||||||
|
|
|
|
|
СЗО |
|
|
|
|
|
|
|
|
|
|
|
р(х) = |
\p(x,y)dy, |
|
|
|
|
|
||||
|
|
|
|
|
— СО |
|
|
|
|
|
|
|
|
|
|
|
|
|
00 |
|
|
|
|
|
|
|
|
|
|
|
р (У) = |
J |
Р (X, |
y)dx. |
|
|
|
|
(3.22) |
||
_______________________ |
|
|
— ОО |
|
|
|
|
|
|
|
|
*) Из того, что x(k) и y(k) порознь непрерывны, еще не следует существо вание р(х, у) в смысле (3.18). Если при любом фиксированном значении x(k) значения y(k) могут сплошь заполнять некоторый интервал, то это на прак тике может служить интуитивным основанием предположения о существова нии р(х, у).— Прим. ред.
Математические основы анализа случайных процессов |
81 |
ІбрСЛИ |
|
р(х, У)=р (х)р(у), |
(3.23) |
то случайные величины х{іг) и у(/г) называют статистически неза висимыми. Отсюда следует, что для статистически независимых случайных величин
Р(х, у)=Р(х)Р(у). |
(3.24) |
Математическое ожидание действительной однозначной не прерывной функции g(x, у) двух случайных величин х(!г) и y(k) определяется выражением1)
№[g(x, у)\ = j j g (х, у)р{х, y)dxdy. |
(3.25) |
Если, например, g(x, у) = (х(/г) — px)(y(k) — pty), где рх и ру —
средние значения случайных величин x(k) и y(k) соответственно, то формула (3.25) определяет ковариацию Сху величин x(k) и y(k)
Сху= М [(х (k)— рх)(у (k) —Ру)\ =
= М [х (k)y (/г)]—М [*(£)] М [у (£)]=
СО |
|
= J J {х— рх){у— ру) Р (х, y)dxdy. |
(3.26) |
Заметим, что Схх — а | |
есть дисперсия случайной величины x{k), |
|
определяемая соотношением (3.11). |
|
|
Между ковариацией и стандартными отклонениями случай |
||
ных величин x(k), y(k) |
существует простое соотношение ^ |
4 |
|
\Сху\ < о хоу.) |
(3.27) |
Таким образом, ковариация случайных величин не превышает произведения их стандартных отклонений. Доказательство при водится ниже в подразд. 3.2.1.
Как следует из приведенного неравенства, значение норми
рованной |
величины |
|
I |
рxy= S « L - , |
(3.28) |
|
а.\9ц |
|
а) Эта формула справедлива и для разрывных функций g(x, у), во всяком вручав если интеграл в правой части (3.25) имеет смысл и
СО
I f g (х> У) I Р (X, y)dxdy < оо. —Прим. ред.
6 -2 2 4 4
82 |
Глава 3 |
называемой коэффициентом корреляции, находится в предела?^ между --1 и +1. Говорят, что случайные величины х(Іг) и y{k) некоррелированы, если их коэффициент корреляции равен нулю. Понятие некоррелированности следует отличать от приведенного ранее понятия независимости случайных величин. Заметим, что если х(Іс) и y(k) независимы, то, как следует из равенства (3.23),
СО |
|
|
м ІХ (k)y (&)] = j* j* xyp (x, |
y)dxdy = |
|
— CO |
|
|
|
CO |
|
\ xp (x)dx |
I yp (,y)dy = M \x (&)|M [у (/г)]. |
(3.29) |
— 00 |
— CO |
|
Поэтому Схи и, следовательно, рГІІ равны нулю, так что независи мые случайные величины всегда некоррелированы. Обратное утверж дение, вообще говоря, неверно: некоррелированные случайные величины не обязательно независимы. Однако в физически важ ных ситуациях, когда рассматриваются две или более случай ных величии с нормальным (гауссовским) законом распределения, некоррелированность означает также и независимость случайных величин. Доказательство приводится ниже в подразд. 3.1.3.
Пример 3.5. Распределение суммы двух случайных величин. Пусть x(k) и y(k) — случайные величины с совместной плот ностью распределения р(х, у). Определим одномерную плотность распределения p{z) суммы случайных величин
z (k)=x (k) + у (li). |
(3.30f |
Если X +■ у — z, то имеем у = z — x. |
Отсюда |
p(x,' y)=p(x, z— x). |
(3.31) |
При любом фиксированном значении г величина х может прини мать значения в пределах от — оо дооо. Следовательно,'
oo |
|
p (z)= \p{x, z—x)dx. |
(3.32) |
— oo
Таким образом, для нахождения одномерной плотности распреде ления суммы случайных величин z{k) достаточно знать совмест ную плотность распределения слагаемых. Если x(k) и y(k) — независимые случайные величины с одномерными плотностями
распределения pt(x) |
и р2(у) |
соответственно то, |
р(х, у) = |
= РЛФгІУ) = Рі (х)Р2 |
z — X) |
и |
* |
|
00 |
|
|
Математические основы анализа случайных процессов |
83 |
іф.І.З. Гауссовское (нормальное) распределение
Говорят, что случайная величина x(k) подчиняется распреде лению Гаусса (нормальному), если ее плотность распределения имеет вид
р (х) = (& / '2 л ) 1е х р w a)- |
(3.34) |
где а и b — действительные постоянные и Ь > 0. Нетрудно убе диться,. что а и b суть соответственно среднее значение и стандарт ное отклонение случайной величины x(k). Действительно,
w
М[х (&)] = I' хр (x)dx = а —рл,
*оо
М[(x(k)—а)2] = J (х—а)2р (x)dx=b2= a x2.
Таким образом, плотность распределения случайной величины, подчиняющейся нормальному закону, есть
Р(*)=(Дс V ^ |
) : ехр |
(х-р*)2 |
|
(3.35a) |
2оі |
|
|||
цПо определению нормальная функция распределения |
имеет вид |
|||
|
X |
|
|
|
Р (x)=(axf/r2л )_1 |
j expj^- |
2аі |
dl. |
(3.356) |
|
|
|
|
С физической точки зрения важность нормального закона рас пределения объясняется, в частности, существованием централь ной предельной теоремы [30, 33], которая утверждает, что сумма большого числа совместно действующих независимых случай ных величин распределена в довольно общем случае по закону,
близкому к нормальному.
Рассмотрим конкретный пример. Пусть хД/г), x2(k), ..., xN(k) представляют собой N взаимно независимых случайных величин
с любыми законами распределения. Пусть далее рг и of |
— сред |
ние значения и дисперсии случайных величин хД/г), t = |
1, 2,... |
Ц.., N. Рассмотрим случайную величину |
|
* ( * ) = 2 аЛ (*)» |
|
і=і |
|
6* |
|
84 |
|
|
|
|
Глава 3 |
|
|
|
|
|
где ог — произвольные фиксированные |
постоянные. Как следуе£ |
|||||||||
из приведенной |
ниже |
формулы |
(3.137),’ |
среднее |
значение р,* |
|||||
и дисперсия |
случайной величины ' x(Ji) |
определяются форму |
||||||||
лами |
|
|
|
|
|
|
|
|
|
|
|А* = М [*(£)] = |
М |
|
|
|
N |
|
|
|
N |
|
2 |
аіхі (/г) |
|
М Г** |
= |
і=і |
|||||
|
|
_ « = і |
_ |
|
/=і |
|
|
|||
|
|
|
|
N |
|
|
|
N |
||
о!= М [(*(£)—(ц)2] = М |
J ] ai(Xi (k) — щ.)5 = £ |
щві |
||||||||
|
|
|
|
|
і = 1 |
|
|
|
|
|
Последнее выражение есть следствие взаимной |
независимости |
|||||||||
величин xt(k) и |
Xj{k) |
при і ф у'1). |
Центральная |
предельная тео |
||||||
рема утверждает, |
что |
при |
достаточно общих |
условиях и при |
||||||
N —- оо распределение случайной |
величины x(k) |
стремится к |
||||||||
нормальному закону |
со |
средним значением рж и дисперсией ст^. |
||||||||
Рассмотрим |
теперь N случайных |
величин |
хг(к), х2{к), ..., |
Xfij(k), вообще говоря коррелированных. Введем следующие обозна чения для соответствующих средних значений, дисперсий и кова риаций:
\Х; (/?)],
а?=М [(х; (А )-^)*],
|
|
S ' |
С и — М ІС*І (70 — !J-i) (Xj |
], |
С „= а ? . |
Говорят, что случайные величины |
имеют |
N-мерный гауссовский |
(нормальный) закон распределения, если соответствующая ^-мер ная плотность распределения определяется выражением
ѳхр |
f |
— - 2- | С | - ^ 2 |
I cij I |
(xi — щ) (*/ — И/) |
р{хъ X2, . . .,Хы)= |
- Ч__________ |
=1__________________=_ |
||
V |
|
(2n)NU I С 11/2 |
||
|
|
|
|
(3.36) |
гдеС — ковариационная матрица с элементами |
Си , определенная |
|||
ниже, I С I— детерминант |
матрицы С |
и |СІ;-| |
— алгебраические; |
') Первая |
из этих формул справедлива в самом общем случае, лишь бы |
|
М \xi\k)} были конечными, |
а вторая 'праведлива для любых попарно некор |
|
релированных |
случайных |
величин.— Прим. ред. |
Математические основы анализа случайных процессов |
85 |
элемента Сг;-. Точнее |
|
Сп С12 . . |
C-IN |
С= С2і С22 . . C2n |
|
_CfjiCN2 .. |
C n n - |
а алгебраическое дополнение | Си | элемента Cij определяется как детерминант (N — 1)-го порядка, полученный вычеркиванием і-й строки и /-го столбца из матрицы С и умножением полученного минора на (— \)i+J.
, Наиболее важная особенность /^-мерного нормального закона
распределения состоит в том, |
что он |
полностью определяется |
|
средними значениями рг и ковариациями Сі}. При N = 1 соотно |
|||
шение (3.36) принимает вид |
) |
|
|
Р(хі) — (°іѴ |
(*і — Ці)2 |
|
|
1 е2ха?Р |
(3.36а) |
||
|
|
|
т. е. сводится к нормальной плотности распределения, определен
ной ранее формулой (3.35а). |
двумерная нормальная |
плотность |
|||||
При |
N = |
2 |
совместная |
||||
распределения |
имеет вид |
|
|
|
|||
Р(х1, *я)= |
|
|
|
|
|
|
|
ехр (—[2(1—р?,)]-> |
,'Хі — Hl \ 2 |
п ' [Хі — Pi \ (Хъ — (Х2 |
|
||||
ц=___!___________ |
|
|
----------------------» |
||||
|
|
|
|
2 7 1 0 ^ 2 / 1—Рі2 |
(3.336) |
||
|
|
|
|
|
|
|
|
где р12 = |
С,2/аі0 2 |
коэффициент корреляции случайных величин |
|||||
Xi(k) и х2(А'). Заметим, что |
при некоррелированных |
величинах |
|||||
Хі(£) и Xo(k), т. |
е. |
при р12 |
= |
0 имеем |
|
||
|
|
|
|
р{хъ |
хг)= р (х 1)р(хг), |
|
откуда следует, что если величины x^k) и х2(£) некоррелнрованы, то они также независимы. Этот результат не остается справедли вым для любого закона распределения.
Можно записать аналогичные формулы и для распределений
более высокого порядка при N — 3, 4, 5, .... |
Нетрудно показать, |
|||
дто при произвольном N и при некоррелированности всех воз |
||||
можных пар |
нормально |
распределенных |
случайных величии |
|
(т. е. при рг;- |
= 0, |
і Ф /) |
эти случайные величины будут и взаим |
|
но независимыми |
в вероятностном смысле: |
|
||
|
р (хг х2, . . . , |
хы)= р (xjp (х2) . . . , p ( x N). |
& |
|
|
Таблица 3.1 |
|
|
Некоторые плотности распределения |
|
||
Виі иаспредслеинп |
|
|
Плотность распределения |
|
Дискретное |
. _ |
|
р (* ) = Л 6 (х — а) + Bö ( х — |
Ь) + • • • |
|
1 2 |
% |
• • • + N Ö { x - n ) , |
1 |
|
--------1L |
где А - г В + ■■• ■+ N — |
Однородное (прямоуголь ное)
1 |
1 |
( (6 — а ) - 1 при а < х ^ |
Ь, |
а |
Ь |
Р (А) = 1 |
|
|
|
1 0 при других X |
ц. |
Гармонического процесса |
|
|Ѵ|7| |
р |
^ _ |
{ |
(~ V X* — |
при \ х \ < Х , |
||||
|
|
||||||||||
|
-X |
X |
|
|
ІО при других X |
|
|||||
Гауссовское (нормальное) |
|
|
» |
|
|
|
|
|
|
|
|
|
|
|
р (X) = |
(а, |
/ Ъ - Г |
1 |
е - |
|
|||
|
|
|
/ f a : |
|
|
|
|
|
|
|
|
Рэлея |
| / |
\ |
_____р ix) = |
1 |
^ |
e |
- v2/2f2 |
при х > 0 , |
^ |
||
|
|||||||||||
|
О |
|
|
|
|
1 |
0 |
при других X |
|
||
Максвелла |
|
|
|
р |
(X) = |
|
|
|
|
|
|
|
/ |
V . |
( |
^ |
f / т |
- |
е-Х 2/2С 2 ПР И |
х > 0 ’ |
|||
|
О |
|
|
|
|
ѵ |
0 при других X |
|
|||
і |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
Усеченное |
. |
Ср ( х ) |
В предположении, |
|
что исходная |
плот- |
|||||
|
|
1 |
\ |
ность распределения рх(х) определена |
|||||||
|
|
на интервале ( — с » , |
оо), усеченное рас- |
||||||||
|
а |
t ~ |
пределение |
имеет |
вид р ( х ) = Срх (х) |
||||||
|
|
|
|
при а < X < Ь; р |
(х) = 0 при другііх х, |
||||||
|
|
|
|
где |
|
|
|
|
|
|
|
|
|
|
|
|
j р (х) dx = С 1 P i (х) dx = 1 |
||||||
|
|
|
|
|
—со |
|
|
|
а |
|
Математические основы анализа случайных процессов |
87 |
|
-Л |
Продолжение табл. |
3.1 |
* |
||
В и д р а сп р е д е л е н и я |
П л о т н о с т ь р а с п р е д е л е н и я |
|
Срезанное
*
В
лДСгМ'
Л
аі~
В предположении, что исходная плот ность распределения рх{х) определена на интервале (—о о , о о ) , срезанное рас
пределение имеет вид
[ р, (X) при а < X < 6 ,
А5 (х— а) при X ~ а,
р (.V) =
Bö (х — Ь) при х = Ь,
О при X < а пли X > Ъ,
где У p ( x ) d x ~ j p1(х) dx + А +- B = I
—со |
о |
Важность /Ѵ-мерного нормального закона распределения в физических проблемах объясняется, в частности, существованием
многомерной центральной предельной теоремы [30, 33J, которая утверждает, что распределение векторной суммы большого числа взаимно независимых N-мерных случайных величин стремится при достаточно общих условиях к іѴ-мерному нормальному заікону.
Ч В табл. 3.1 дана сводка некоторых частных плотностей рас пределения, используемых в теории вероятностен.
3.2. Стационарные случайные процессы
Случайный процесс |
\xk(l)}, — оо < і < оо (называемый так |
же временным рядом1), |
или стохастическим процессом) представ |
ляет собой ансамбль действительных (или комплексно-значных) функций, который задается посредством описания его вероят ностной структуры. Из соображений удобства в последующем под. переменной t будет подразумеваться время. Каждая отдельная функция xk (/), где t — переменная, а индекс k фиксирован, назы вается выборочной функцией. На практике выборочная функция (или какая-либо конечная реализация выборочной функции вре мени) может рассматриваться как результат наблюдения в от-, даіьном опыте. Возможное множество результатов опытов опре деляет выборочное пространство индекса k\ это пространство мо-*)
*) Под временным рядом обычно принято понимать случайный процесс с дискретным временем.— Прим, перге.