МАРТ_Часть2_Главы4_8
.pdf
|
|
|
|
t1 |
k t1 dt1 |
Если |
потребовать |
выполнения |
условия |
K t1, t2 |
|
|
|
|
|
t0 |
|
k k t2 , |
|
t1 |
l t2 dt2 |
|
|
то COV k , l k k t2 |
и с учетом ортогональ- |
||||
|
|
t0 |
|
|
|
ности системы функций |
k t |
получим |
требуемый |
результат |
|
COV k , l 0. |
|
|
|
|
|
Таким образом, координаты функции являются собственными функциями оператора Фредгольма, ядром которого является корреляционная функция процесса. Нетрудно показать (сделать это читателю предоставляется самостоятельно), что с учетом свойств корреляционной функции этот оператор является симметричным и положительно определенным. Следовательно, его собственные функции действительно образуют ортонормальную систему (см. ч. 1, гл. 5). Разложение реализаций СП по собственным функциям оператора, порождаемого корреляционной функцией, называется разложением Карунена Лоева.
Рассмотрим более общий случай, когда огибающая и фаза являются стационарными СП, т. е. (t) = (t) cos( t – (t)). Будем также считать, что
(t) – стационарный СП с произвольной ПВ отсчетов W (x), которая в силу стационарности (t) не зависит от времени. Такие процессы называют квазигармоническими. При сделанных предположениях можно показать [8], что в совпадающие моменты времени огибающая (t) и фаза (t) независимы, а отсчеты фазы распределены равномерно в интервале [– , ] вне зависимости от вида распределения отсчетов процесса (t), т. е. W( , ) = W( )W( ) = = W( )1 2 . Можно также показать [8], что ПВ огибающей W( ) определяется
через характеристическую функцию (v) процесса (t) следующим образом:
W( ) = vJ 0 (v ) (v)dv .
0
Записанное соотношение является обратным преобразованием Ганкеля, связывающим для данного процесса его ХФ (v) и ПВ огибающей W( ), в со-
ответствии с которым (v) W ( ) J 0 (v ) d . При этом ПВ отсчетов
0
- 82 -
процесса (t) должна быть четной функцией W (–x) = W (x), вещественной. Например, если
|
|
1 |
|
|
x |
2 |
|
|
W (x) = |
|
|
exp |
|
|
, |
||
|
|
|
|
|
|
|||
|
2 2 |
|
|
2 |
2 |
|
||
|
|
|
|
|
|
|||
т. е., отсчеты (t) подчинены нормальному закону, то (v) =
для ПВ огибающей W( ) мы получим распределение Рэлея
а ХФ (v) –
|
|
2 |
2 |
|
exp |
|
v |
|
и |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
W( ) = |
|
|
2 |
|
. |
2 2 2
Сдругими примерами квазигармонических процессов можно ознакомиться с, 0exp
помощью [8].
Во многих задачах статистической радиотехники приходится сталки-
ваться с квазипериодическими стационарными СП вида |
|
(t) = (t) F( t + (t)), |
(5.3) |
где (t) и (t) – стационарные СП, (t) > 0, | (t) | < , F( ) – периодическая функция с периодом 2 , т. е. F( t + (t) + 2 k) = F( t + (t)), k = 0, 1, 2, … .
Кроме того, предполагается, что Fmax = 1. При F( ) = cos получаем рассмотренный ранее квазигармонический СП. Процесс вида (5.3) описывает, например, случайную последовательность импульсов, форма которых определяется функцией F, а процессы (t) и (t) описывают флуктуации амплитуды и параметров импульса (длительность, временное положение) соответственно.
Можно показать [8], что, как и для квазигармонического СП, у данного процесса в совпадающие моменты времени отсчеты СП (t) и (t) независимы, а «фаза» распределена равномерно в интервале [– , ].
Рассмотренные примеры СП относились к четвертому разделу введенной в начале гл. 2 классификации СП (случайные процессы с непрерывными состояниями и непрерывным временем). Приведем примеры СП, относящихся к другим разделам.
Телеграфный сигнал. Случайный процесс, называемый телеграфным сигналом, может принимать только два значения: +u и –u, а моменты перемены знака – случайны. Пример реализации такого процесса приведен на рис. 5.4.
- 83 -
(t)
U
t
–U
Рис. 5.4
Будем считать, что число перемен знака в единицу времени подчиня-
ется распределению Пуассона с параметром : |
P M |
M |
. Напо- |
e |
|||
|
|
M! |
|
мним, что – среднее значение дискретной СВ , т. е. среднее число перемен знака процессом (t) в единицу времени. Процесс (t) будет стационарным, если не зависит от времени. ПВ отсчетов данного процесса будет иметь вид:
W (x; t) P (t) 0 (x Um ) P (t) 0 (x Um ) .
Так как состояния Um и –Um равновероятны, что вытекает из стационарности
процесса, то выполняются равенства
P (t) 0 P (t) 0 1 2 ; W (x; t) 1 2 (x Um ) (x Um ) .
Среднее значение и дисперсия процесса (t) будут равны
|
|
|
|
M (t) |
xW ( x) dx 1 2 |
x ( x Um ) ( x Um ) dx 0.5 Um Um 0 ; |
|
|
|
|
|
|
|
|
|
|
D (t) |
x M (t) 2W ( x) dx 1 2 Um2 Um2 Um2 . |
|
|
|
|
|
Этот результат очевиден, так как квадрат любой реализации есть U m2 для всех t, и результат усреднения с учетом равенства нулю среднего значения дает полученное выше значение U m2 .
Найдем теперь корреляционную функцию процесса (t).
K( ) M (t) (t ) M (t) (t ) ,
поскольку 0. Произведение (t) (t ) может принимать два значения.
Это U m2 , если на интервале длиной произошло четное число перемен знака
(М = 0, 2, 4, …) и –U m2 , если на этом интервале имело место нечетное число
- 84 -
перемен знака |
(М = 1, 3, 5, …). При этом – среднее число перемен знака |
||||||||||||||||||||||||
процессом (t) |
на интервале длиной . Вероятности этих событий равны со- |
||||||||||||||||||||||||
|
|
|
|
( ) |
2M |
|
|
|
|
|
|
( ) |
2M 1 |
|
|
|
|
|
|
|
|
||||
ответственно |
|
|
|
e |
и |
|
|
|
|
|
e |
. Таким образом, |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
M 0 (2M )! |
|
|
M 0 |
|
(2M 1)! |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
( ) |
2M |
|
|
|
|
|
|
|
|
|
( ) |
2M 1 |
|
|||||
|
M (t) (t ) U m2 |
|
|
|
|
e –U m2 |
|
|
|
|
e = |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
M 0 |
(2M )! |
|
|
|
|
|
|
M 0 |
|
(2M 1)! |
|
||||||||
|
|
|
|
|
=U m2 e |
|
|
1 M |
( ) |
M |
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
M 0 |
|
|
M! |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
1 M |
( ) |
M |
|
и |
|
|
|
|
|
|
|
|
|||||||
Учитывая равенство |
|
|
|
e |
четность |
корреляционной |
|||||||||||||||||||
M! |
|
|
|||||||||||||||||||||||
|
|
|
|
M 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
функции стационарного процесса, |
получим окончательно K( ) = U 2 |
e 2 | | . |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
m |
|
Как и следовало ожидать, K(0) = D (t) = U |
2 . |
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
m |
|
|
|
|
|
|
|
|
|||
Спектральная плотность, соответствующая полученной K( ), |
равна |
||||||||||||||||||||||||
|
U 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
S( f ) |
m |
|
. Графики K( ) и S(f) приведены на рис. 5.5. |
|
|
|
|||||||||||||||||||
2 2 f 2 |
|
|
|
||||||||||||||||||||||
Как уже отмечалось, КФ не является исчерпывающей характеристикой СП. СП, реализации которых резко различаются по форме, могут иметь одинаковые КФ. Так, например, при прохождении «белого» шума с S( ) N0
2
|
|
|
|
|
1 |
|
|||
через интегрирующую RC-цепь с коэффициентом передачи K ( j ) |
|
, |
|||||||
1 j T |
|||||||||
где T RC – постоянная времени цепи, СПМ на выходе будет иметь вид |
|||||||||
S ( ) S ( ) |
|
K j |
|
2 |
|
N0 2 |
|||
|
|
||||||||
|
|
|
|
. Если положить T 1 2 , то СПМ рассмот- |
|||||
|
|
|
1 T 2 |
||||||
|
|
|
|
|
|||||
ренного ранее телеграфного сигнала и данного СП совпадут по форме. Следовательно, одинаковую форму будут иметь и КФ этих СП. Выбором значения коэффициента N0
2 можно обеспечить полное совпадение КФ и СПМ рассматриваемых процессов. Этот результат понятен. Ведь K и S( ) являются энергетическими характеристиками СП, а форма реализации определяется, в том числе, и фазовым спектром.
- 85 -
K ( ) |
S(f) |
Um2 |
2 |
|
Um |
|
|
0 |
|
0 |
f |
|
Рис. 5.5
Еще раз подчеркнем, что K характеризует скорость изменения СП, которую можно также охарактеризовать средним числом перемен знака СП в единицу времени. Подумать о том, какой вид будут иметь реализации процесса на выходе интегрирующей RC-цепи при подаче на вход «белого» шума, мы предлагаем читателю.
По нашей классификации рассмотренный процесс относится к классу процессов с непрерывным временем и дискретными состояниями. Часто приходится сталкиваться с кусочно-постоянными процессами, у которых изменение состояния происходит в заранее заданные моменты времени, например в моменты t = kT, где k = 0, 1, 2, …, а Т – интервал дискретизации. Такой процесс (t) можно представить как результат дискретизации произвольного процесса (t) с запоминанием отсчетов, взятых в моменты времени t в тече-
ние интервала Т. Реализации процессов (t) |
и (t) приведены на рис. 5.6. |
||||
|
|
(t) |
|
(t) |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
t
Т
Рис. 5.6
Найдем характеристики процесса (t), считая (t) стационарным СП. Очевидно, что математическое ожидание и дисперсия у процессов (t) и (t) одинаковы. Найдем корреляционную функцию процесса (t), считая, что T k и отсчеты процесса (t) можно считать независимыми СВ. Для простоты рассуждений будем также считать, что М{ (t)} = 0. Рассмотрим произведение (t) (t – ). Если > T, то перемножаются независимые СВ и сред-
- 86 -
нее значение произведения равно произведению средних, а K ( ) = 0 при | | > T. Если же < T, то на интервале Т – значения (t) и (t – ) совпадают, и
усреднение дает средний квадрат отсчета, т. е. D (t) = 2. На интервале же длиной перемножаются соседние независимые отсчеты, и после усреднения получается произведение средних значений или 0, в силу условия М{ (t)} = 0.
Таким образом, для корреляционной функции процесса (t) получим:
|
2 |
1 |
|
|
|
T , |
|
|
|
T ; |
|||
|
|
|
|
||||||||||
K ( ) = |
|
|
|
|
|
|
|||||||
|
|
|
|
T . |
|||||||||
0, |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Можно показать (мы предоставляем это читателю), что для общего |
|||||||||||||
случая корреляционная функция процесса K ( ) представляет собой резуль-
тат линейной аппроксимации K ( ) с шагом Т, как это показано на рис. 5.7.
На рис. 5.7. K ( ) изображена основной линией, а K ( ) при разных значениях Т – пунктирной.
Спектральная плотность мощности для T k имеет вид
S( f ) |
2 |
sin fT |
2 |
||
|
T |
|
. |
||
|
fT |
||||
|
|
|
|
|
|
Найти спектральную плотность для общего случая мы предлагаем читателю.
|
K ( ), K ( ) |
|
Заканчивая первое знакомство с |
||
|
|
видами случайных процессов, рас- |
|||
2 |
|
|
|||
|
|
|
смотрим импульсные случайные про- |
||
|
|
|
цессы, которые часто выступают в |
||
|
|
|
роли помехи для радиотехнических и |
||
|
|
|
телекоммуникационных |
систем. |
|
0 |
T1 |
|
Остановимся, следуя [8], на несколь- |
||
ких моделях импульсных случайных |
|||||
|
|
|
|||
|
T2 |
процессов. |
|
|
|
Рис. 5.7 |
|
1. Одиночный импульс со слу- |
|
|
|
|
|
чайным временным положением яв- |
ляется случайным процессом вида (t) = s(t – t0), где s(t) – заданная функция,
а t0 – случайная величина с ПВ W(t0).
- 87 -
|
N 1 |
2. Случайная импульсная последовательность (t) = |
Ak s(t tk ) , где |
|
k 0 |
s(t), как и ранее, заданная функция, а Аk и tk – случайные величины с сов-
местной ПВ W(А0, А1, …, AN–1; t0, t1, …, tN–1). Такой СП является моделью многолучевого канала распространения сигнала в точку приема. Случайные
величины Аk и tk характеризуют ослабление сигнала и его запаздывание при распространении по k-й траектории (лучу). Для локационных задач данный процесс моделирует отражение зондирующего сигнала от естественных или искусственных отражателей. При этом число отражателей N, как и число лучей, являются дискретными случайными величинами.
Рассмотрим статистические характеристики перечисленных СП. Для одиночного импульса со случайным временным положением среднее значение и корреляционная функция равны соответственно
M (t) s(t t0 ) W (t0 ) dt0 ,
K (t, ) M (t) M (t) (t ) M (t ) M (t) (t )
M (t) M (t ) s(t t0 ) s(t t0 ) W (t0 ) dt0 M (t) M (t ) .
Как видно из приведенных выражений, среднее значение является сверткой импульса и ПВ случайной величины t0, определяющей временное положение импульса, а корреляционная функция K (t, ) представляет собой свертку
ПВ W(t0) и функции s(t0) s(t0 – ).
Если функции s(t0) и W(t0) резко различаются по длительности, то одну из них по отношению к другой можно считать дельта-функцией. Используя ее фильтрующее свойство, получим следующие приближенные выражения для среднего значения и КФ:
M (t) W (t t ) s(t) dt ,
K (t, ) W (t t ) s(t) s(t ) dt M (t ) M (t )
- 88 -
при s << W , где s и W – |
длительности функций s(t0) и W(t0) соответ- |
|
ственно, а t и t – абсциссы максимумов функций s(t) и s(t) s(t – ). |
||
При W << s получим: |
|
|
|
|
|
M (t) s(t t0 ) , K (t, ) s(t t0 )s(t t0 ) , |
||
где t – абсцисса максимума ПВ W(t0).
0
Как видно из приведенных выражений, процесс (t) в общем случае является нестационарным. Он становится приближенно стационарным в широ-
ком смысле, если плотность вероятности W(t0) =1/T для |
|
t0 |
|
T |
2 и Т >> s . |
|
|
|
|||||
Тогда |
|
|
|
|
|
|
|
|
|
|
|
|
|
M (t) 1 T s(t) dt const , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
K (t, ) 1 T s(t) s(t ) dt M (t) 2 K ( ) . |
|
|||||
|
2 . |
|
||||
Естественно, этот результат справедлив для t T 2, T |
|
|||||
|
|
|
|
|
N 1 |
|
Для случайной импульсной последовательности (t) = |
Ak s(t tk ) |
|||||
k 0
будем считать, что число импульсов N, моменты появления tk и их ампли-
тудные множители Ak независимы, т. е.
W(N, А0, А1, …, AN–1; t0, t1, …, tN–1) =
= W(N) W(A0) … W(AN–1) W(t0) W(t1) … W(tN–1),
причем все распределения W(Ak) одинаковы и равны W(A), а распределения
моментов появления W(tk) также одинаковы и равны W(t0). Тогда, полагая, что моменты появления импульсов подчинены равномерному распределению
W(t0) = 1
T и Т >> s , после усреднения по N, Ak и Tk с учетом полученных выражений для среднего значения и корреляционной функции, получим:
|
|
M (t) M (N ) M ( A) 1 T s(t) dt , K ( ) =M(N) M(A2) 1 T s(t) s(t ) dt + |
|
|
|
+ M ( A) 2 M (N 2 ) M (N ) M (N ) 2 1
T s(t) s(t ) dt .
- 89 -
Если число импульсов на интервале длиной Т подчиняется распределению Пуассона с параметром = M (N), а распределение W(A) является симметричным относительно нуля и М (А) = 0, то t 0 и
|
M ( A2 ) |
|
|
|
|
|||
K ( ) = |
|
|
s(t) s(t ) dt . |
|
|
|
||
|
T |
|
|
|
||||
|
|
|
|
|
|
|||
|
|
|
|
|
|
|||
Полученной корреляционной функции соответствует СПМ вида |
||||||||
|
|
|
|
|
|
|
|
|
S ( ) K ( ) exp j d k |
|
~ |
|
2 |
, |
|||
|
|
|||||||
|
s ( ) |
|
|
|||||
|
|
|
|
|
|
|
|
|
где множитель k 2 M ( A2 ) T |
характеризует интенсивность данной мо- |
|||||||
дели помехи как по уровню (множитель М(А2)), так и по средней относительной частоте появления импульсов (множитель
T ).
Контрольные вопросы
1.Как связаны между собой КФ и СПМ стационарного СП?
2.Сформулируйте свойства КФ и СПМ стационарного СП.
3.Назовите числовые характеристики КФ и СПМ.
4.Сформулируйте соотношение неопределенности для стационарно-
го СП.
5.Какой процесс называют «белым» шумом, финитным «белым» шумом? Какой вид имеют K и S для названных процессов?
6.Приведите графики (качественно) K и S для СП с гауссовским (колокольным) энергетическим спектром.
7. Случайный процесс x t um cos 2 f0t , где |
– СВ с ПВ |
|||
W |
1 |
2 , 0,2 |
|
|
|
0, 0,2 |
, является стационарным. Останется ли он стацио- |
||
|
|
|
|
|
|
1 |
, 0, |
|
нарным, если ПВ будет равна W |
|
? Ответ обоснуйте. |
|
|
0, 0, |
|
|
8. |
Дайте определение узкополосного процесса. |
||
9. |
Какой вид в общем случае имеет КФ узкополосного СП? |
||
10. При каких условиях КФ |
узкополосного СП имеет вид |
||
K2 cos 2 f0 ?
11.Дайте определение комплексного СП.
- 90 -
12. Какой вид имеет K стационарного СП с дискретным энергетическим спектром?
13. Как выглядит каноническое представление СП? Как могут быть найдены координатные функции канонического разложения СП?
14. Каким образом можно сформировать стационарный СП с
S S0 |
sin T |
2 |
|||
|
|
|
? |
||
T |
|||||
|
|
|
|
||
15. Проанализируйте свойства импульсной случайной последователь-
N 1
ности t Ak s t tk , где N, Ak и tk – независимые СВ.
k 0
Глава 6. НЕПРЕРЫВНОСТЬ, ИНТЕГРИРУЕМОСТЬ И ДИФФЕРЕНЦИРОВАНИЕ СЛУЧАЙНЫХ ПРОЦЕССОВ
Для детерминированных функций важнейшими свойствами являются непрерывность, дифференцируемость, существование интеграла. Рассмотрим, как они определяются для случайных процессов. Перечисленные свойства связаны с операцией предельного перехода и понятием предела последовательности СВ, которое в отличие от обычного понятия сходимости последовательности, рассматриваемого в математическом анализе, зависит от принятого критерия сходимости. Поясним сказанное более подробно.
Говорят, что последовательность СВ 1, 2, …, п сходится по вероят-
ности к СВ , если для любого > 0 lim P n 0 или n
lim P n 1. Если последовательность n сходится по вероят- n
ности к СВ , то последовательность ПВ W n(x) сходится при п в обыч-
ном смысле к W (x). Если = а – детерминированная величина, то последо-
вательность {W n(x)} сходится к (х – а). Сходимость по вероятности обо-
значают как |
по вер. |
|
|
|
lim n (по вер.) или n . |
|
|
||
|
n |
|
|
|
Другой вид сходимости последовательности СВ – это сходимость в |
||||
среднеквадратическом, определяемая как |
lim M n 2 |
0 |
или |
|
|
|
n |
|
|
|
- 91 - |
|
|
|
