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

Nugmanov

.pdf
Скачиваний:
10
Добавлен:
10.02.2015
Размер:
1.6 Mб
Скачать

 

 

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ω ,

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

-∞

 

 

 

сделаем

зеркальное

отображение

 

 

 

спектральной плотности мощности по

 

 

0

 

 

 

 

отношению к точке ω0 и переместим

 

Рис. 2.1

 

 

 

начало координат в точку ω0

при

 

 

 

 

помощи преобразования

ν ω0 ω , (Рис. 2.1). Это преобразование

позволяет рассматривать спектральную плотность мощности в области низких частот.

В результате получим

 

 

 

 

 

 

 

 

 

 

1

 

ω0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B(τ)

 

F0 (ν) сos (ω0 ν)τ dν

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω0

 

 

 

 

 

 

 

 

 

 

ω0

 

 

 

 

 

 

 

1

 

F

(ν) сosντ dν

 

сosω

 

τ

1

 

F

(ν) sin ντ dν

sin ω

 

τ .

 

 

 

0

 

0

 

0

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Введем обозначения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

ω0

 

 

 

 

 

 

 

 

 

1

 

ω0

 

 

 

a c (τ)

 

F0 (ν) cos ντ dν;

 

a s (τ)

 

 

F0 (ν) sin ντ dν .

 

(2.30)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Функции 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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]