Nugmanov
.pdf
|
|
w ξ (x1, , x n , t1, , t n ) |
|
|
|||||
|
n / 2 n |
|
|
n |
n |
|
|||
exp |
2 Di j (xi m ) (x j m ) . (2.17) |
||||||||
|
|
1 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2 ) |
|
D |
|
2 D i 1 |
j 1 |
|
Как видно из (2.6) и (2.17) совместная плотность распределения вероятности w ξ (x1, , x n , t1, , t n ) зависит от нормированной
ковариационной матрицы R. Для стационарного процесса в широком смысле она будет равна
1 |
||
|
|
t1) |
R(t 2 |
||
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
t1) |
R(t n |
R(t1 t 2 )
1
R(t n t 2 )
R(t1 |
|
t n ) |
|
|
|
|
|
|
|
R(t 2 |
t n ) |
. (2.18) |
||
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
1 |
|
|
|
|
|
Если сдвинуть все точки отсчёта по времени на одну и ту же величину , (t1 , , tn ) , то элементы нормированной
ковариационной матрицы не изменятся, и многомерная плотность распределения не будет зависеть от сдвига . Отсюда следует, если нормальный случайный процесс стационарен в широком смысле, то он будет стационарным и в узком смысле.
Как видно из выше изложенного, нормальный случайный процесс содержит в качестве параметров математическое ожидание и ковариационную матрицу.
2.4 Каноническое разложение случайного процесса
Анализ случайных процессов при их прохождении через радиотехнические цепи бывает сложным из-за математических трудностей. На практике пользуются методами, упрощающими вычислительный процесс. Одним из этих методов является каноническое разложение случайного процесса.
Вводится элементарный случайный процесс ξ0 (t) V φ(t) , где
V - случайная величина, φ(t) - не случайная функция времени, называемая координатной функцией. Определим математическое ожидание и корреляционную функцию процесса ξ0 (t) :
53
Mξ (t) m |
|
|
φ(t) , B |
ξ |
|
(t |
1 |
, t |
2 |
) M(ξ |
0 |
(t |
1 |
) ξ |
0 |
(t |
2 |
)) M(V2 ) φ(t |
1 |
) φ(t |
2 |
), |
||||||||||||||
0 |
|
|
|
|
v |
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
ковариационная функция равна |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
o |
|
|
|
|
|
|
o |
|
|
|
)) (M(V2 ) m |
2 ) φ(t |
|
|
|
|
|
|
|
|
|||||||||
B o (t |
1 |
, t |
2 |
) M(ξ |
0 |
(t |
1 |
) ξ |
0 |
(t |
2 |
1 |
) φ(t |
2 |
) . |
|
|
|
|
|||||||||||||||||
|
ξ0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
v |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Как видно из этих выражений, в общем виде элементарный случайный процесс не является стационарным хотя бы в широком смысле.
Представим случайный процесс ξ(t) через элементарные случайные процессы
N |
|
ξ(t) mξ (t) Vi φi (t) , |
(2.19) |
i 1 |
|
где mξ (t) - не случайная функция, Vi - коэффициенты разложения – |
некоррелированные случайные величины с математическим
ожиданием, равным нулю и дисперсией σ 2 |
: |
|
|
|
vi |
|
|
0 |
при j i , |
|
|
|
2 |
|
(2.20) |
M(Vi Vj ) |
|
||
D(Vi ) σ |
vi |
при j i . |
|
|
|
|
Множество координатных функций может быть как конечным, так и бесконечным. Представление случайного процесса в виде (2.19) с ограничениями (2.20) называется каноническим разложением
случайного |
процесса. |
Для |
центрированного |
процесса |
о |
|
|
|
|
ξ (t) ξ(t) mξ (t) имеем |
|
|
|
|
|
|
о |
N |
|
|
|
Vi φi (t) . |
|
|
|
|
ξ (t) |
(2.21) |
|
|
|
|
i 1 |
|
Возможно интегральное представление |
|
|||
|
(t) m (t) V( ) (t, ) d, |
(2.22) |
||
|
|
|
|
|
|
|
|
|
где V(τ) - белый шум с математическим ожиданием M (V(τ)) 0, φ(t, τ) - неслучайная функция времени t и параметра τ .
Такое представление случайного процесса называется
интегральным каноническим разложением.
54
Будем полагать в дальнейшем в этом разделе, что процесс ξ(t) -
|
o |
|
|
|
центрированный, т.е. ξ(t) ξ (t) , и |
выразим |
корреляционную |
||
функцию процесса ξ(t) , используя каноническое разложение: |
||||
|
|
N |
N |
|
|
|
|
|
|
Bξ (t1, t2 ) M ξ(t1)ξ( t2 ) M |
Vi φi (t1) Vj φ j(t2 ) |
|||
|
i 1 |
i 1 |
|
|
N |
N |
N |
|
|
M Vi Vj φi (t1) φ j(t2 ) M Vi2 φi (t1) φi (t2 ) |
||||
i 1 |
j 1 |
i 1 |
|
|
|
N |
|
|
|
|
D Vi φi (t1) φi (t2 ) . |
(2.23) |
i 1
Разложение (2.23) называется каноническим разложением ковариационной функции. Дисперсия процесса ξ(t) , выраженная через координатные функции, имеет вид
N |
φ i2(t1 ) . |
|
D(ξ) Bξ (t1 t 2 ) D Vi |
(2.24) |
i 1
Доказано [7], если ковариационная функция случайного процесса ξ(t) имеет вид (2.23), то случайный процесс ξ(t) может быть представим в виде (2.21). На практике часто применяют это утверждение. Выбор координатных функций φi (t) зависит от свойств случайного процесса ξ(t) и набора известных случайных величин V1, V2 , , VN . Более подробно теория и практика канонического разложения случайной функции изложена в [7].
Рассмотрим стационарный случайный процесс ξ(t) в интервале
(0, Tн ). Корреляционная функция |
Bξ (τ) может быть представима |
|||
разложением в ряд Фурье на интервале ( Tн τ Tн ) : |
|
|||
|
|
|
||
Bξ (τ) Bξ (t 2 t1 ) D0 Di сos ωi τ , |
|
|||
|
i 1 |
|
||
|
|
|
|
|
D0 Di сos (ωi (t 2 t1 )) , |
ωi i ω 0 i |
π |
. |
(2.24) |
|
||||
i 1 |
|
Tн |
|
|
|
|
|
|
Коэффициенты разложения определяются как
55
Tн
D0 1 Bξ (τ) dτ , (2.25) 2 Tн Tн
Tн
Di T1 Bξ (τ) сos ω i τ dτ , i 1, 2, , (2.26)
н Tн
Перепишем выражение (2.24) в виде
|
|
|
|
Bξ (t 2 t1 ) D0 Di сos (ωi (t 2 t1 )) |
|
||
|
i 1 |
|
|
|
|
|
|
D0 Di |
сos ωi t1 сos ωi t 2 Di |
sin ωi t1 sin ωi t 2 . |
(2.27) |
i 1 |
i 1 |
|
|
В результате |
получили разложение |
ковариационной |
функции |
случайного процесса ξ(t) по координатным функциям cos ωi t и sin ωi t на интервале ( Tн , Tн ) . Ввиду взаимной однозначности представления канонического разложения случайного процесса ξ(t) и его ковариационной функции имеем
|
|
|
ξ(t) U0 Ui cos ωi t Vi sin ωi t , |
(2.28) |
|
i 1 |
i 1 |
|
где U0 , Ui , Vi - не коррелированные случайные величины с математическими ожиданиями и дисперсиями соответственно
M(U0 ) M( Ui ) M(Vi ) 0 ,
D(U0 ) D0 , D( Ui ) D(Vi ) Di .
Представление (2.28) называется спектральным разложением случайного процесса ξ(t) на интервале (0, Tн ) . Наличие случайной величины U0 свидетельствует о том, что процесс ξ(t) не является эргодическим, так как нарушается условие эргодичности Слуцкого. Если коэффициент D0 в (2.28) будет равен нулю, то процесс ξ(t) будет эргодическим.
56
2.5 Квазидетерменированный случайный процесс
Случайный процесс ξ(t, A, , ) называется квазидетерменированным, если он описывается неслучайной функцией времени, содержащий в качестве параметра одну или несколько случайных величин с известными совместными плотностями распределения вероятности.
Случайность проявляется в том, что в момент начала отсчета случайные параметры принимают одно из своих возможных значений и в дальнейшем не изменяются во время наблюдения процесса, т.е. процесс в дальнейшем будет детерминированным. Многомерная плотность распределения вероятности зависит от распределения вероятности в момент t t1. Выразим многомерную плотность распределения вероятности квазидетерменированного процесса через одномерную плотность и условную плотности распределения вероятности:
wξ (y1, , yn , t1, t n )
wξ (y1, t1) wξ (y2 , , yn , t 2 , t n / y1, t1) (2.29)
Положим, что найдена одномерная плотность распределения w ξ (y1, t1) . Ввиду того, что значения случайных величин в моменты
времени t 2 , t n известны и неизменны, то можно использовать -
функцию для |
описания условной плотности |
распределения |
wξ (y2 , , yn , t 2 , t n / y1, t1) : |
|
|
|
n |
|
w ξ (y2 , , yn , t 2 , t n / y1, t1) δ(yi Q(ti / y1, t1)) , |
||
|
i 2 |
|
где Q(ti / y1, t1) |
- функция, описывающая процесс |
ξ(t) в моменты |
времени t t 2 |
при условии, |
что значение процесса в момент t t1 |
|
известно. |
|
|
|
Пример 2.1. |
Случайный |
процесс ξ(t) A сos ( t ) содержит |
|
случайные величины амплитуду A , частоту и |
фазу , не |
||
зависящие от времени. Для простоты будем считать |
амплитуду A |
случайной величиной с известной плотностью распределения wA (x), x1 x x2 , а частоту и фазу будем считать известными, т.е.
57
ω, 0, |
ξ(t) A cos ωt . |
|
Случайный |
|
процесс |
принимает |
||||||||||||||||||
значения |
y(t) x cos ωt . |
В |
|
момент |
времени |
|
|
|
|
t t1имеем |
||||||||||||||
y(t1) x cos ωt1, |
и амплитуда |
|
сигнала будет |
|
равна |
x |
y(t1) |
. |
||||||||||||||||
|
|
cos ωt1 |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Одномерная плотность распределения wξ (y1, t1) будет равна |
|
|
||||||||||||||||||||||
w ξ (y1, t1) w A ( |
y1 |
|
) |
|
dx |
|
w A ( |
y1 |
|
) |
|
1 |
|
|
|
. |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
cos ωt |
1 |
|
dy |
|
cos ωt |
1 |
|
|
cos ωt |
1 |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Функция Q(ti / y1, t1) равна: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
Q(ti / y1, t1) |
|
|
y(t1) |
cos ωti . |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
cos ωt1 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таким образом, многомерная плотность распределения примет вид w(y1, , yn , t1, , t n )
|
y(t1) |
|
|
|
|
n |
y(t1) |
|||
w A ( |
1 |
|
|
δ(yi |
||||||
|
) |
|
|
|
|
|
|
cos ωt i ). |
||
cos ωt1 |
|
|
cos ωt1 |
|
|
cos ωt1 |
||||
|
|
|||||||||
|
|
|
|
|
|
|
|
i 2 |
|
|
2.6 Узкополосный случайный процесс
По свойствам спектральной плотности мощности случайные процессы могут быть разделены на широкополосные и узкополосные.
Случайный процесс ξ(t) называется узкополосным, если его спектральная плотность мощности сосредоточена вблизи какой-
либо частоты, (Рис 2.1). Другими словами, узкополосный процесс характеризуется полосой частот ω0 , в пределах которой сосредоточена основная мощность процесса по отношению к частоте ω0 , т.е. признаком узкополосности процесса является соотношение
ω ω0 .
58
Определим ковариационную функцию процесса, исходя из
|
F( ) |
|
определения. |
Считаем, |
что |
|||
|
|
спектральная |
плотность |
мощности |
||||
|
|
|
|
|||||
|
|
|
|
F(ω) известна |
|
|
||
|
|
|
|
B(τ) 21π |
|
|
|
|
|
|
|
|
F(ω) exp( jω τ) dω. |
|
|||
|
|
|
|
|
|
|
|
|
- 0 |
0 |
0 |
|
Перейдем |
к полосе |
частот |
в |
|
|
|
|
пределах (0, ) |
|
|
|||
|
|
F( ) |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
B(τ) |
F0 (ω) сos ωτ dω , |
|
||
|
|
|
|
|
2π |
0 |
|
|
|
|
|
|
|
|
|
|
|
-∞ |
|
|
|
сделаем |
зеркальное |
отображение |
||
|
|
|
спектральной плотности мощности по |
|||||
|
|
0 |
|
|||||
|
|
|
отношению к точке ω0 и переместим |
|||||
|
Рис. 2.1 |
|
||||||
|
|
начало координат в точку ω0 |
при |
|||||
|
|
|
|
|||||
помощи преобразования |
ν ω0 ω , (Рис. 2.1). Это преобразование |
позволяет рассматривать спектральную плотность мощности в области низких частот.
В результате получим
|
|
|
|
|
|
|
|
|
|
1 |
|
ω0 |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
B(τ) |
|
F0 (ν) сos (ω0 ν)τ dν |
|
|
|||||||||||
|
|
|
|
|
|
|
2π |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
ω0 |
|
|
|
|
|
|
|
|
|
|
ω0 |
|
|
|
|
|
|
|||
|
1 |
|
F |
(ν) сosντ dν |
|
сosω |
|
τ |
1 |
|
F |
(ν) sin ντ dν |
sin ω |
|
τ . |
|||||||
|
|
|
0 |
|
0 |
|||||||||||||||||
|
0 |
|
|
|
|
|
|
0 |
|
|
|
|||||||||||
|
2π |
|
|
|
|
|
|
|
|
|
2π |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Введем обозначения |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
1 |
|
ω0 |
|
|
|
|
|
|
|
|
|
1 |
|
ω0 |
|
|
|
a c (τ) |
|
F0 (ν) cos ντ dν; |
|
a s (τ) |
|
|
F0 (ν) sin ντ dν . |
|
(2.30) |
|||||||||||||
|
2π |
|
2π |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Функции a c (τ) и as (τ) - медленно меняющиеся функции времени, так как их спектры принадлежат низкочастотной области.
59
В новых обозначениях имеем
B(τ) a c (τ) cos ω0τ as (τ) sin ω0τ a(τ) cos (ω0τ φ(τ)) , (2.31)
где ac (τ) a(τ) cos φ(τ) , |
|
as (τ) a(τ) sin φ(τ),. |
|
|
|
|||||||
|
|
|
|
|
|
|
|
a |
|
(τ) |
||
|
a |
|
(τ) 2 |
a |
|
(τ) 2 , |
s |
|||||
a(τ) |
с |
s |
φ(τ) arc tan |
|
|
. |
||||||
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
a |
с (τ) |
Как видно из (2.30) и (2.31), корреляционная функция узкополосного процесса является гармоническим колебанием с частотой ω0 , модулированным низкочастотным сигналом a(τ), фаза
F( ) |
|
гармонического колебания φ(τ) зависит от τ . |
|||||
|
|
|
|
|
|
|
|
|
|
На практике обычно спектральная плотность |
|||||
|
|
мощности F(ν) бывает симметричной функцией |
|||||
|
|
(Рис. |
2.2). В |
этом |
случае, |
если пренебречь |
|
|
|
мощностью |
процесса в области |
(ω0 , ) , |
|||
0 |
0 |
синусная составляющая as (τ) |
обратится в нуль |
||||
Рис. 2.2 |
|
и |
получим |
простое |
|
выражение |
|
|
B(τ) |
a c (τ) сos ω0τ . |
|
|
|
||
|
|
|
|
|
Вчастности, если спектральная плотность мощности
определяется как (Рис 2.3а , f0 50 Гц , 30 1/ Сек )
|
2 α (α2 ω2 ω2 ) |
|
|
|
|
|
||||
|
|
|
|
|
|
|||||
F(ω) |
|
|
|
0 |
|
|
, |
ω |
, |
|
(α2 (ω ω |
0 |
)2 ) (α2 |
(ω ω |
0 |
)2 ) |
|||||
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
F( ) |
B( ) |
|
|
|
- |
0 |
0 |
|
Рис. 2.3а |
Рис. 2.3б |
60
корреляционная функция случайного процесса имеет вид
B(τ) e α τ cos ω0 τ ,
которая изображена на Рис. 2.3б.
Представление корреляционной функции в виде (2.31) позволяет записать модель узкополосного процесса в виде
ξ(t) A(t) сos(ω0 t (t)) , |
(2.32) |
где A(t) - огибающая узкополосного случайного процесса, является также случайным процессом, (t) - случайная фаза узкополосного случайного процесса.
Представим случайный процесс (2.32) в виде
ξ(t) A(t) сos (t) сosω0 t A(t)sin (t)sin ω0t
Ac (t) cos ω0 t As (t) sin ω0 t , |
(2.33) |
|
где |
Ac (t) A(t) cos (t) , |
|
y(t) |
As (t) A(t) sin (t) |
(2.34) |
|
t |
Одна |
из |
возможных |
реализаций |
|
узкополосного |
случайного |
процесса |
|||
|
|||||
|
представлена на Рис.2.4. |
|
Составляющие Ac (t) и As (t)
Рис. 2.4 |
являются медленно меняющимися |
|
|
|
функциями по сравнению с функциями |
сosω0 t и sin ω0 t . |
|
Рассмотрим плотность распределение вероятности огибающей A(t) и фазы (t) узкополосного процесса. Положим косинусная Ac (t) и синусная As (t) составляющие являются стационарными
случайными процессами хотя бы в широком смысле, взаимно независимы и распределены по нормальному закону с соответствующими математическими ожиданиями (mc , ms ) и
дисперсиями σс2 σs2 σ2 :
61
|
|
|
|
|
|
|
1 |
|
|
|
w |
|
(x |
|
, x |
|
) |
|
exp |
|
|
Ac ,As |
c |
s |
|
|
||||||
|
2 |
|||||||||
|
|
|
|
2 π σ |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
где |
xc (t) x(t) cos φ(t) f1, |
|
Ввиду того, что процессы
x c mc 2 |
|
xs ms 2 |
, (2.35) |
|||
|
|
|
|
|
||
2 σ |
2 |
|
2 σ |
2 |
|
|
|
|
|
|
|
xs (t) x(t) sin φ(t) f2 .
Ac (t) и As (t) - стационарные и
рассматриваются |
одномерные |
плотности |
распределения |
вероятности, в дальнейшем обозначим x(t) x, |
φ(t) φ. |
Согласно правилам преобразования координат выразим совместную плотность огибающей и фазы как
w A, (x, φ) w Ac ,As (f1, f2 ) J ,
где якобиан преобразования равен J x . После соответствующих подстановок получим
|
|
|
|
x |
|
|
|
|
|
|
1 |
|
|
x |
|
|
|
|
|
2 |
x |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
w |
A, |
(x, φ) |
|
2 |
exp |
|
|
|
|
2 |
c |
cos φ m |
c |
|
s |
sin φ m |
s |
|
. |
||||||||
|
|
2 π σ |
|
|
|
2 σ |
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
Произведем упрощение показателя экспоненты при помощи |
||||||||||||||||||||||||||
подстановок |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
mc m cosβ, |
ms msin β , |
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
m2 |
m2 |
m2 , |
|
tg β |
ms |
. |
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
c |
|
|
s |
|
|
|
|
|
|
|
|
|
mc |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В результате получим совместную плотность вероятности
w A, (x, φ) 2 πxσ2
exp
1 |
|
|
2 |
|
2 |
|
|
|
|
|
|
|
|||
|
|
x |
|
2 x m cos (φ β) m |
|
|
. (2.36) |
|
2 |
|
|
||||
2 σ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Проинтегрировав по φ , получим одномерную плотность распределения вероятности огибающей узкополосного процесса
|
x |
|
|
x2 |
m2 |
|
1 |
π |
|
x m cos (φ β) |
|
|
|
2 σ2 |
|
|
|
||||||
w A (x) |
e |
|
e |
σ2 |
dφ = |
||||||
|
|
||||||||||
σ2 |
|
|
|
|
2 π |
|
|||||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
π |
|
|
62