МАРТ_Часть2_Главы4_8
.pdf
l. i. m. n , где обозначение l. i. m. – аббревиатура английского названия этого предела (limit in the mean square). Если = а, то учитывая, что
lim |
M n a 2 = lim |
M n M n M n a 2 |
||||
n |
n |
|
|
|
|
|
|
= lim D |
n |
(M |
n |
a)2 |
|
|
n |
|
|
|
||
|
|
|
|
|
|
|
и тот факт, что предел суммы двух неотрицательных слагаемых равен нулю, только если предел каждого слагаемого равен нулю, получим
lim M n a , |
lim D n 0 . |
n |
n |
Таким образом, предел средних значений СВ п есть а, а предел их дисперсий равен нулю.
Для последовательности ПВ W n (x) полученное условие означает схо-
димость W n (x) к (х – а), т. е. сходимость по вероятности. Обратное утверждение в общем случае неверно, так как среднеквадратическая сходимость предполагает конечность M 2n , а это требование не всегда выполняется.
Теперь мы подготовлены к тому, чтобы дать определение непрерывности и дифференцируемости случайного процесса.
Случайный процесс (t) называется непрерывным в точке t в средне-
квадратическом, если lim M (t t) (t) 2 0. Когда это условие вы-
t 0
полняется для всех t [a, b], процесс (t) называют непрерывным на этом интервале. Необходимым и достаточным условием непрерывности процесса
в точке t является непрерывность его корреляционной функции K (t1, t2) при t1 = t2 = t. Для стационарного процесса это означает непрерывность корреля-
ционной функции K ( ) в нуле.
Случайный процесс дифференцируем в точке t в среднеквадратическом и имеет производную (t), если
lim |
|
|
(t t) (t) |
|
|
2 |
0 . |
|
|
||||||
M |
|
t |
(t) |
|
|
||
t 0 |
|
|
|
|
|
|
Учитывая связь между сходимостью в среднеквадратическом и сходимостью по вероятности, можно утверждать, что СП, непрерывный и диффе-
- 92 -
ренцируемый в среднеквадратическом, будет непрерывным и дифференцируемым и по вероятности.
Необходимым и достаточным условием дифференцируемости случайного процесса будет существование и непрерывность второй смешанной
производной корреляционной функции СП 2K (t1, t2 ) , которая определяет
t1 t2
корреляционную функцию процесса (t), т. е.
K (t1, t2 ) 2K (t1, t2 ) .
Для стационарного процесса
K ( ) 2K ( ) .2
За счет чего в данной формуле появился минус, предлагаем догадаться чита-
~ ~
телю. Учитывая, что F s (t) j s ( ) , где s ( ) F s(t) , а F – оператор Фурье, получим, что спектральная плотность производной
S ( ) 2S ( ) .
Из рассмотренных выше примеров случайных процессов «белый» шум является разрывным, недифференцируемым процессом (см. рис. 5.3, а); процесс с лоренцевским энергетическим спектром (рис. 5.3, е) является непрерывным, но недифференцируемым. Остальные рассмотренные процессы непрерывны и дифференцируемы.
Рассмотрим вопрос о взаимной корреляционной функции процесса (t) и его производной (t). Считая процесс дифференцируемым в среднеквадратическом, можно показать, что
K (t1, t2 ) K (t1, t2 ) .
t2
Для стационарного процесса
K ( ) K ( ) K ( ) .
Преобразование Фурье взаимной КФ даст взаимный энергетический спектр процесса и его производной, который в силу свойств преобразования Фурье будет равен
S ( ) j S ( ) .
- 93 -
Так как корреляционная функция стационарного СП является четной, то K ( ) – нечетная функция, т. е. K ( ) = – K ( ) и K (0) = 0. Таким образом, значения дифференцируемого СП и его производной в совпадаю-
щие моменты времени (t1 = t2 = t, = 0) некоррелированы. Рассмотрим вопрос об интегрировании случайных процессов.
T
Изучим сначала линейный функционал от СП (t) вида s(t) (t) dt ,
0
где s(t) – детерминированная функция. Записанный интеграл есть случайная величина, к которой при max tk 0 в среднеквадратическом сходится последовательность СВ
n |
n |
n s(tk ) (tk ) tk , |
tk T . |
k 1 |
k 1 |
Среднее значение и дисперсия этой случайной величины равны соответственно
T |
T T |
|
M s(t)M (t) dt , D |
s(t1)s(t2 )K (t1, t2 ) dt1dt2 . |
|
0 |
0 |
0 |
Для стационарного случайного процесса |
|
|
М{ (t)} = const; K (t1, t2 ) K (t1 t2 ) , |
||
и, следовательно, |
|
|
T |
T T |
|
M M (t) s(t)dt ; D |
s(t1)s(t2 )K (t1 t2 ) dt1dt2 . |
|
0 |
0 |
0 |
Если М { (t)} = 0, то M 0.
При k << T в силу центральной предельной теоремы распределение
W (x) будет близко к нормальному с параметрами M и D . Для «белого» шума, у которого М (t) = 0, а K (t1 t2 ) N0
2 t1 t2 , нормаль-
ная СВ будет иметь M = 0 и D = N0 |
2 |
T |
|
s2 (t)dt . |
|||
|
|
|
0 |
Часто приходится рассматривать неопределенный интеграл от стацио- |
|||
|
t |
|
|
нарного случайного процесса (t) |
( )d . Записанное выражение можно |
||
|
T |
|
|
|
- 94 - |
|
|
воспринимать как результат действия оператора Вольтерра с ядром, тождественно равным единице. Более общим случаем будет реакция линейной системы с импульсной характеристикой h(t, ) на стационарный случайный
|
t |
|
|
процесс (t), т. е. (t) |
h(t, ) ( )d . Среднее значение и корреляционная |
||
|
T |
|
|
функция данного СП будут равны соответственно |
|
||
|
|
t |
|
|
M (t) M (t) h(t, )d , |
|
|
|
|
T |
|
|
t1 |
t2 |
|
K (t1, t2 ) |
|
h(t1, 1)h(t2 , 2 )K ( 1 |
2 )d 1d 2 . |
T T
Как видно из приведенных выражений, процесс (t) нестационарный, так как среднее значение зависит от времени.
Если линейная система стационарна (имеет постоянные параметры), то h(t, ) = h(t ) и в установившемся режиме (Т = ) процесс (t) будет ста-
ционарным в широком смысле. Его математическое ожидание
M (t) M (t) h(t)d
0
не зависит от времени, а корреляционная функция K (t1, t2 ) будет зависеть от разности t1 t2 . Убедиться в этом мы предлагаем читателю. Для “белого” шума M (t) 0 , а КФ определяется выражением
|
t1 |
t2 |
|
2 ) N0 2 ( 1 2 )d 1d 2 = |
K (t1, t2 ) |
|
h(t1 |
1)h(t2 |
|
|
|
|
||
t2 |
|
|
|
|
N0 2 h(t1 |
2 )h(t2 2 )d 2 N0 2 h( x)h( x t1 t2 )dx |
|||
|
|
|
|
0 |
|
|
N0 2 h( x)h( x )dx . |
|
0 |
|
Дисперсия процесса (t) будет равна в этом случае |
|
|
|
D t = K (0) = N0 |
2 h2 (t) dt . |
|
0 |
- 95 - |
|
|
t |
|
|
|
Для процесса (t) |
(t) dt, t T |
среднее значение и дисперсия бу- |
||
|
T |
|
|
|
дут равны соответственно |
|
|
|
|
t |
|
|
t |
t |
M (t) M (t)dt, t T ; D (t) |
|
K (t1 t2 ) dt1dt2 . |
||
T |
|
|
T T |
|
Для «белого» шума M (t) 0 , а дисперсия D t = N0
2 t T .
В установившемся режиме для стационарной линейной системы при отыскании корреляционной функции на выходе можно воспользоваться частотным методом, заменив в выражении
|
t1 |
t2 |
|
|
K (t1, t2 ) |
|
h(t1 |
1) h(t2 |
2 ) K ( 1 2 ) d 1 d 2 |
корреляционную функцию K ( 1 2 ) ее представлением через спектраль-
|
|
|
|
ную плотность S ( ) : K ( 1 2 ) = 1 2 S ( ) exp j 1 2 |
d . |
||
|
|
|
|
Таким образом, |
|
|
|
t1 |
t2 |
|
|
K (t1, t2 ) |
h(t1 1) h(t2 2 ) 1 |
2 S ( ) exp j 1 2 d d 1d 2 |
|
|
|
|
|
Меняя порядок интегрирования по 1, |
2 и и учитывая, что |
|
|
|
|
|
|
h(t) exp j t dt K ( j ) ; |
h(t) exp j t dt K * ( j ), |
||
0 |
|
0 |
|
где K(j ) – коэффициент передачи линейной системы, получим окончательно:
|
|
K (t1, t2 ) 1 2 S ( ) K ( j ) |
2 exp j t1 t2 d , |
откуда видно, что K (t1, t2 ) действительно зависит от разности t1 t2 , а про-
изведение S ( ) K ( j ) 2 является спектральной плотностью процесса (t) :
S ( ) = S ( ) K ( j ) 2 .
Контрольные вопросы
1.Какие виды сходимости последовательностей СВ Вам известны?
2.Дайте определение непрерывности СП в точке и в области.
-96 -
3.Сформулируйте необходимые и достаточные условия дифференцируемости СП.
4.Какой вид имеют взаимная корреляционная функция и взаимный энергетический спектр процесса и его производной?
5.Как связаны между собой корреляционная функция стационарного
СП и его производной? Какова связь между энергетическими спектрами СП x t и x' t ?
T |
|
6. Запишите ПВ случайной величины s t x t dt , где |
x t – нор- |
0 |
|
мальный «белый» шум с S N0
2 , s t детерминированная функция.
7. На вход интегрирующей RC-цепи с постоянной времени T |
в момент |
|
времени |
t 0 подается сумма постоянного сигнала и «белого» |
шума с |
S N0 |
2 . Найдите среднее значение и дисперсию СП на выходе цепи. |
|
Будет ли этот процесс стационарным? |
|
|
8. Как связаны в установившемся режиме линейной цепи с коэффици- |
||
ентом передачи K j СПМ на входе и выходе цепи? |
|
|
Глава 7. НОРМАЛЬНЫЙ СЛУЧАЙНЫЙ ПРОЦЕСС
Следствием центральной предельной теоремы является чрезвычайно широкое распространение в радиотехнике и других областях науки и техники нормального случайного процесса. Если выполняется условие F fэ , гдеF – полоса пропускания линейной системы, а fэ – эффективная ширина спектра входного случайного процесса, то значения выходного СП в произвольный момент времени можно приближенно рассматривать как взвешенную сумму fэ
F независимых случайных величин. При этом независимо от распределения отсчетов входного СП предполагается, что условия центральной предельной теоремы выполняются. Тогда выходной процесс будет приближенно нормальным. Очевидно, что степень приближения зависит от числа независимых слагаемых fэ
F и распределения отсчетов входного процесса. Это явление называют эффектом нормализации.
Совокупность отсчетов нормального СП образует нормальный слу-
чайный вектор 1, 2 , , n т , статистические свойства которого были
подробно рассмотрены ранее. Следствием этих свойств является то обстоятельство, что нормальный случайный процесс полностью определяется ма-
- 97 -
тематическим ожиданием M (t) m1(t) |
|
|||
xW ( x; t) dx и корреляционной |
||||
|
|
|
|
|
функцией |
|
|
|
|
|
|
|
x1 m1(t1) x2 m1(t2 ) W ( x1, x2 ; t1, t2 ) dx1dx2 , |
|
K (t1, t2 ) M11(t1, t2 ) |
|
|
||
т. е. его полное описание дается в рамках корреляционной теории.
Если нормальный СП подвергается линейному преобразованию, то в результате мы получаем нормальный случайный процесс, если преобразование осуществляется с помощью оператора, либо нормальную случайную величину, если речь идет о линейном функционале от нормального СП.
Например, производная дифференцируемого стационарного нормального СП есть стационарный нормальный процесс с нулевым средним значе-
нием и корреляционной функцией K ( ) d 2K ( )
d 2 , где K ( ) – КФ исходного процесса.
Спектральная плотность процесса (t), как уже отмечалось, равна
S ( ) 2S ( ) , где S ( ) – спектральная плотность исходного процесса.
Результат действия оператора Вольтерра с ядром h(t, ) на стационар-
t
ный нормальный процесс даст нормальный СП (t) h(t, ) ( )d , который
t0
в общем случае уже не будет стационарным. Его среднее значение
t
M (t) M (t) h(t, )d
t0
и дисперсия
t t
M (t) M (t) 2 h(t, 1)h(t, 2 )K ( 1 2 )d 1d 2 .
t0t0
зависят от времени.
Если t0 = – , а h(t, ) = h(t ) , что соответствует установившемуся режиму отклика стационарной линейной системы с импульсной характеристикой h(t) , то процесс на выходе можно считать стационарным нормальным процес-
сом со средним значением M (t) M (t) h(t)dt const и дисперсией
0
- 98 -
2 |
t |
t |
|
|
|
h(t 1)h(t2 2 )K ( 1 2 )d 1d 2 |
|||
|
||||
|
|
|
||
|
|
|
|
|
|
h(t1)h(t2 )K (t1 t2 )dt1dt2 |
const , |
||
|
|
0 |
|
|
не зависящими от времени. Корреляционную функцию можно найти с помощью обратного преобразования Фурье спектральной плотности выходного
|
|
|
|
|
|
|
|
||
процесса S ( ) = S ( ) |
|
K ( j ) |
|
2 |
, где K ( j ) h(t) exp j t dt – коэффици- |
||||
|
|
||||||||
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
||
ент передачи линейной цепи, т. е. K ( ) 1 2 S ( ) |
|
K ( j ) |
|
2 exp j d . |
|||||
|
|
||||||||
Часто приходится иметь дело с линейным функционалом вида
T
(t) s(t) (t)dt ,
0
где s(t) – детерминированная функция, (t) – стационарный нормальный СП. Как результат линейного преобразования нормального СП, является
нормальной |
случайной |
величиной |
со |
средним |
значением |
|
T |
T T |
|
|
|
M M s(t)dt и дисперсией D |
s(t1)s(t2 )K (t1 t2 )dt1dt2 . Если |
||||
|
0 |
0 |
0 |
|
|
корреляционная функция процесса имеет вид K (t1 t2 ) N0 |
2 t1 t2 , то |
||||
T
D = N0
2 s2 (t)dt .
0
Как уже отмечалось, процесс с такой КФ называется “белым” шумом по аналогии с белым светом, имеющим постоянную во всем частотном диапазоне спектральную плотность мощности S ( ) = N0
2 .
Важной моделью СП с независимыми приращениями (см. классификацию СП) является винеровский процесс, для которого разность отсчетов(t + ) – (t) при любом имеет нормальное распределение с нулевым средним значением М[ (t + ) – (t)] = 0 и дисперсией М{[ (t + ) – (t)]2} = 2 .
Корреляционная функция Винеровского процесса K (t1, t2 ) 2 min(t1, t2 ) ,
ее график приведен на рис. 7.1.
- 99 -
Винеровский процесс нестационарен, его часто называют процессом, описывающим броуновское движение частицы под действием ударов моле-
кул жидкости. Производная винеровско- |
|
|
|
|
|
|
|||||||
го процесса является гауссовским (нор- |
|
K (t1, t2) |
|
|
|
||||||||
мальным) «белым» шумом. |
|
|
|
|
|
|
|
|
|
|
|||
Поскольку нормальный СП полно- |
|
|
|
|
|
|
|||||||
стью определяется заданием среднего |
|
|
|
|
|
|
|||||||
значения и КФ, которые могут быть вы- |
|
|
|
|
|
|
|||||||
числены с помощью двумерной ПВ, то |
|
|
|
|
|
|
|||||||
задание двумерной ПВ при произволь- |
t2 |
|
|
|
|
t1 |
|||||||
ном расположении отсчетов полностью |
|
|
|
|
|
t1 = t2 |
|||||||
|
|
|
|
|
|
||||||||
определяет нормальный СП. Для произ- |
|
|
Рис.7.1 |
||||||||||
вольной двумерной |
|
ПВ справедливо |
|
|
|
|
|
|
|||||
представление вида |
|
|
|
|
|
|
|
|
|
|
|
|
|
W |
x , x |
|
p (x ) p |
|
|
|
|
|
|
|
|
|
|
2 |
2 |
(x |
2 |
) C |
Q |
(x )Q |
(x |
2 |
) , |
||||
|
1 |
1 1 |
|
k 0 l 0 |
kl 1k |
1 2l |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
||
где Q1k (x1) и Q2l (x2 ) – ортонормальные многочлены, соответствующие ве-
совым функциям p1(x1) и p2 (x2 ) . Коэффициенты Стп находятся с помощью умножения обеих частей записанного равенства на Q1m (x1) Q2n (x2) и двойного интегрирования по области, определенной весовыми функциями, что с учетом ортонормальности Q1m (x1) и Q2n (x2) дает:
|
|
|
|
|
|
x , x |
|
Q |
|
|
|
|
|
|
C |
mn |
|
|
W |
2 |
( x ) Q |
2n |
( x |
2 |
) dx dx |
. |
|||
|
|
|
|
1 |
1m |
1 |
|
1 2 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Так как произведение полиномов Q1m (x1) Q2n (x2) есть выражение вида
m n
aib j x1i x2j , i 0 j 0
то коэффициенты Стп выражаются через корреляционные моменты mij. В роли весовых функций обычно выступают одномерные ПВ W 1 x1 и
W 2 x2 . Для стационарного нормального СП с нулевым средним значением и единичной дисперсией выполняются соотношения
|
|
|
x = |
|
1 |
|
|
x2 |
|
|
|
|
x |
= |
|
1 |
|
|
x2 |
|
p (x ) =W |
|
|
|
|
exp |
1 |
, p (x |
) = W |
|
|
|
|
exp |
2 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|||||||||||||
1 1 |
|
1 |
|
2 |
|
|
2 |
|
2 2 |
|
|
2 |
|
2 |
|
|
2 |
|||
|
|
1 |
|
|
|
|
|
2 |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
- 100 - |
|
|
|
|
|
|
|
|
|
|
|
и разложение двумерной ПВ имеет вид:
W |
x , x |
= |
|
|
1 |
|
|
|
|
x2 |
2rx x |
2 |
x2 |
|
|||||||||||
|
|
|
|
|
|
exp |
|
|
1 |
|
1 |
|
|
|
2 |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|||||||||||
|
|
1 |
|
2 |
2 1 |
r2 |
|
|
|
|
|
|
|
2(1 r |
) |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
1 |
|
|
|
x2 |
x2 |
|
|
|
r n |
|
|
|
|
|
|
|
|
|
|
||||
|
= |
|
|
exp |
|
1 |
2 |
|
|
|
|
|
H |
n |
(x )H |
n |
(x |
2 |
) , |
|
|||||
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
2 |
|
|
|
|
|
2 |
|
|
|
n! |
|
1 |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
n 0 |
|
|
|
|
|
|
|
|
|
|
||||||
где r – коэффициент корреляции отсчетов процесса, а Hn (x)
=


n! – нормиро-
ванные полиномы Эрмита, с которыми мы познакомились в гл.6 первой части учебного пособия.
Для общего случая двумерного нормального распределения выражение для двумерной ПВ примет вид:
|
|
|
|
1 |
|
|
2 |
|
2 |
|
|
W |
x , x |
|
= |
exp |
|
x1 M 1 |
|
x2 M 2 |
|
|
|
2 |
|
||||||||||
|
1 |
|
2 1 2 |
|
|
2 |
|
2 |
|
|
|
|
|
|
|
|
|
2 1 |
|
2 2 |
|
|
rn |
x M |
|
x |
2 |
M |
|
||||||
|
|
|
1 |
1 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|||||
n! |
Hn |
|
1 |
|
Hn |
|
2 |
|
. |
|||
n 0 |
|
|
|
|
|
|
|
|
||||
Узкополосный нормальный случайный процесс. Корреляционная функция и связанная с ней спектральная плотность СП были рассмотрены ранее. Для узкополосных СП, как и для детерминированных сигналов, весьма
продуктивным является представление исходного СП (t) в |
виде |
|
(t) a(t) cos 0t b(t) sin 0t , где процессы |
a(t) (t) cos (t) |
и |
b(t) (t) sin (t) называют квадратурными компонентами процесса (t). Переход от процесса (t) к огибающей и фазе осуществляется с по-
мощью сопряженного процесса (t), получаемого из исходного с помощью оператора (преобразования) Гильберта:
1 ( )
(t) = Н (t) = t d .
Для существования процесса (t) достаточно потребовать равенства нулю среднего значения процесса (t), т. е. М (t) = 0 при любых t.
С помощью процесса (t) можно построить аналитический СП
(t) = (t) + j (t) = (t) exp j (t) exp j 0t ,
где комплексная случайная функция (t) exp j (t) называется комплексной огибающей СП (t).
- 101 -
