МАРТ_Часть2_Главы4_8
.pdfГлава 4. КЛАССИФИКАЦИЯ СЛУЧАЙНЫХ ПРОЦЕССОВ
Опираясь на определение СВ, можно дать следующее общее понятие случайного процесса (СП) [1]: случайным процессом (t) называется функция двух аргументов (t, ), где , – множество элементарных событий; t T, T – область определения функций (t, ). При фиксированном значении t (t, ) является случайной величиной, а для каждого фиксированного(заданного элементарного события) (t, ) зависит только от t и определяет реализацию СП (траекторию, выборочную функцию). Как и для СВ, область значений (t, ) может быть счетным (в том числе и конечным) и континуальным множеством Ф.
Классификацию СП можно осуществить по характеру множеств Т и Ф и виду статистической связи между значениями СП, соответствующими различным моментам времени t T.
В зависимости от характера множеств Т и Ф случайные процессы можно разделить на четыре класса (примеры их реализаций приведены на рис. 4.1,
1– 4 соответственно):
1)процессы с дискретными состояниями и дискретным временем (дискретные случайные последовательности);
2)процессы с дискретными состояниями и непрерывным временем (дискретные СП);
3)процессы с непрерывными состояниями и дискретным временем (случайные последовательности);
4)процессы с непрерывными состояниями и непрерывным временем. Для классификации СП по виду статистической связи между значения-
ми СП нам потребуется ввести полное статистическое описание случайного процесса. Говорят, что имеется полное статистическое описание СП (t), если для любых п и t1, t2, …, tn T можно задать функцию распределения
|
x , x |
|
, , x |
|
|
|
, , t |
|
P |
n |
(t |
|
F |
2 |
n |
; t , t |
2 |
n |
|
) x . |
|||||
|
1 |
|
1 |
|
|
i |
i |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i 1 |
|
|
При этом должны быть выполнены условие симметрии, состоящее в том, что
F x1, x2, , xn; t1, t2, , tn F xi1 , xi2 , , xin ; ti1 , ti2 , , tin ,
где i1, i2, …, in – перестановка чисел 1, 2, …, п; и условие согласованности:
F x1, , xi 1, , xi 1, , xn ; t1, t2 , , ti 1, ti , ti 1, , tn
F x1, , xi 1, xi 1, , xn ; t1, t2 , , ti 1, ti 1, , tn для i .
-62 -
Столь же полное вероятностное описание дадут многомерная ПВ:
W x1, x2 , , xn ; t1, t2 , , tn n F xi1 , xi2 , , xin ; ti1 , ti2 , , tin
x1 xn
и многомерная ХФ:
|
|
|
|
v1, v2 , , vn ; t1, t2 , , tn W x1, x2 , , xn ; t1, t2 , , tn |
|||
|
|
|
|
|
n |
|
|
|
|
|
dxn , |
exp j vk xn |
dx1dx2 |
||
|
k 1 |
|
|
или кумулянтная функция:
v1, v2 , , vn ; t1, t2 , , tn ln v1, v2 , , vn ; t1, t2 , , tn .
(t) |
(t) |
|
|
|
t |
|
t |
(t) |
1 |
(t) |
2 |
|
|
t |
t |
3 |
4 |
Рис. 4.1
Для дискретных случайных последовательностей (цифровых сигналов), когда множество возможных значений дискретного сигнала конечно или счетно, вероятностное описание дается распределением вида
P k1, k2 , , kn ; t1, t2 , , tn Bep t1 xk1 , t2 xk2 , , tn xkn ,
где х1, х2, …, хп – возможные (квантованные) значения процесса (сигнала), а k1, k2, …, kп – целые числа.
- 63 -
Как и для описания многомерных случайных величин, для описания случайных процессов могут быть использованы моменты. Начальные и центральные моменты определяются, соответственно, выражениями:
|
|
|
|
|
|
|
|
|
mk , ...,k |
n |
(t1, |
, tn ) |
xk1 xnkn W x1, |
, xn ; t1, , tn dx1 dxn , |
|||
1 |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x1 M 1 k1 |
xn M n kn |
M k , ...,k |
n |
(t1, , tn ) |
|
|||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
W x1, , xn ; t1, , tn dx1 dxn , |
||||
где k1, k2, …, kп – неотрицательные целые числа.
Приведенные выражения зависят от моментов времени t1, t2, …, tn и поэтому называются моментными функциями. Наиболее часто для описания СП используют математическое ожидание (среднее значение) случайного процесса
|
|
|
|
|
(t) m1(t) |
xW (x; t) dx , где t T и корреляционную функцию (КФ): |
|||
|
|
|
|
|
|
|
|
|
x1 m1(t1) x2 m1(t2 ) W (x1, x2; t1, t2 ) dx1dx2 , |
K (t1, t2 ) M11(t1, t2 ) |
|
|
||
где t1, t2 T . Таким образом, корреляционная функция определяет зависи-
мость ковариации случайных величин 1 = (t1) и 2 = (t2) от моментов вре-
мени t1 и t2, в которых берутся отсчеты СП (t). Как и для случайных величин, можно ввести коэффициент корреляции
r t1, t2 |
|
K t1, t2 |
|
, |
|
|
|
|
|
||
|
|
M 2 t1 M 2 t2 |
|
|
|
где М2(t1) и М2(t2) – дисперсии отсчетов процесса в моменты времени t1 и t2.
Напомним, что |
|
r t1, t2 |
|
1. |
|
|
Для того чтобы функция K(t1, t2) была корреляционной функцией случайного процесса, она должна быть ядром симметрического положительно определенного интегрального оператора, т. е. K(t1, t2) = K(t2, t1) и
K (t1, t2 )x(t1)x(t2 )dt1dt2 0
для любых x(t) L2, откуда следует, что K(t, t) > 0.
- 64 -
Если имеется два случайных процесса (t) и (t), то можно ввести в
рассмотрение взаимную корреляционную функцию:
|
|
|
x1 M (t1) y2 M (t2 ) W ( x1, y2; t1,t2 ) dx1dy2 . |
K (t1, t2 ) |
|
|
Эта функция также положительно определена , но не симметрична, так как
K (t1, t2) = K (t2, t1).
Теперь мы подготовлены к тому, чтобы дать классификацию СП по их вероятностным характеристикам.
Наиболее простым с точки зрения вероятностного описания будет СП, у которого любая совокупность значений (отсчетов), взятых в произвольные моменты времени, принадлежащие множеству Т, независима, т. е.
n |
xi ; ti , |
|
F x1, x2 , , xn ; t1, t2 , , tn F i |
|
|
i 1 |
|
|
n |
|
n |
W x1, , xn ; t1, , tn W i xi ; ti , 1, , n ; t1, , tn |
i i ; ti . |
|
i 1 |
|
i 1 |
Существуют СП, у которых статистические связи распространяются лишь на два соседних отсчета. К их числу относится марковский процесс, для которого
W xi ; ti / xi 1, , x1; ti 1, , t1 W xi ; ti / xi 1; ti 1 .
Пользуясь представлением многомерной ПВ через условные ПВ, можно по аналогии записать
|
|
|
|
|
|
|
|
|
|
n |
|
|
ti / xi 1; ti 1 , |
W x1, , xn ; t1, , tn W x1; t1 W xi ; |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
i 2 |
|
||
а так как |
W |
x |
; t |
|
/ x |
i 1 |
; t |
i 1 |
|
W xi 1, xi ; |
ti 1, ti |
, то ясно, что для мар- |
|
i |
|
|
|
||||||||||
|
|
i |
|
|
|
|
W xi 1; |
ti 1 |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|||
ковского процесса многомерная ПВ выражается через двумерные ПВ. Другим примером СП такого вида является процесс с независимыми
приращениями, для которого СВ 2 – 1, 3 – 2, …, n – n 1 независимы для любых п и t1, t2, …, tn T. Напомним, что i – значение СП в момент ti. Для того, чтобы записать ПВ любого порядка для процесса с независимыми
Положительная определенность не означает, что функции K (t1, t2) и K (t1, t2) не могут принимать отрицательные значения.
- 65 -
приращениями, достаточно знать лишь одномерные ПВ (t) и (ti-1) – (ti), а так как ПВ разности двух СВ выражается через их совместную ПВ, то для полного вероятностного описания СП с независимыми приращениями достаточно знания двумерной ПВ.
Стационарные СП. Случайный процесс (t) называется строго ста-
ционарным (стационарным в узком смысле), если его ФР, ПВ или ХФ при любых , п и t1, t2, …, tn T инвариантны к сдвигу моментов t1, t2, …, tn на величину , т. е. W x1, , xn ; t1, , tn W x1, , xn ; t1 , , tn .
Аналогичные выражения можно записать для ФР и ХФ.
СП называется стационарным в широком смысле (по Хинчину), если выполняются более скромные требования:
1) среднее значение процесса m1(t) не зависит от времени, т. е. m1(t) = const; 2) КФ K (t1, t2) зависит лишь от разности t1– t2 = , т. е. K (t1, t2) =
= K (t1– t2) = K ( ).
Из второго условия следует, что дисперсия стационарного процесса
D (t) = K (t, t) = K (0) = const.
В силу симметрии K (t1, t2) корреляционная функция стационарного
СП является четной функцией: K (– ) = K ( ).
При экспериментальных исследованиях СП, когда мы имеем дело с одной реализацией, важным является понятие эргодичности случайного процесса. Дело в том, что для стационарных процессов можно ввести операцию временного усреднения на основе одной реализации x(t) случайного процесса (t). Для обозначения операции временного усреднения часто используют волнистую черту над усредняемым выражением:
|
|
|
~ |
|
|
|
|
|
|
|
|
|
1 T |
|
|
||||
|
|
|
|
|
|
f [x(t)] |
lim |
|
f [x(t)] dt . |
|
|||||||||
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
T T 0 |
|
|
|||||||
Например, для среднего значения |
~ |
|
|
|
|
|
|
|
|||||||||||
x , корреляционной функции B и дис- |
|||||||||||||||||||
персии D t получим: |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
~ |
|
|
|
|
|
|
1 T |
~ |
|
|||||||
|
|
|
|
|
|
x (t) |
lim |
|
|
|
x(t) dt x , |
|
|||||||
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
T T |
0 |
|
|
|
|||||
~~~~~~~~~~~~~~~~~~ |
|
|
1 |
|
T |
|
|
|
|
||||||||||
|
|
~ |
|
|
|
~ |
|
|
|
|
|
|
|
~ |
~ |
|
|||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
x(t) x |
|
|
|
x(t ) x |
|
|
|
lim |
|
|
|
|
x(t) x |
x(t ) x |
dt B( ) , |
||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
T T |
0 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
- 66 - |
|
|
|
|||||
~~~~~~~~~ |
|
|
1 T |
|
|
|
~ 2 |
|
|
~ 2 |
|
||
x(t) x |
|
lim |
|
|
x(t) x |
dt D t B(0) . |
|
||||||
|
|
T T |
0 |
|
|
|
Естественно, возникает вопрос, как соотносятся между собой характеристики процесса, полученные путем статистического усреднения (с помощью ПВ), и средние по времени. В связи с этим дадим определение [6]. Стационарный в узком смысле процесс называется эргодическим, если любые его вероятностные характеристики, найденные на основе статистического усреднения, по множеству реализаций с вероятностью, сколь угодно близкой к единице, совпадают с соответствующими средними по времени. Иными словами, если известна одна-единственная реализация процесса для t , , то путем сдвигов по времени может быть получен бесконечный статистический ансамбль реализаций. Следовательно, по одной реализации можно узнать всевозможные вероятностные характеристики [7].
Часто понятие эргодичности связывают с конкретной вероятностной характеристикой, например со средним значением, называя эргодическим та-
~
кой СП, для которого x m1. При таком определении СП может быть эргодическим, но нестационарным [8].
Пусть, например, (t) m1 cos t , где – СВ, у которой М = 0, а m1 M t . В смысле данного определения процесс является эргодическим, поскольку
~ |
|
|
1 T |
|
(t) m1 |
|
lim |
|
cos t dt m1 . |
|
||||
|
|
T T |
0 |
|
Но он не является стационарным в силу зависимости от времени его дисперсии
M (t) m1 2 M 2 cos2 t .
Мы в дальнейшем будем опираться на первое определение. Эргодичность процесса можно определить на основе условной ПВ
W x2; t2 x1; t1 . Для стационарного в узком смысле процесса одномерная ПВ W (x0; t) не зависит от времени, а двумерная W (x1, x2; t1, t2) зависит только от разности t1 – t2 = , поэтому условная ПВ
W x2 ; t2 |
|
x1; t1 |
|
W x1, x2 ; t1, t2 |
|
|
|||||
|
W x1; t1 |
|
|||
|
|
|
|
||
для стационарного процесса будет равна
- 67 -
|
W x2 |
|
x1; |
W x1, x2 ; |
|
||
|
|
|
|||||
|
|
|
W x1 |
|
|
||
и будет зависеть только от t1 – t2 = . |
|
|
|||||
|
|
|
|
||||
Если |
lim W x2 / x1; W x2 |
и не зависит от х1, |
то процесс (t) |
||||
|
|
|
|
|
|
|
|
называется эргодическим. Иными словами, требуется, чтобы отсчеты стационарного процесса в любые моменты времени t1 и t2 при t1– t2 становились бы независимыми. Для нормального случайного процесса, ПВ отсчетов которого задается многомерным нормальным распределением (см. гл. 3), это условие выполняется, если корреляционная функция стремится к нулю, когда
, т. е. lim K ( ) = 0.
Для эргодического процесса имеет место важное свойство, используемое для экспериментального определения ПВ и ФР. Оно связано с временем пребывания эргодического процесса между двумя уровнями х1 и х2
Tx1x2 ti (рис. 4.2).
x(t)
x2 x1
t
t1 |
t2 |
t3 |
ti |
|
|
|
T |
Рис. 4.2
Предел отношения Tx1x2
T при Т равен вероятности попадания
отсчетов СП в промежуток (х |
, |
х ), т. е. |
P x |
x x |
2 |
lim |
Tx1x2 |
. Если |
|||
|
|||||||||||
1 |
|
|
2 |
|
1 |
|
T |
T |
|||
|
|
|
|
|
|
|
|
|
|||
х1 – х2 = х достаточно мало, то Р(х1 < x < х2) = W (x) x и |
|
|
|||||||||
|
|
|
1 |
|
Tx x |
2 |
|
|
|
|
|
W (x) |
|
lim |
1 |
. |
|
|
|
|
|||
|
|
|
|
|
|
|
|||||
|
|
|
x T |
T |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
- 68 -
При х1 = – получаем возможность для оценки ФР, так как Р( < x < х2) = = F(x2). Важную роль в радиотехнике играют квазидетерминированные про-
цессы. Квазидетерминированный процесс определяется как совокупность
функций времени t заданного вида s(t; ) , зависящих от случайного парамет-
ра 1, 2, , n – случайного вектора. Иногда такие СП называют элементарными случайными функциями.
Например, s(t; A, , ) Asin t , где А, , – случайные величины с совместной ПВ W ( А, , ).
Другим примером квазидетерминированного процесса, часто встречающимся в задачах радиотехники, является сигнал вида
N 1
s(t; A0 , A1, , AN 1, 0 , 1, , N 1) Ai S0 (t iT ) cos 0t i , i 0
представляющий собой группу (пачку) из N радиоимпульсов с несущей частотой 0, сдвинутых друг относительно друга на величину Т, имеющих де-
терминированную огибающую S0(t) и случайные амплитуды Ai, Ai > 0 и фа- |
|||||||||
зы i, i [– , ], определяемые ПВ |
|
||||||||
W |
(t, A , A , , A |
1 |
, |
0 |
, , , |
N 1 |
) . |
8 |
|
|
0 1 |
N |
|
1 |
|
|
|||
|
При прохождении сигнала s(t) мно- |
7 |
|||||||
голучевого |
канала |
на входе приемника |
|
||||||
|
|
|
|
|
|
M |
|
|
|
наблюдается процесс вида Ai s(t i ) , |
4 |
||||||||
|
|
|
|
|
|
i 1 |
|
|
5 |
где {Ai; i}, i = 1, 2, …, M – совокупность |
3 |
||||||||
|
|||||||||
случайных амплитуд Ai, определяемых |
2 |
||||||||
6 |
|||||||||
ослаблением сигнала в соответствующем |
1 |
||||||||
луче, и случайных задержек i, |
связан- |
Рис. 4.3 |
|||||||
|
|
|
|
|
|
|
|
|
|
ных с запаздыванием сигнала при распространении. Рассмотренная классификация может быть образно представ-
лена с помощью «галактики» случайных процессов [8], приведенной на рис. 4.3. На рис. 4.3 приняты следующие обозначения: 1 – нестационарные СП; 2 – стационарные в широком смысле СП; 3 – стационарные в узком смысле
СП; 4 – эргодические СП; 5 – квазидетерминированные СП; 6 – СП с незави-
- 69 -
симыми значениями; 7 – марковские СП; 8 – СП с независимыми приращениями.
Контрольные вопросы
1.Дайте определение случайной функции, случайного процесса.
2.По каким признакам осуществляется классификация СП?
3.Что дает полное статистическое описание СП?
4.Запишите выражения для математического ожидания m1 t и корре-
ляционной функции K t1,t2 СП.
5.Как определяется взаимная корреляционная функция СП t и n t ?
6.Дайте определение СП с независимыми отсчетами, марковского СП, СП с независимыми приращениями.
7.Дайте определение СП, стационарного в широком и узком смыслах.
8.Как определяется операция усреднения по времени для стационарного СП?
9.Запишите корреляционную функцию стационарного СП на основе усреднения по времени.
10.Дайте определение эргодического СП.
11.Как для эргодического СП можно получить оценки ПВ и ФР?
12.Что называется квазидетерминированным СП? Приведите примеры.
Глава 5. КОРРЕЛЯЦИОННЫЕ И СПЕКТРАЛЬНЫЕ ХАРАКТЕРИСТИКИ СЛУЧАЙНЫХ ПРОЦЕССОВ
Для детерминированных сигналов исчерпывающим описанием является задание функции s(t), определяющей поведение сигнала во временной области или преобразования Фурье функции s(t), при условии, что оно суще-
|
|
|
ствует, |
~ |
s(t) exp j2 ft dt , называемое спектром сигнала. Для слу- |
s ( f ) |
||
|
|
|
чайных процессов аналитическое задание каждой реализации как функции времени, за исключением квазидетерминированных СП, невозможно, но некоторое представление о поведении реализации СП во времени можно получить с помощью КФ рассматриваемого процесса. Напомним, что КФ показывает, как меняется корреляционная связь между отсчетами СП, расположенными в моменты времени t1 и t2. Если эта связь быстро изменяется, то соот-
- 70 -
ветственно быстро будет меняться и КФ. Реализации такого процесса будут быстро меняющимися функциями времени. Для иллюстрации на рис. 5.1 приведены реализации (слева) и КФ (справа) для двух стационарных СП.
|
K 1( ) |
|
t |
|
|
K 2( ) |
|
|
|
|
|
t |
|
|
Рис. 5.1 |
|
|
Дисперсии процессов 1(t) и 2(t) предполагаются одинаковыми. Пре- |
||
|
|
|
~ |
x(t) exp j2 ft dt , если оно су- |
|
образование Фурье реализаций СП x ( f ) |
||
|
|
|
ществует, является случайной функцией переменной f |
и, как и сама реали- |
|
зация x(t), не может быть использована для описания СП. |
|
|
Для стационарного СП в качестве характеристики в частотной области |
||
используется преобразование Фурье от КФ |
|
|
|
|
|
S( f ) K ( ) exp j2 f d , |
(5.1) |
|
называемое спектральной плотностью СП. Если СП (t) описывает случайный ток, случайное напряжение или случайную напряженность поля, то S(f) называют спектральной плотностью мощности (СПМ) или энергетическим спектром флуктуаций процесса (t). Флуктуациями называют отклонение значений процесса (t) от среднего значения.
Корреляционная функция выражается через спектральную плотность мощности с помощью обратного преобразования Фурье:
|
|
|
K ( ) |
S( f ) exp j2 f df , |
(5.2) |
Формулы (5.1) и (5.2) определяют содержание теоремы Винера–Хинчина, утверждающей, что КФ и спектральная плотность связаны друг с другом преобразованием Фурье.
- 71 -
